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IlfTRODUCTION 

Introduction 

An algebraic language translator-compiler is a computer program 
which reads and translates mathematical algorithms written in a language 
close to mathematical notation, and produces ^ i-e.^ compiles^ a computer 
program which when executed will perform the described algorithms . 

MAD_, the Michigan Algorithm Decoder, is an algebraic language 
translator=compilerj programmed at the Computing Center of the University 
of Michigan. This manual is a description of the MAD language, the 
language in which algorithms must be written In order to allow MAD to 
translate them. It is also a guide to the use of MAD. 

This mahual was designed and organized as a reference manual. It 
is not intended by itself to meet the requirements of a text book. The 
book, THE LANGUAGE OF COMPILERS; AN IHTRODUCTION , by Bernard A. Galler, 
to be published by McGraw Hill^ does, very satisfactorily, meet the re- 
quirements of a text book. The combined use of this manual and that text 
book should provide the user with a satisfactory introduction to MAD- 

This goal of producing an easily used reference manual rather than 
a text has dictated several of the noticeable features of this book. Wo 
index is provided, but a complete Table of Contents labeled by a nested 
decimal classification system is used. Ho page numbers are used in the body 
of the work, but the material is labeled in parallel with the Table of Contents . 
Gaps provided in the Table of Contents and the use of loose-leaves will 
enable the manual to be conveniently revised and augmented . 

In keeping with the reference nature of the manual, an attempt 
at precise definition of the terms and elements of the MAD language has 
been made. even at the expense of often being quite formal. 
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1 A Preliminary Survey 

1.1 The Hollerith Character Set 

The characters to "be used in writing in the MAD language are as 
follows : 

(1) Alphabetic Characters^ or Upper Case Letters: 

A, B, C, D, E, F, G, H, I, J, K, L, K, N, 0, P, Q, R, S, T, 
U, V, W, X;, Y, Zj 

(2) Wumeric Characters, or Decimal Digits: 
0, 1, 2, 3s ^, 5, 6, 7, Q, 91 

(3) Special Characters (The characters ", quotation marks, are 
not characters of the language. They are used here to set off 
the characters): "+", "-" (minus), '"" (apostrophe), "=", "*" , 

tl/ll IT/l! llNri II II II II /____„ \ ll(tll II II /VT„_1, „_„„^ 

/, \ , ) , • , , ^corama;, $ , (.Dlank space, 
sometimes represented by the character "b"). 
Groups (1) and (2) are referred to collectively as Alphanumeric 
Characters . 
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1.2 The Sequence of Events in Using MAD 

(1) The problem is expressed in the MAD language as described in 
this manual. This expression of the problem is a "source 
language program" . 

(2) The source language program is ingested by the MA.D translator- 
compiler (see section 14^ Mechanics of Using MAD). The MA.D 
translator reads and translates the source language program 
and from it compiles a "machine language" computer program. 
This computer program is the "object program". 

This phase is called "compilation"". The time at which the 
program is compiled is said to be "compilation time". 
(5) The object program is executed on the computer. This is the 
"execution phase" and the time at which the object program is 
executed is said to be "execution time". 
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1.3 Statements 

A source language program is an expression of a problem in the form 
of algorithms which, when executed, provide the desired solution. 

Algorithms are expressed in MAD by writing a series of "statements". 
Statements may be either executable statements or non-executable state- 
ments, known also as "declarations". It is possible further to sub- 
divide statements into classes such as "control statements", "testing 
statements", "iteration statements", "input-output statements, etc. 
All statements which are available in the MAD langioage are explicitly 
defined in sections 5> ^, 5j 6 and "J. 

Every statement must begin on a card, column 11 of which is blank. 
There may not be a remnant (a tail end) of a preceding statement on the 
same card. (This information appears elsewhere in this manual j it is 
given here, also, for emphasis to aid readers in overcoming preconceived 
notions due to familiarity with other translators.) Any card is recognized 
as the last card comprising a statement when the card which follows it has 
a blank in column 11. 
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l.h MAD and the Operating System 

All of the properties which are strictly properties of MAD and which 
are necessary to enable a user to write a MAD program are described in 
this manual. There are certain other facilities which might be desired 
andj which^ although not described herein^ may exist. Generally, if such 
options exist^ they exist as properties of the operating system of which 
MAD is a part^ and the user must consult the description of the properties 
of the pertinent operating system in order to determine what other useful 
options are available. 

By the same token, it is true in most computing installations that, 
in order to run a MAD program to solve a problem, the existence of a MAD 
source program designed to do just that is not sufficient. Again, the 
operating system of which MAD is a part must be considered. The very 
properties of the operating system which make it flexible and useful 
impose requirements on a problem in addition to those necessary to satisfy 
MAD. 

The intention of any operating system is to provide greatly increased 
flexibility and many useful options with but few accom^panying extra re- 
quirements . 

In particular it is in the province of MAD to translate the source 
language program and compile an object program. After this, MAD has no 
more control over the destiny of the particular problem. It is in the 
province of the operating system to load the MAD object program into the 
computer, and transfer control to it so that it may execute. After 
execution of the object program is completed, it is the operating system 
which assumes control of the computer and proceeds to the next problem. 

Examples of properties which may be available to users of MAD l^ut 
which are in reality properties cf the containing operating system are= 

(l) Compilation with the option of executing the object pro- 
gram on the same computer run; 
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(2) Execution of a problem which requires more than one "core 
load "j each of which is an object program resulting from 
separate compilations . This option is variously known 

as "ping-pong" or "chaining". 

(3) Translation and compilation processing of more than one 
MAD source program in a single computer run; 

(k) Use of an assembly language translator and the MAD compiler 
together on the same problem. 
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2 Elements of the Language 

In this section the basic building blocks of the language are de- 
sc3?ibed . These elements are used in the construction of the statments 
described in chapters 5, ^, '?, 6 and 7- 

A brief discussion of the concept of "mode" is in order. It is a 
basic concept of the language. The "mode of an element" (and specific 
mention will be made throughout the remainder of the manual regarding the 
specific elements to which the concept is applicable ), is determined by 
the qiiantities which the element in question may represent. The mode of 
any element in turn influences the manner in which operations are performed 
on the element. 

There are five and only five modes in the MAD language. They are 
the following: 

Integer 

Floating Point 

Boolean 

Statement Label 

Function Name 

(a) Any element of integer mode will be interpreted as representing 
a value eqoal to an integer i such that i satisfies 

-9999999999 < i < +9999999999 
The decimal point is assumed to be immediately to the right of 
the rightmost digit and no fractional part may appear. 

(b) Any element of floating point mode will be interpreted as 
representing a value equal to any real number r such that r 
satisfies : 

10-58 < j^i ^^q38 

Thus such an element is a "mixed number". 
Example to cover (a) and (b); 

Although the integer 1 is a real number 1 and the real number 

2 is an integer, the integers 1 and 2 are quite distinct from the 

floating point . pumbers . 1' and, 2 , respectively . 
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In particular J using octal digits (O through 7) to re- 
present groups of three binary bits (O^ l) the internal 
representations of 1 and 2 are as follows in actual internal 
IBM 7090 representation. 

Number Floating Point Integer 

1 20lif00000000 000000000001 

2 202400000000 000000000002 

(c) Any element of Boolean mode will be interpreted as representing 
a value equal to one of: 

"true", 
"false". 

These values will be represented in MAD in the following way: 

"true" by IB^ 
"false" by OB. 

(d) Any element of statement label mode will be Interpreted as 
representing a value equal to a statement label. See section 
2. 3 J Statement Labels, for a complete description. 

(e) Any element of function name mode will be interpreted as re- 
presenting a value equal to a function name. See section 2.h, 
Functions, for a complete description. 

Although, abstractly, the value of an element of any given mode may 
be "equal" to the value of an element of one of the other four modes, the 
format of the representations internal to the computer is different and 
operations on the elements will be performed differently. 
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r\ -[ n ^^^ — 1-« 

£:.j. uuiiHoaiibS 

The following types of constants may be written in the MAD language. 

2.1.1 Integer Constants 

Integer constants may be one to ten decimal digits preceded by 
sign (- or +). Hence, integers can range from -9999999999 to +9999999999 
with the decimal point always assimied to be immediately to the right of 
the rightmost digit but always omitted . Integer constants are of integer 
mode. 

Signs 

(1) Sign if present is always placed to the left of the integer. 

(2) Negative sign is never omitted. 

(3) Positive sign may be (and usually is) omitted. 
E3CAMPLE': 2, -2, 0, +0, 100 are all integers. 

Leading Zeros 

Leading (but not trailing) zeros may be omitted. 
EXAMPLE? 5 and OO5 are the same but 3 and 300 differ. 

2.1.2 Floating Point Constants 

There are two basic forms, with and without exponent. In all cases 
a floating point constant F must satisfy: 

-zQ -zQ 

10"^" < |f| < 10^" 

Floating point constants are of floating point mode. 

(a) Without exponent the constant contains one to eight digits and 
a decimal point ( . ) which must be written but which may appear 
anywhere in the number . Examples ; 

0., 1,5, -0.05, +100-0, -k. 

(b) With exponent the constant contains from one to eight digits 
with or without a decimal point y followed by the letter E, 
followed by the exponent. Exponent is one or two digits pre- 
ceded by sign and represents the power of ten by which the 
number to the left of the exponent is to be multiplied. 
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Examples ; 




MAD 


NOTATION 




.05E-2 


- 


.05E2 




5E02 


5 


.E2 



MEANING 

,05 X 10'^ 
-.05 X 10^ 
5 X 10^ 
5 X 10^ 

Note thatt 

(1) Plus signs may be optionally omitted in front of the exponent 
as well as in front of the number itself. 

(2) Decimal point may be omitted^ in which case it is asstmied to be 
immediately to the left of the letter E. 

O 1 X ft 1 -y-»ViciV»£a-f--? r» r^i^-no-f-Qn-f-o 

Alphabetic constants consist of from one to six characters from the 
following set of admissible characters; 

"A", "B", ...J "Z": alphabetic characters. 

"0", "1", ..., "9": decimal digits „ 

"+", "-" (minus), •"" (apostrophe), "=", "*% 

"/") "('% ")"} "'"} ">" (comma)^ " " (blank space, sometimes represented 

by the character %"). 

Each alphabetic constant must be delimited on the left and on the 
right by the special character "$". The mode of alphabetical constants 
used in arithmetic expressions is integer. 

Notice especially that although blank spaces are ignored elsewhere in 
the language, except where it is specifically indicated that they are not, 
they count as characters in alphabetic constants . 

EXAMPLES : $ABCr$ , $T0 BE$ , $:EC . h^ ,$5 +5=8$ . 

NOTE: An alphabetic constant is stored internally as an integer. An^'" 

alphabetic constant which is written with fewer than six characters 

will be stored internally left ,1ustified with blanks appended on 

the rights thus,$ABCD$ will be stored internally as if it had 

been written $ABCDbb$ . 
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Since the character "$" serves to delrnit an alphabetic 
constant^ it is not possible to write an alphabetic constant 
with "$" as one of the internal characters . An alphabetic 
constant with the character "$" as the left-most of the 
internal characters may be created, however, by writing the 
following alphabetic constant: 

2.1A Boolean Constants 
There are two Boolean constants: 

IB is written for "true", 
OB is written for "false". 

2.1.5 Octal Constants 

Octal constants may be written in two ways, with or without scale 
factor. The mode of octal constants is integer. 

Without - scale factor an octal constant is written with from one to 
twelve octal digits followed by the letter "K". 

EXAMPLE : OJK, 12 5K, 777777777777K 

With scale factor one or more decimal digits follows the "K". This 
decimal integer is interpreted as the exponent of the power of eight by 

n i.J...t^^J,J. UAA^w V.'V- ^^l-^U_ J-O-Jt l-"^-.^ v.. J. j^J. V^V.- V^\<t. J_AA^ VAa\— J.L _l_l-' Vl»/ t-'V* UJ«X_1- l-'^_^.X.^V_\A • 

EXAMPLE: 127K2 becomes OOOOOOOI27OO, 
IKLO becomes 010000000000- 

2.1.6 Statement Label Constants 
See section 2.3, Statement Labels. 

2.1.7 Function Name Constants 
See section 2.^, Functions. 

2.1.8 Redefinition of The Mode of Constants 

Any constant whether Integer, Floating Point, Boolean, Alphabetic 
or Octal, may be declared to be of a mode other than its normal mode. 
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When a constant which is to have its mode redefined is used it is written 
in the form 

cMi 

where "c" is the constant in q.uestion written as described above for each 
of the five types, "M" is the character "M" written following the last 
character of "c" to indicate mode redefinition and i is one of the integers 

for Floating Point 

1 for Integer 

2 for Boolean 

3 for Function Name 

k for Statement Label 

The conversion performed on the constant in reading it and the form 
in which it appears internally is that associated with Its original mode. 
The mode of the constant is reassigned to be the mode whose code is i 
after the standard processing of the constant is completed. 
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2 .2 Variables 

The name of any variable consists of from one to six alphabetic 
or numeric characters j the first of which must be an alphabetic character. 
The name used for any variable must not be identical to any statement 
label or to the alphanumeric part of any function name. 

EXAMPLES : tEMP,KLMZjF55,P32K,RESULT,X,ALPHA. 

2.2.1 Simple Variables 

A simple variable is a variable whose name refers to a unique 
element. It may be thought^of - as a degenerate vector with one element,, or as 
a matrix of order (l,l). 

Simple variables are written in the form 

mm 

where NAME adheres to the form of a variable name as defined above. 'Each 
of the names in the immediately preceding example is written as a simple 
variable . 

2 .2 .2 Arrays 

An array is a set of elements. The number of elements in the set 
is greater than one. Every variable name which is the name of an array 
must be so declared as explained in section 3.6, DIMENSIjz5n Declarations. 

Arrays of order (m^l) or (l,n) where ra and.n. > 1/ are "Vectors", "Linear 
Arrays", or "One Dimensional Arrays". Vectors are stored internally with 
the elements in adjacent locations . 

Two dimensional arrays, or matrices, of order (m,n) are stored 
internally by rows, i.e., in the order: 

^11^^12' • "'^ln-'^21-'^22^ • • •'^2n' * * *•' V'^m2^ ' ''^mn " 

In general, n-dimensional arrays, where n > 2, are stored in the order de- 
termined by varying the rightmost subscript first, then the next rightmost, 
etc. Thus, a four dimensional array of order (m,n,p,q.) would be stored in 
the order: 
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^111^^1112-' " * "'^111^^^1121^^1122^ • • •'^112q' * ' •'^llpq^^l211, 

' Inpq.^ ■* mnpq 

It is possible in MAD to regard any array of dimension n, even 
though n may he greater than one, as a linear array. The order of the 
elements of an array when regarded as a linear array is as described in 
the preceding paragraph. Any element of an array may thus be referred 
to by a subscripted variable with a subscript of one element rather than 
n. This single element subscript of an array element is called a "linear 
subscript" of the array element. The value of the linear subscript is 
equal to the position of the given element relative to the element with 
all subscripts zero^ the linear subscript of which is also zero. 

See sections 2.2.k, Subscripts of Subscripted Variables, and 3 •6. 2, 
Ifetrix DIMENSI'^N Declaration- for additional comments concerning linear 
subscripts . 

2.2.5 Subscripted Variables 
A subscripted variable is written in the form 

NAME (S^,S2,...,S^) 

where MME adheres to the form of a variable name as defined above, and 
must be the name of an array variable, and each S , i = 1, 2, ..., n, is 
a permissible subscript as described in sections 2.2.4, Subscripts of 
Subscripted Variables, and 2.2.9, Statement Label Variables. The characters 
"(" and ")" delimit the subscript and the character ", " separates the 
elements of the subscript. Notice that the subscripts written in the MAD 
language are written on the same line as the name they are subscripting 
and are not written below the line as in usual mathematical notation. 
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AX sny given ulsus s suuscnpoSu. vsriauxe is 'tiie nsnie of a single 
specific element of an n-dimensional array of elements . 

At any given time the value of a subscripted variable is a single 
specific element of an n-dimensional array of elements. 

It is a convention of the language that if a name^ NAMEA^ which is 

uj.±v^ Aii^iuw wj- f-i.j.± Qx J. a v« _Lo wj. iULicrii cij_i_jiir; _i_xi ilii, CJf 'I."GSJi J.yXJ. ciis J.1 IX* "V-tj^^t^o 

a simple variable, i.e., 

MMEA 

its value is equivalent to the value of the subscripted variable 

NAMEA(0,0, ...,0). 

Examples of subscripted variables whose values are array elements: 

BETA (I ) 

X(J,5) 

Y(7) 

N(B+J4-*F,Q) 
M(k(2+5,)+T.(1)+6) 
BETA = BETA(O) 

x~ x(o,o) 

2.2.4 Subscripts of Subscripted Variables 

With one exception, the following discussion applies to all types 
of subscripted variables. See section 2.2.9, Statement Label Variables, 
for the exception. 

As far as the translator is concerned, any arithmetic expression 
(see section 2.6, Arithmetic Expressions) may be used as a subscript when 
the subscript is a linear subscript. This is true whether the dimension 
of the array is one or greater than one. 

The qualification "as far as the translator is concerned" arises 
from the fact that the value of the expression must be greater than or 
equal to zero. This restriction applies when the expression is evaluated 
at execution timej violations are not discernible by the translator at 
compilation time. 
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If the value of the subscript expression used as a linear subscript 
is in floating point mode (see section 2.6.k, Mode and Arithmetic Ex- 
pressions )j it will be truncated to an integer before it is used as a 
subscript . 

When, in the general form of a subscripted variable, the dimension 
n > 2, each subscript element may be an arithmetic expression but the 
expression must be of integer mode. Moreover, the use of subscript ex- 
pressions having values of other than integer mode will go undetected 
by the translator when n > 5- Unlike the case of linear subscripts, it 
may be meaningful for a subscript expression appearing as an element in a 
non-linear subscript, to have a value less than zero- Care must be exer- 
cised. See section 3.6.2, Matrix DIMENSION Declaration, for a complete 
discussion. 

2.2.5 Integer Variables 

All of the preceding general discussion of variables applies to the 
integer variables. Values represented: any integer number I satisfying 

-9999999999 < I < +9999999999 
These variables are of integer mode. 

2.2.6 Floating Point Variables 

All the preceding general discussion of variables applies to the 
floating point variables. Values represented; any real number F satisfying 

10-5S< |F| <105^. 

These variables are of floating point mode. 

2.2.7 Boolean Variables 

The values represented by Boolean variables are 

IB = true, 
OB = false. 

These variables are of Boolean mode. 
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The discussions under the sections Simple Variables, Array Variables, 
Subscripted Variables and Subscripts of Subscripted Variables apply to 
Boolean variables . 

2.2.8 Function Name Variables 

The values represented by function name variables are any function 
name, see section 2.4, Functions. These variables are of fimction name 
mode. 

The discussions under the sections Simple Variables, Array Variables, 
Subscripted Variables and Subscripts of Subscripted Variables apply to 
function name variables • 

2.2.9 Statement Label Variables 

The values represented by statement label variables are any statement 
labels, see section 2.5, Statement Labels. These variables are of state- 
ment label mode. 

The discussion of Simple Variables applies to statement label 
variables . Subscripted statement label variables are more restricted 
than general subscripted variables . Because of the way in which elements 
of statement label vectors are defined, "arrays" of statement label mode 
are restricted to one dimensional arrays. Thus, subscripts of subscripted 
statement label variables must satisfy n = 1, of the general form of a 
subscripted variable described in section 2.2-3, Subscripted Variables. 
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2.3 Statement Labels 

Statements are defined in sections 3 and 4. A statement must be 
labeled when it is desired to have another statement refer to it . 
A statement label may take on one of two forms: 

(a) C 

where C consists of from one to six alphabetic characters or 
decimal digits j the first of which must be alphabetic; 

(b) C(n) 

where C is as in (a), n is an integer constant and the special 
characters "(" and ")" delimit n on the left and right re- 
spectively. A statement label of this form is an element 
of a statement label vector. 
A statement label is a constant of statement label mode. 

to the alphanumeric part of any function name. 

A statement label appears in the label fields columns 1 through 
10, of the statement it identifies . When a statement extends to additional 
cards (see section 14.1.2.2, Continiiation Cards), the statement label 
need not be piinched on the additional cards • 
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2.k- Functions 

The name of a function is written in the form 



where C may consist of from one to six alphabetic characters or decimal 
digits, the first of which must be alphabetic. The special character 
" . " is part of the name and must follow the last character of C . A 
function name is a constant of function name mode. The alphanumeric 
portion of the name must not be identical to any variable used in the 
program, or to any statement label. 
Examples of function names t 

SIN. 

SORT. 

ELjfe. 

2.i)-.l Single Valued Functions 
The value of a single valued function is represented, in general, 
by the forms 

In the form, "C." is a function name and each A., i = l,2,...,n, is an 
argument of the function (see section 2.k.^, Arguments of Functions). 
The special characters "(" and ")" delimit the field of arguwentfe-and ' 
the special character "," separates pairs of arguments. The single valued 
functions may be used by specifying the appropriate value form as an operand 
in any arithmetic expression (see section 2.6, Arithmetic Expressions). 

The value form of a function is also often referred to as the "call" 
for the function, i.e., the, value form "calls" for the function to be 
evaluated . 

Examples of .values of single valued functions? 

BIN. (3 -1^1396) 

SQRT.(X+Y*Z-10.) 

MLTRGF. (A,B-C*D,2 .*X,5) 
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2 A. 2 Non-Single Valued Functions 
Since the result of the operation of a single valued function 
upon its arguments is by definition a single value, it made sense in 
section 2A.1 to define the value of a single valued function. In MAD 
the facility exists for allowing arguments to be operated upon by 
"functions" which may not be used in any expression. These are called 
non-single valued functions . The notation to indicate such an operation 
is identical to that used for the value of a single valued function, namely. 



C ■ V-^^-l }"-0} ' ' ' f-^yy ' 



where the elements of the form are exactly as described in section 2.4.1, 
Single Valued Functions. The only distinction, which is, however, a 
non-trival one, is that the name used for any non-single valued function 
must be different from the names of all single valued functions and 
conversely. 

The additional mechanism required to execute non-single valued 
functions is explained in section 4. 8 EXECUTE Statement. 

Example of the specification of a non-single valued function which 
is to sort the n elements of some list according to prescribed rules: 

sjz5rt.(list,w) 

2.4.3 Translator Defined Functions 

Certain commonly used functions are already available to the MAD 
translator. The programmer may use these functions merely by writing 
them in the appropriate manner. A list of the ciirrently available 
translator defined functions may be found in the manual describing library 
subroutines . 

2.4.4 Programmer Defined Functions 

The facility exists in MAD to permit the programmer to define- for 

his use functions of any complexity. Two types of functions, "internal 

functions" and "external functions", may be defined. The proced\ire for 

defining functions is explained in section 5-8 Single Statement INTERNAL 

FUNCTION Definition, and in Section 6, Function Definitions. 
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The form to indicate the use of any function of either of these 
types is the same as that described in section 2.i4-.l, Single Valued Functions. 
Programmer defined functions may be either single valued or non-single 
valued functions . 

2.4.5 Arguments of Functions 

In general, as far as the MAD translator is concerned, any argument 
of any function may be any expression. For a given function however, 
the function itself imposes restrictions on the arguments. Therefore, 
in determining the restrictions upon the arguments for a given function, 
it is necessary to consult the function description. 

2.4.6 The Mode of a Single Valued Function 

By "the mode of a fvmction" is meant "the mode of the value of a 
function", e.g., if the mode of the value of a function is integer, then 
the mode of the function is integer. 

The mode or modes of various argiunents of a function do not prejudice 
the mode of the function. There may be functions all of whose arguments 
may have modes different from the mode of the value of the fxmction. 
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2.5 Arithmetic Operations 

A number of unary and binary arithmetic operations are available. 
The arithmetic of the binary operations is integer arithmetic or floating 
point arithmetic according to the modes of the operands • If both binary 
operands are integer then integer arithmetic is used. If both operands 
are floating point, floating point arithmetic is used. If the operands 
are of mixed modes, the operand in the integer mode is first converted 
to floating point; the arithmetic operation is performed as a floating 
point operation, yielding a floating point result . 

2.5-1 Unary Operations 
These operate on the single operand immediately to the right of 
the operator. 



Cfperation Symbol 


Definition 


Example 


.ABS. 

+ 


Absolute 

Value 
Identity 

Negation 


.ABS.(X+Y) means (X+Y) 

+ (x/y) means (x/y) 
-(Q*P) means -Q*P 



2.5 -2 Binary Operations 
These operations involve the operand to the left and the operand 
to the right of the operation. 
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Oper. 
S3niibol 


Definition 


Example 


+ 


Addition 


Zi|-+QU0 


- 


Subtraction 


A-XYZ 


* 


Multiplication. Juxtaposition 
may not he used to signify 
multiplication . 


2.5*1 


/ 


Division. If both operands are 
integers the result is made an 
integer. The fractional part of 
the true quotient is truncated 
(not rounded). 


D5B/(C-5) 


.P. 


Exponentiation. Raise the left 
operand to the power which is 
the value of the right operand. 


VAR.P.2 
means 

(VAE^ 
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2.6 Arithmetic Expressions 

2.6.1 Def4.nition of Arithmetic Expressions 

(a) Every unsigned constant^ whether integer, floating point 

or alphabetic, is an arithmetic expression. Every individual 
variable, whether integer or floating point. Is an arithmetic 
expression. Every subscripted array variable, whether integer 
or floating point, is an arithmetic expression. Every value 
of a single valued function, whether integer or floating 
point, is an arithmetic expression. 

(b) If V is any arithmetic expression of a type in (a) above 
then both 

+ V 
- V 
are arithmetic expressions . 

(c) If E is any arithmetic expression, then 

(E) 

is an arithmetic expression, and also both 

+ (E) 
- (E) 

are arithmetic expressions . 

(d) If F is any arithmetic expression, then all of the following 
are arithmetic expressions: 

.ABS. F 

F + V 

F + (e) 

F - V 

F - (e) 

F * V 

F * (e) 

F / V 

F / (E) 

F .P. E 
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(e) The only arithmetic expressions are those arising in (a) 
through (d) which do not exceed a length which can he 
contained on ten cards . 

2.6.2 Hierarchy of Operations in Arithmetic Expressions 
The order in which individual terms of an expression are to be 
evaluated and combined must be unique and in the MAD language is made 
so by defining the hierarchy of the arithmetic operations . Unless 
altered by parentheses the order of arithmetic operations performed 
within an arithmetic expression, is given by the following list: 

Operation 
Symbol 

.ABS. , + (unary operations) (eqiial hierarchy) 

.P. 

(unary operation) 

* , / (equal hierarchy) 

+ , - (binary operations) (equal hierarchy) 

Within an expression, operations of equal hierarchy are performed 
from left to right unless otherwise indicated by parentheses . 
Example :1: The expression 

A + b/C + D .P. E*F - G 

means 

A + l+D^xF-G 

Example 2: A * B/C * D/E * F 
means 

J± x F 

E 
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Example J: C(k) + A(3) * B(j)/9.7 + 5-5 * P 

means 

Example ki A+B-C+D-E 
means 
(((A + B) - C) + D) - E 

Example 5? X/Z * y/R * S 
means 
(((x/z)* Y) /r) * S 

r^ C -z Ti -v-J-T -1" A — r 4-1 1-4 „ m--^ ; 

d^yj.j r£ii'cni-.nt;tjt:fc> xli au rtx x uxjiucoj.u jlajjx ec)tiiUJ:i 

Parentheses as used in the usual algebraic sense may be used to 
override the usual rules of precedence for a given expression and are 
frequently used for the purpose of sim,plifying mathematic expressions . 
Example 1: A * B = C/D + E * F 

means 

n 
AB " ^ + EF 

but 

A * (B - C)/(D + E) * F 

means 

A(B " C) 
(D + Ej ^ 

Example 2s (C(k) + A(5) * B(j))/(9.7 + 3-5 * P) 

means 



C(K) + A(.^) * Bfj) 
9.7 + 5.5 * P 
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2.6.k Mode and Arithmetic Expressions 

As just defined, it is possible for an arithmetic expression to 
contain operands of either integer mode or floating point mode. When 
an arithmetic expression is evaluated, the result, which is the "value 
of the arithmetic expression", must also be either of integer or floating 
point mode. Thus it is proper to speak of the mode of the value of an 
arithmetic expression^ a somewhat looser terminology, "the mode of an 
arithmetic expression", is sometimes -used in place of this. 

In an arithmetic expression an alphabetic constant is a constant 
of integer mode. 

An arithmetic expression is considered to be in the floating point 
mode if any operand of any arithmetic operation in the expression is in 
the floating point mode. If all operands are integer (or alphabetic ), 
then the expression is considered to be in the integer mode. In this 
determination arguments, though not values, of functions are ignored. 

Thus, if Y, Z, and W are floating point variables, while the 
function GCD. and the variables I and J are in the integer mode, then 
the expressions 

Y + GCD.(I,J) 

Y + Z - I 
I + 1. 
GCD.(I,J)/Z 

are all floating point expressions while the expressions 

I + GCD.(I,J) 

(I + J)/3 
I + 1 
GCD.(l,j)/l 

are all integer expressions . 
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If an arithmetic expression has subexpressions of different modes, 
a conversion may be necessary before some of the operations can be per- 
formed. Thus, in the expression 

Y + 5 

if Y is in the floating point mode it cannot be added directly to the 
integer 5' But for one precaution the user need not be concerned with 
this since the instructions necessary for the conversion of the integer 
to floating point form before adding are automatically inserted by the 
translator during the translation process . The precaution is that if 
the integer being converted is greater than 13^,217,726 (i.e., 2 ) 
then an improper conversion will take place . 

In some cases, however, the user must understand the sequence in 
which the conversions will be made. Consider the expression 

(Y + 7/5) + (I * J/k) 

where Y is in the floating point mode, and I, J, and K are in the integer 
mode. According to the parenthesizing conventions the computation will 
proceed in the following order (where the T's are temporary locations): 

T^ = I * J 

T^ = T^/K 

T^ = 7/3 
T^ = Y + T^ 
T^ = T^ + Tg 

and T will be the value of the expression. 

Now, since both I and J are integers, the first multiplication 
will be integer multiplication, and T, will be an integer. Since 
the following involves two integers, it will be integer division, and, 
if K happens to have a larger value than T, the quotient is 0. Similarly, 
T^ will have the value 2 becauce cf the division cf two integers. In the 
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computation of T^, however, we have "mixed modes," since Y is floating 
point and T is integer. Here T will be automatically converted to 
floating point before adding. Likewise, in the next step, the integer 
T„ will be converted to floating point before adding to the floating 
point number T^. 

In other words, although the mode of the expression is floating point 
because of the presence of the floating point variable Y, some of the 
computation (imtil Y is involved) is performed in integer arithmetic, and 
this may occasionally cause the final value to be different from the 
value one might expect from a different analysis . 

In the example above, the divisions would be performed in the floating 
point mode if the expression were written: 

(Y + 7-/5) + (I * J)/(K+ 0.) 

Of course, many times the expression will be written as originally 
stated just to achieve the "truncation" effect. 
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2.7 Mathematical Relations 

In order to permit comparison of the algebraic values of pairs of 
arithmetic expressions (see section 2.6, Arithmetic Expressions)^ state- 
ment label expressions (see section 2.10, Statement Label Expressions) 
or fimction name expressions (see section 2.11, Function Kame Expressions), 
the following mathematical relations are provided; 



MAD Symbol 


Mathematical Symbol 


Meaning 


.X. 


< 


E.L.F means "E is less than 


.LE. 


< 


E.LE.F means "E is less than 






or equal to F" 


,E. 


= 


E.E.F means "E is equal to F" 


.m. 


^ 


E.NE.F means "E is not equal 
to F" 


■.G. 


> 


E.G.F means "E is greater 
than F" 


.GE. 


> 


E.GE.F means "E is greater 






than or equal to F" 
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2.8 Boolean Operations 

The following Boolean, or logical, operations are available. Let 
M and P "be Boolean Expressions as defined in section 2-9, Boolean 
Expressions . 



MM) 


Corresponding 




Operation 


Logical 




Symbol 


Symbol 


Definition 


.nj6t. 


^ 


The value of .'N0I M is IB if and only 
if the value of M is OB. 


.OR. 


V 


The value of M.^E.F is OB if and only 
if both M and P have the value OB. 


.EX^R. 


& 


The value of M.EX^R.P is IB if and 
only if exactly one of M or P has 
the value IB. 


.AND. 


A 


The value of M.AKD.P is IB if and 
only if both M and P have the value 
IB. 


.TTTKN. 


z? 


The value of M.TFRW.P is OB if and 
only if the value of M is IB and 
the value of P is OB. 


.EQV. 


_^ 


The value of M.JIQU.P is IB if and 






only if the value of M is equal to 
the value of P. 
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2 .9 Boolean Expressions 

2.9.1 Definition of Boolean Expressions 

(a) Boolean constants. Boolean simple variables. Boolean subscripted 
array variables and Boolean-valued single valued functions are 
Boolean expressions (see section 2.1.^, Boolean Constants, 2.2.7, 
Boolean Variables, and 5"2, Mode Declaration). 

(b) If E and F are arithmetic expressions, then the following are 
Boolean expressions: 

E .L. F 
E .LE. F 
E .E. F 
E .GE. F 
E .G. F 

(c) If G and H are both statement lable expressions or both function 
name expressions, then 

G .E. H 
G .NE. H 

are Boolean expressions . 

EXCEPTION: If G or H are elements of a vector preset by a VECTOR 
VALUES statement (see section 5.7, VECT0R VALUES De- 
claration), then G .E. H and G -NE. H are not expressions > 

(d) If M and P are Boolean expressions, then the following are Boolean 
expressions: 

(M) 

.N0T. M 
M 4^. P 
M .EXJZ^R. P 
M .AND. P 
M .THEN. P 
M .EQV. P 
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(e) The only Boolean expressions are those arising from (a) through 

(d). 
Examples of Boolean expressions: 

(X .G. 5 .AMD. Y ,LE.2).j6r.(GAMMA.L.EKIM) 

/ J, ^rt ft-r-. -.^^ \ /-.r-. -r-ri ^-r^^-^^ ^y.\ . ^-rr^ f^ /,,-. \ ^ „., , ^-^ „ \ 

I £A K>% .ixi_x'^i/xi i.h: h:)-'>s i i jm i ti i«i i I H* I y I i i . m*i-'s< i : mi 

y •^ikX.^ti^ * ^^a_ju ^<^t_ y / ^^.^ • J.J1 M ■ jujj. *»f ^.i.iiAi / «f 1.X1J-' " \.'- * ^^^^ / « J-i ■ ■ ■■ ■^ ' ■ ' "■ y 

((P -AM). M) .THEN.M) .EQV. (P.0R. .N0T.P) 

2.9-2 Hierarchy of Operations in Boolean Expressions 
As for arithmetic expressions, the hierarchy of operations which may 
appear in a Boolean expression must he defined in order to establish a 
unique order in which the terms of the expression are evaluated and com- 
bined. Unless altered by parentheses, the order of operations p-rformed 
within a Boolean expression is given by the following list: 

Operation Symbols 

•ABS., + (equal hierarchy) 



*) /(eqiial hierarchy) 

+t - (squal hierarchy) 

.rj., .I'ici',, *(i., .u-jii., .jj.j .ijjCi. ^eq^uHx nierarcny^ 

.N0T. 

.AM). 

.flR., .EXjjSr. (equal hierarchy) 

.THEN. 

• EQV. 

Within an expression, operations of equal hierarchy are performed 
from left to right imless otherwise indicated by parentheses . 

Examples : 

(1) .ABS.(B - C) means |b - G|, while .ABS.B - C means |b| - C 

(2) - B + C means (-B) + (C), while -(B + C) means the negation 
of the sum. 
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(5) B.P. - X + Y means B""^ + Y, while B.P.(-X + Y) means B"'^"^^. 

(k) K2/Z - 5 means (K2/z) - 3, while K2/(Z - 3) implies that Z - 5 
is the denominator. 

(5) A * B + C means (A * B) + C 

(6) A. P. 3/ J means ik^)/j. 

(7) X.L. Y + 5 means (X) .L. (Y + 3). 

(8) P.AND..Nj6t.P .EQV.Q means (P.and. (.Nj6t.P)).EQV.Q 

2.9-3 Parentheses in Boolean Expressions 
Parentheses are used in the usual way in Boolean expressions . 

2.9-^ Mode and Boolean Expressions 
If an expression is a Boolean expression (see section 2-9.1, De- 
finition of Boolean Expressions), then the expression has a value which 
is of the Boolean mode: a somewhat looser terminology is to say that 
"the mode of the expression is Boolean" . 
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2-10 Statement Label Expressions 

A statement label expression may be any of the following: 

(a) h. statement label constant} 

(b) a statement label variable; 

(c) a function value of statement label mode. 
There are no other statement label expressions. 

If an expression is a statement label expression, then the expression 
has a value which is of statement label mode; a somewhat looser terminology 
is to say that "the mode of the expression is statement label". 

Examples ; 

(a ) If each of the following : 

LABEL 

C(3) 

XYZPDQ 

appears in the label field of some statement, then each one 
is a constant of statement label mode and hence each is a 
statement label expression. 

(b) If each of the variables 

X 

YZABG 

ijk(q) 

has been declared to be of statement label mode (see section 
3.2, Mode Declaration), then each is a statement label ex- 
pression. Notice that the entire array UK must be of statement 
label mode. 

(c) If the function value given by 

PAIR. (X,YZABC,IJK(Q) ) 

is of statement label mode, then this is a statement label 
expression. 
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2.11 Function Namg, Expressions 

A function name expression may be any one of the following; 

(a) a function -Harae, eonstant| 

(b) a function name variable^ 

(c) a function value of function name mode. 
There are no other fiinction name expressions. 

If an expression is a function name expression, then the expression 
has a value which is of function name mode| a somewhat looser terminology 
is to say that "the mode of the expression is function name". 

Examples s 

(a) If each of the following? 

SIN. 

ATAN. 

FUNCTN. 

is a function name, then each is a function name constant 
and hence each is a function name expression. 

(b) If each of the variables 

A125 

C2BX11 

LjfeA(l2,P) 

has been declared to be of function name mode (see section 3.2, 
Mode Declaration), then each is a function name expression. 

(c) If the function value given by; 

XYZF.(A1,A2) 

is of fimction name mode, then this is a function name expression, 
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2 .12 Summary of Expressions 

Pour kinds of expressions have been defined and may occur in the 
MAD language: 

Arithmetic Expressions, section 2.6. 

Boolean Expressions, section 2.9' 

Statement Label Expressions, section 2.10. 

Function Name Expressions, section 2.11. 

Of these, the Boolean expression is the most general since all of 
the others may appear as sub-expressions of a Boolean expression. 

When the terminology "any expression" is used, it is to mean any 
of the above four types of expressions . 
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3 DECLAHA.TIONS (Non-executable statements) 

Ttie ptirpose of declarations is to furnish information to the translator. 
With the exception of function declarations, all such statanents may occur 
anywhere in a programo Declarations may have statement labels, but^they are 

1 anoreri bv the t.ranslat.oi" anri mav not. he -ref AfpnceH hv n-i-.Vtoy e+Q + omon+cs 

3«1 Rsaark Declaration 

Each card containing a Remark declaratic© must have an "R" in colinm 11. 
The statement itself in coliaans 12-72 is any string of allowable characters 
(see seco 1.1, The Hollerith Character Setji and is ignored t^ the translator. 
It is reproduced where it occurred in the printed Msting of the source 
language program, thus furnishing information to the reader of the program. 
Example: 



Cols. 
1-10 



Colo 

11 



R 



Cols. 
12-72 



THE FOLLOWING IS THE 3RD OF A SET OF 5 » 3/12/62 



Colso 
73-80 
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3o2 Mode Declaration 

Variable and function values may be declared to be one of the following 
modes s 

FL0ATING P0IMT 
INTEGER 
B00mM 
FUNCTI0N NAME 
STATEMENT LABEL 

The form of the mode declaration is 

TT) u^, u^, » . o , u^ 

where /^ is one of the five modes listed above and each U^ , separated from the 
next by "j" is a variable name or function nameo In parti cular, no sub- 
scripted variable fonn may appear in the mode declaration. No U. may be 
declared to be of more than one mode throtaghout a program » 
Example St INTEGER ADDN» , Z5X, ALPHA 

B00LSAN BAR, ANAo, N32 

FUNCTT0KT NME BETA, CLT. 
MjiD assumes the mode of all variables and function names !:o be r'Li,/,,Tj.'lG 
POINT unless othervdse declared ■- This normal mode may be altered by a 
statement of the form 

KOMAL MODE IS /T) 

where 7>^ is one of the five modes listed above. Only one such statement may 
oocm- in a given program and this declared normal mode is iii, ef-'ect through- 
out the program, regardless of where it appears in the program-. 
Also: 

(1) All constants are assigned a mode consistent mth tr.axr 
respective forms- (See sen-. 2ol, Constants.) 

(2) A vector which is a dimension vector of some array in a 
DIMENSION' declaration is assigned INTEGER mode (see sec 3 = 6 = 2, MiiTRIX 
BIMMSI0N;-, 
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(3) A irector which is preset ixy a VEGT0R VALUES declaration is 
assigned a mode consistent with the assigned values (see -sec^ 3«7 VECT0R 
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3c 3 EQUIVALENCE Declaration 

The EQUIVALENCE declaration is of the forms 



EQUIVALENCE (V^, V^, . . . . V^). (V^^^, V^^^^ . . . , V ), 



^ n+p+q+1' n+p+a+2* '•'■'» ^+^^^4,^) 



where each ¥. is a variable name or a subscripted variable, the subscript of 
which is linear-. All the elements indicated within a pair of parentheses 
will be assigned to occupy the same storage location throughout the program. 
An array name V written without subscript is by convention taken to mean the 
element V(0, 0, » o , j 0)c. Ifcne array element from each of two different 
arrays appears within a pair of parentheses these are made equivalent and 
a one-to-one equivalence is thereby induced upon the overlapping rsnaining 
elements of the two arrays. One element from each of any number of different 
arrays may appear within a pair of parentheses. 

must be established by the appropriate mode declaration for each variable 
namee 

Any number of groups of equivalences may occur in any one EQUIVALENCE 
declaration (up to ten cards/statement), and any nimber of EQUIVALENCE 
declarations may appear at any place within a program. 
Example ss 

(a) A single statenren:. 
EQUIVALENCE (A, B, G2), (XLF, TSH), CP43A, XXX, ZZZ), (Q, R), (SIX, I6) 
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(b) A portion of a program: 

EQUIVALENCE (A, B, C2), (XLF, TSH) 



EQUIVALENCE (P43A, XXX), (Q^ R) 
EQUIVALENCE (XXX, ZZZ) 



EQUIVALENCE (SIX, 16) 



The resiilts of examples (a) and (b) are identical. 
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3 --4 PR0GHriM C0m0N Declaration 

The FR0GRM C0«40N declaration has the form 

where each V. is a variable name separated from the next by '%•♦, This 
declaration causes the specified simple variables and entire arrays to be 
stored in an area separate from the usual storage and separate from ERASABLE 
storage 5 section 3'5« These variables are not stored overlapping in storage 
as in EQUIVALENCE^ They are stored successively in order of their appearance^ 

Reoccurrences of the PR0GRAM C0M0M declaration do not erase the 
variables already assigned to PR0GRAK CfkM.0}i^ the new assigrments are appended 
to the previous list. 

(1) One use of this declaration is the provision for several 
sections of a program to refer to variables and arrays by the same naae , thus 
allowing the sections of the program to be translated and checked out 
independently of each other. A program segnented in this fashion wo*uld have 
the form of a main program plus several EXTERNAL FUNCTI0N programs, section 
6o5^ with the main program being used primarily to call on each of the 
EXTERNAL FUNCTI0N programs. Both the main program and the se^ents must 
contain the necessary PR0GRAM C0M0 declarations. Although variables and 
arrays to be used jointly by several EXTERNAL FUNCTIjZi^s can be communicated 
as arguments to the functions, assigning them to PR0GRAM C0OI0N makes them 
available with greater ease to all of the functions = The reservation of 
PR0GRAM Ci^»IM0N storage is performed only once and for the main program only; 
the PR0GRM 0^101^ declaration in the EXTERNAL FUNCTI0N segment allows the 
proper storage references to be made. If PR0GRAM C^I0N is to be used in 
this way (main program and EXTERNAL FUNCTI0Ns) the main program must be 
loaded into the computer xirat ^ 

(2) Another use for the PR0GRAM C0M0N declaration is in the 
situation where a program is so large that it cannot all be contained in 
the computer at onceo The program must then be written in segments, and if 
one segment is to use the results of a previous se^ent's computation^ the 
variables involved should be declared in PR0f;RAM 00^0^^ storage and will be 
retained throughout c The PR0GRM C0i4M0N and DIMENSI^ (sec, 3,6) declaration 
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which set up the storage allocation must be identical in all segnents which 
use these variables and arrayr. Additional variables and arrays may be added 
to the end of the PRjfcRAM C^0N list by any segiient. 
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3»5 ERASABLE Declaration 

The ERASABLE declaration has the form 

ERASABLE V, , V„, < . . , V 

. ^1 2' ; . n ■ 

where each V. is a simple variable or an array, separated from the next by 
","0 The variables and arrays assigned to ERASABLE are not overlapping as 
in EQUIVALENCE J but are put in a storage separate from the usual and separate 
from PR0GRAM C0M0No Each ERASABLE declaration deletes the effect of any 
previous ERASABLE declarations thus allowing variables and arrays to occupy 
the same storage at different times = 

Notice that external functions and translator defined functions may, 
and do, make use of erasable storageo Therefore care must be exercised if 
the programmer wishes data in erasable storage to rraiain intact after the 
operation of such a functiono The data which is to be preserved must appear 
as an entry in an ERASABLE declaration so that entries which appear to the 
left of it renresent at least an nnirh fj-rasabT <=» fif.nfa erf- a= t « +" Ke ncoH K-tr 
the function which is to execute in the interim and which requires the 
greatest amount of erasable storage » The entry or entries to the left of the 
critical entry in the iHASABLE declaration may be a dummy or dummy entries 
solely for the purpose of skipping over the non-safe erasable storage « 
Examples Suppose the following program segment appears? 



Z = 12o7 



R = SQRTc (A) 



^ = R + Z 
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where SQRT. is assumed to use two locations of erasable storage and X and Y 
are assumed to be expendable at least for this segnent. The variable Z is 
calculated before the SQHT. function is executed and is used afterwards. The 
expendable variables X and I, by virtue of appearing to the left of Z in the 

destroyed by the operation of SQET., thus preserving Z until required, since 
Z is assigned to a safe location one beyond those destroyed by SQRT». 
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3-6 DIMENSI0N Declarations 

Vectors and matrices may be declared in the same DIMENSI0N declaration^ 
The separation of their descriptions here is done for the sake of clarity-. 
If a variable is to be an array it must be declared in a D3MENSI0N decla- 
ration! this declaration need not appear in the program before the first 
use of the naaeo 

3=6=1 Vector DIMENSI0N Declaration 
The form of a DIMENSI0N declaration for vectors (one-dimensional arrays) 

is 

DIMENSI0N V^(p^), V2(P2), --, V^(p^) 

where each V. is a variable name and is followed by an integer constant 
enclosed in parentheses^ This integer constant is the largest value that 
the subscript of V. will assume during execution of the program= The size 
of the region reserved for the array V. will be this integer +1= The "»•• 
separates each declared vector from the nexto The subscript of an element 
of a vector array should not attain a value less than zero during execution 
of the programo 

3»6,2 Matrix DIMENSI0iJ Declaration 
Tie form of a DIMENSI0N declaration for matrices (n - dimensional 
arrays where n > 2) is 

DIMENSI0N ¥^(p^,D^(x^)), V^CPg.D^Cx^)), .«», VPm^V^^^ 

where each V. is a variable name (other than Statement Label) and is 
followed by two arguments enclosed in parentheseso (The general form of a 
subscripted variable which is an element of V. is 

i ^^1' ^2' ' ° - s ®j, ) 

,.rt^ere n.. is the number of dimensions of V.)o Each p., the first argument, 

as in vector DIMENSI0N above, is an integer constant whose value is the largest 

that the linear subscript of V. cam assumeo Each D.(x. )jthe second argument j 

defines ziie »'dimension vector" for the associated V. array in the foi Lowing 

way; 
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D. is a variable name| 

X. is an integer constant designating a specific element 



in the D. vector; 
D. (x. ) contains an integer constant, n. , \idiose value is the 

M J. J_ 

TTllilTin^''' '^T mtrion CT rvna r\-p if • 

— ■ ■—n i ..,,, . . ■,■■■ 1» 

D^(x^+1) contains the linear subscript (integer) of the V. 

array which is also to be the base element of the 
V. matrix, i.e., V(l, 1, > <, , , 1); 
D.(x.+2) contains an integer constant which is the largest 
value that the 2nd subscript of V. may assume; 



D.(x,+n.) contains an integer constant which is the largest 
value that the nth subscript of V. may assuneo 

(The lowest value of each subscript of V. is assumed to be 1). 

The dimension vector, Dj^(x^)> is autcanatically of INTEGER mode and may 
not be declared other than INTEGER. D. must itself be dimensioned, either 
by DIMENSI0N or by I:. T0R VALUES (sec. 3»7)5 D. must be dimensioned to have 
its largest linear subscHot > v.+n „ 

As is mentioned above, it is possible to declare an array in such a 
way that the linear subscript of the base element (the element with all 
subscripts = l) is gre-iter t^sr- '^ -^e-J^f -- ■'his is done it becomes meaningful, 
and it is permitted , to refer ;,: ...a elements with linear subscripts less 
than the linear subscript of tbe base element but not less than the linear 
subscript zero, using the multi-element" subscript (non-linear subscript) 
fonn with the values of the subscripts equal to zero and even less than 
zero. See example 3 below. 
Examples; 

(1) D2MENSI0N XA(400,V), V(3) 

XA is a 401-elanent array ndiose elements may be referred to as XA(0), 
XA(1), ..«, XA(400). V contains the 3-^enient dimension vector for the 
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matrix XA (V is itself a 4-element array and is automatically of INTEGER 
mode) I if V has been preset as follows ; 

V(0) = 2| 
V(l) - 6| 

v(2) = r^ I 

V{3) may De used as an integer variable elsewhere | 

and if XA(Ij J) is any element in the matrix Ikt 

(a) XA has 2 dimensions^ since V(0) ~ 2| 

(b) the base element of XA, XA(l^l)j is XA.(6), since V(l) = 6| 

(c) the allowable range of J is -5 < J < 13, since ¥(2) = 13o 

So that the elements XA(0) through XA(5) are not "within" the matrix XA, 
and XA(6}, XA(7), =-, XA(395), »-, XA(400) may be referred to as XA(l,l), 
XA(1,2), 0,, XA(30, 13), coo, XA(31, 5)o Although XA(0) through XA(5) are 
not "within" the matrix XAj the following is trues 

XA(5) = XA(0jl3) = XA(1,0) 
XA(4) ^ XA(0,12) ~ XA{1,»1) 
XA(3) f XA(0,11) S XA(l,-2) 
XA(2) E XA(OAO) = XA(l,-3) 
XA(1) i XA(0,9) = XA(1,.4) 
XA(0) £ XA(0,8) = XA(l,-5) 

(2) DIMENSI0J Y5(250s J(5)), J(15) 

Y5 is a 251-element array whose elements may be referred to as Y5(0), 

Y5(l), =c.», Y5(250) = 

J contains the 5-element dimension vector for the matrix Y5| if J has 

been preset as follows? 

J(0) through J(5) may be used as integer variables elsewhere| 

J(5) = 41 

J(6) = 20| 

J(7) = 4? 

J(8) = 6| 

J(9) = 3| 

J(10) through J(15) may be used as integer variables elsewhere | 
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and if Yi>(l, Kj Lj M) is any eletaent in ifie Y? matrix? 

(a) Y5 has 4 dimensions, since J(5) = .4| 

(b) the base element of 15, Y5(l,ls^.l), is 15(20)^ since J(6) = 20j 

(c) the allowable range of K is 1 < K < 4? since J(7) '- 4s 

(d) the allowable range of L is 1 < L < 6, since J(8) = 6: 

So that the elements I5(0) through Y5(19) are not "vdthin" tiia matrix Y5> 
and 

Y5(20) £ Y5(lA,lA)l 

Y5(22) E Y5(1,1A,3)| 

Y5(23) £ Y5{1,1>2,1)| 

Y5(35) 5 Y5(1A,6,1)S 

Y5(37) r Y5(l,l,6,3)| 

Y5(38) H Y5(l,2,l,l); 

Y5(91) r Y5(l,4,6,3)? 

Y5(92) = Y5(2AA,1)5 

Y5(235) ^ Y5(3, 4,6,3)5 

and Y5(250) 5 Y5(4A,533)^ 

NOTE li The general formiola for computing a linear subscript L of the 

array X from a known matrix subscript (w, , w , «o», w. , <.»», w ), 



where the upper values of w^ through w^ are W^, W^, ° -> » j W^ ) = " ° > 

W , and wi 
n' 

is B, is; 



W , and where the linear subscript of the base element X(l,l,<.<. » ,1) 



n-1 i-1 

L = B + (w -1) + I C(w . - 1) TT W .] 

" i=i ""-^ j=o "-J 

E«g. , given Y5(3s4c6,3)^ Y5(w, , w^, w„5 w, ) above, 

X /d ^ 4 

L = 20 + (3-1) + [(6-l)(3)3 + C(4-l)(3)(6)3 
+ [(3-l)(3)(6)(4)] = 235 
and Y5(3,456,3)5Y5(235)» 

NOTE 2: An algorithm for computing for the matrix X a subscript 

(w, , w , »o = , w ) which corresponds to the linear subscript L, 
where the linear subscript of the base element is B, and the upper 



bounds on w„, w_, <.<.o, w are W_, Wo» •»■>, W iss 
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(1) L - B = R 



(2) -I = Q^H.-!i^,Q^^l=,^ 



TT w . TT w . 

i=2 1=2 



(3) i_,Q +-!2__ J Q2'°^ = ^2 
^^' n ^2 n 

TTW.; TTW. 

i=3 1=3 



r r 

(4) -^ — = Qt + ■— ^ — ; G.. i- 1 = w„ 

^ ' n 3 n ^ 3 

77" w^ TT w^ 

i=4 1=4 



r ^ r , 

^ ■'■ n-1 ;*' ' ri-1 ri-1. 



n n 



(n+1) r -, + 1 = w 

E.g., given Y5(250) and W^ = 4, W^ = o, w^^S W^ = 3 abov^: 

(1) ■1^0..2'.: = 230 

(^) 230/72 = 3 + ^ ; 4 - v;^ 

(3) 14/18 = + i| 5 1 -^ w^ 

(n) 14/3 = 4 + I ; 5 = w„ 

(n-hl) 3 = v'„ 



3.6.3 Automatic DIMENS10N 
Thers are two cases in which MAD performs automaiieally the nec<3scarv 
D]iJiEI\iSI0Jlng : 

(1) If L is a trbatoment label name, and a is the Lar/^cst si)>;,;cript 
v/hich appears v/ltn L ;i.a tiiu statbment label I'juid vS6c= 2.3? 
Statement Labels) then an automatic DB4ENSI0N L(n} occurs, i.e., 
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n+1 locations are reserved for the L vector. (No harm is done if 
L is also dimensioned by the programmer). 

Examples If the following appear as labels of various statements 
throughout the program; 

LABELV(l) 
LABELVO) 
LABELV(7) 

and LABELV does not appear in a D1MENSI0N declaration, then 8 
locations will be reserved for the IjABELV vector aut<iimatlcally. 
(2) If part or all of a linear array is preset by VECT0R VALUES 

(sec. 3»7), the array need not be declared by DIMENSI0N unless 
the size of the array automatically reserved by VEbT0R VALUES 
is not sufficiently large. 
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3' 7 VECT0R VALUES Declaration 

Any vector or portion of a vector (or array, when using its linear 
subscripts) may be preset by a VECT0R VALUES declaration of the foraiJ 

VECT0R VALUES (/'= c^, c, , ,.., c 

1^ ' n 

where ^is any variable V or linearly subscripted variable V(m)| each 
c^ is any constant, and all c^ are of the same mode. The elements of the 
vector V through V(n), or V(m) through V(m+n), are preset with the values 
Cq through c^ at compilation time~ The mode of the vector (/^is auto- 
matically assigned by MAD to be that of c. , and may not be declared to be 
a mode other than that of c^. A region of n+1 (or m+n+1) storage locations 
is reserved for ^^ but this region may be enlarged by a DIMENSION declar- 
ation or by another VECTOR VALUES declaration. 

For alphabetic constants, the effects of VECT0R VALUES may be extended 
to the following: 

VECT0R VALUES i/= $c^c, c_. . .cj, $c ^, c ^_c ^^..»c $, 

'' L 2 n n+1 n+2 n+3 n+p ^ 

°°'> ^Vp+q+lVp+q+2"^-Vp+q+r* 
where there may appear between $'s a string of characters (see sec. 2«1.3 
for allowable characters in alphabetic constants) of any length. Each 
c. is considered to be a group of six characters o If the last c. between 
$'s does not contain six characters, blanks are appended to the right to 
make a group of six<= 

It is allowed to intermingle integer constants and alphabetic constants 
in any VECT0R VALUES declaration, e»g,-, 

VECTCfe VALUES BETA = 42, $THISbISbNUMBERb42,$, 3 

and BETA will have reserved for it five locations which will be preset as 
follows s 

BETA(O) = 42 

BETA(l) = THISbl 

BETA(2) = SbNUMB 

BETA (3) = ERb42 

BETA(4) = 3 
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Example: 

VEGT0!R values ALPHA(32) - $33HOTHISbSITUATt^ 

CRSATESbANbERR0R»$ , 
$l6HbRBADlM0REbDATA*$ 
The vector ALPHA, elements 32 through /*!, are being preset to be used 
aS & tuSSSags at ex6cutl.en ujjues The chaPactsi's uebwesn ^'a u@fxn& 
print f omat s ( sec . 5 » 5 • 3 , PRINT F0M AT Stat^ent ) . 

Forty-two locations are reserved for ALPHA, ALPHA(O) through ALPHA(4l)» 
Values will be preset as follows? 

ALPHA(32) » 33HOTH ALPHA(33) " iSbSlT 
ALPHA(34) - UATI0N ALPHA(35) •= bCEEAT 
ALPHA(36) = ESbANb ALPHA(37) = ERR0R* 
ALPHA(38) - 16HbRE a1PHA(39) - ADbM0R 
ALPHA(40) = EbDATA ALPHA(41) « «bbbbb 
In order to print these messages, the program would have to coritain 

PRINT F0HMAT ALPHA{32) 
PKTNT F0KIAT ALPHA(38) 

RESTRICfTIONS: 

(1) Vectors which have been declared in ERASASyS storage may not be 
preset by VECT0R VALUES. 

(2) Vectors which have been declared in PR0GRAM C^4M0H storajge may not 
be preset with statement labels or function names; also, these v«ctors 

may be preset only in a 1-seetion program, or, in an n-section program if 
the PR0GRAM C04M0J region is identical in al3 n sections. 
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3^8 Single Statement INTERNAL FUNCTI0N Definition 

A description of the procediores necessary to enable the user to define 
more general t3rpes of functions appears in sec= 69 Function Definitions^ 

The single statement internal function definition is the simplest 
type of ftmction definition = Since this is an internal function j it is 
translated as part of the main program =.• This single statement has the 

fOIHlS 

INTERNAL FUNCTION f. (A^jA^s- » = j\) = E 

where F- is a function name (see sec. 2^k^ Functions), the A. 's are "dummy 
arguments" and E is any expression which is consistent with the declared 
mode of F- = 

The "dummy arguments" are used in the expression E to indicate the 
correct correspondence between a variable name used in E and a position in 
the list of arguments.- When the internal function F. is used , every 
occurrence of the "dummy variable** A. in the expression E will be replaced 
by the value of the arg'ument which occupies the position in the list of 
argisnents of the "disnmy argument" A. -- 

The form of a diasmy argianent may bei 

(1) a name of a simple variable, 

(2) a name of an array variables, 

(3) a name of a function-. 

Notice that neither constants nor subscripted variables may appear as 
dxggmg argiments and that arguments of a function do not appear following 
the function name in the dummy argument list.- 

Names used as dijmmy arguments must be distinct from all other names in 
the program^ Names vdiich appear as dtamay argiments may not appear in any 
PR0GRAM C0®10N, ERASABLE or EQUIVALENCE declaration. 

The modes of the dummy arguments must be declared as for other variables 
if of other than normal mode: 

Dxmmy arguments which are array names need not be dimensioned = 

The alphan'oaeric portion of the name F= of the defined function must be 

distinct .from all ether name? v,s9d in the program, and from the names of the 

functions already available t.o the translator. For a list of these names see 

the manual of library subroutines: 
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Ordinarilys all the names of the dummy argnnent list will appear in the 
expression E (otherwise, the unused ones need not appear as dummy argvanents)., 
Names of variables or functions which do not appear as dummy arguments, but 
which are defined elsewhere in the program, may also appear in the expression 
E. The value of such a variable or function in the expression is the current 
»cu.uc exo one ujliuc one ^kuiuoj.uu -T • 0.0 uopu , oxwicr oj one p.APf .11 1 P. stattaaent 
or as a term in an expression. 

In the use of a function an argiment may be any expression which. agrees 
in mode with the corresponding dummy argument <. 

A single statement internal function definition may aj^ear any\ifhere 
within a program except in another internal function (see sec« 6^4, Internal 
Function Definitions). 

In the example: 

INTERNAL FUNCTION P0LY. (N^ X, FN.) = FN- (J»X).PcN - X/XBAR 
Tirfiich might be used in the statement (the statement label is BETA): 
BETA ZQ = P0LY. (M + 1, Y, SIN» ) + P0LYo(M - 1, Z, C0SO 
it is understood that if N is in the integer mode, then so is.M, and if X 
is in the floating point mode, then so are Y and Z. Both M and N would 
have had to be declared to be in the integer mode, of course. Similarly, 
the values of SINo and C0S<. must be the same mode as the values of FN. 
Moreover, in the use of functions, this mode correspondence cannot be 
checked by the translator • 

The function P0LY. has as one of its arguments the name of a function. 
In the statement BETA the function used in the first term to the right of 
the "=" sign is SIN. and in the second term C0S.. Hence, statement BETA 
is then equivalent to: 

ZQ = SIN.(J*Y).P.(M+1) - Y/XBAR+C0S.(J*Z).P,(M-1)-Z/XBAR* 
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4 Basic Executable Statements 

4ol Substitution Statement 
This statement has the form: 



variable or a subscripted array variable, and the right side, F, consists 
of any expression of the same modeo The only exceptions to this mode 
requirement are the cases: 

(1) If the variable on the left is of integer mode then the value 
of a floating point expression on the right will be converted to integer 
mode • 

(2) If the variable on the left is of floating point mode then the 
value of an integer expression on the right will be converted to floating 
point modeo 

These conversions are the only automatic ones and any other disagree- 
ment of mode is not permitted. 

The substitution statement is to be interpreted in the f ollovdng way: 

(1) compute the value of the expression on the right, 

(2) convert it, if necessary and possible, to the mode of the variable 
on the left of the "=" sign, and 

(3) give the variable on the left the value which results frcm steps 
(1) and (2). 

Thus, if Y is a floating point variable, then the statement 

Y = 1 

will cause the integer 1 to be converted to floating point and then stored 
in the location called "Y"! ioe., Y will now have the value lo (as a 
floating point number). If the statonent were written 

Y = 1. 

then the floating point number lo woiild be stored in the location "Y"; 

i.e., Y would again have the floating point value 1., but in this case 

the conversion of the integer is unnecessary, thus speeding up the 

ccanputation. 
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When a floating point number is to be converted to an integer, it is 
first expressed as a number with both an integer part and a fractional part, 
and then the fractional part is trimcated. Thus, the following floating 
point ntmbers: 

3E5, .3E0, .34568127E2, - .345681E10 
would convert to the following integers, respectively: 

300000, 0, 34, - 3456810000 
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4.2 TRANSFER T0 Statement 

This statonent has the form: 

TRANSFER T0 <^ 

Here^jTmay be any expression in statement label mode and in particular 
may be any statement label. Execution of this statonent causes the 
computation to continue from the statement vdiose label is the value of (3; . 
Examples: (1) TRANSFER T0 SUMX 

(2) TRANSFER T0 SWTCH {K+2) 

If E=»4 then the valw of SWTCH (K+3) is SWTar<6i>. 
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403 Conditional Statements 

There are two types of conditional statements-. 

4*3<'l Simple Conditional 
This statement has the forms 

WHENEVER B^ Q 

Here B is a Boolean expression and Q any executable statement except 
the follovrings END 0F PR0GRAM, another conditional, THR0UGH or EXECUTE. 
If at the time of execution of this statement, the expression B has the 
value IB, the statement Q is executed^ If, however, B=OB, then Q is skipped 
and ccanputation continues from the next statement in sequence. The comma 
in this statement, as in other statements containing punctuation, must be 
written. 
Examples? 

WHENEVER 3MoLEol, TRANSFER T0 END 
WHENEVER I.Cffi.NcAND.J^NE.I-l, I-O 

4»3=2 Compound Conditional 
This type of conditional has the forms 

-^ WHENEVER B^ 



J^ 0R WHENEVER B, 



^2 



(J^ 



WHENEVER B, 
k 

"1- "k 



k+1 END 0F C0NDITI0NAL 
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Often the last condition B, expressed is one for which the condition is 
always true.. This may be expressed by the statement 

0R WHENEVER IB 

or alternately the statement 

0THERWISE 

The . are statement labels ^diich need not be used unless desiredj k 
may eqiial 1 (if no "0R WHENEVER. •<■ " fetatanents are used). Here B^, c, 
B, are Boolean expressions, and the execution of this block of statements 
is as follows: 

Each B. is tested in turn, starting with B-. . If B.. has the value OB, 
then Bp is tested, etc» As soon as some expression, say B., has the value 
IB, then the statements between (but not including) «/. and^. ,(i.e. 0.) 
are executed. At this point the execution of the entire block is con- 
sidered ended, and computation continues from the first statement after 
the declaration labelled ^ ■u±.y I^i other words, at most one of the 
alternative computations is performed; e.g., that one which immediately 
follows the first expression B. -which has the value IB. If none of the 
B. has the value IB, none of the ccanputation in the scope of these state- 
ments is perfonaed. 

Among the, statements of 0. for any i there may appear other conpound 
stataiients. The maximum pennissible nesting depth of compound conditional 
statements is 30« 

rlxarriple: The evaluation of the function vrfiose graph is 
141 








\ 

1 
I 

I 
I 



rrdght be given by the section of the program: 

WHENEVER X »IiE. 0. .0R. 1. .LE. X .AND. X .L. 2. .0R. X oGEO 

Y = 0. 

0R WHENEVER 0. .LE. X .AND. X.L.I. 

Y = X 
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0R WHENEVER 2. .LE. X =AND. X .L« 3. 

Y = 1. 

END 0F C0NDITI0NAL 
This section of program could be rewritten in another way. 

WHENEVER O..LE. X .AND. X .L. 1. 

Y = X 

PR WHENEVER 2..LE. X .AND, X .L. 3- 

Y = 1. 
jDTHSRWISS 

Y = 0. 

END 0F C0NDITIS2^AL 

The indentation of the statements between the conditional statoaents is 
not reqiiired but contributes to the readability. 
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4.4 C0NTINUS Statement 

This statement has the simple form: 

CONTINUE 
It serves as an entry point in the program, and causes no computation 

THR0UGH Statonent, illustrates the use of CONTINUE. 
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4»5 THH0'UGH (Iteration) Statement 

This statement causes the block of statements which follows immediately 
afterwards to be repeatedly executed, each time varying the value of some 
variable, until the specified list of values for that variable is exhausted, 
or until some specified condition is satisfied. The THI^UGH statement may 
take either of the following two forms^ 



4«5»1 THR0UGH J, F0R VALUES 0F V = E, , E^, .«», E 

D 12m 

Here ^ is the statement label of the last executable statement in the 

block to be repeated^ The block of statements following (and not including) 

the THROUGH statement, up to and including the statement whose label is,g/s 

will be called the "scope" of the THR0UGH Statement, Following the word 0F 

appears the name of the iteration variable Vj which may be either an 

individual variable or subscripted array variable of any mode. To the right 

of the "=•• sign may appear any number of expressions E, , o... E . The modes 

1 m 
of the E. bear the same relationship to the mode 1 as they would in the 

Substitution statement V = E. (see sec 4.1), 

The execution of this statement causes the statements within its 

"scope" to be executed, first with V = E, , then again with V = E , and so 

on, imtil the list of expressions^ is exhausted. Ccanputation then proceeds 

with the statement immediately following statement ca/ • At this time the 

iteration variable will have the value of the expression E unless its 

m 

value was changed during the final iteration. Should a transfer be made to 
another part of the program at any time during the iteration, V will have 
its current value = 

An example of this type of statement i-ss 

THR0UGH A, F0a VALUES 0F BETA = 3, 4, X5, Y(6 + I) +.3 
J = 5 * BETA + 6 
Jl = J .P.. 5 - 1 
A X(BETA) = Jl * C0Sa(2. * THETA) 

4"5»2 THR0UGH J , F0R V = E, , E^, B 

Here g^is a statement label which defines the "scope" exactly as in 

(sec. 4=5'1) above (with the exceptions ifg/is the label of the THR0UGH 
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statanent itself, the iteration will proceed as described below, but the 
scope will be empty, and the iteration will consist only of the incrementing 
of V by E^ ^d the testing)^ Following the word F0R is the name of t?ie 
iteration variable V which may be an individual variable or subscripted 
array variable of any mode. The modes of E, and E^ are subject to the rules 
which would apply to the substitution statements V = E^ and V = ? + E^. B 
is a Boolean Expression. 

The execution of the statement proceeds as follows: The variable V is 
set equal to E, • If the value of B = IB, the scope of the THR0UGH statanent 
is not executed. If the value of B = OB, the scope is executed. V is then 
incremented by £„, and B is tested again. In general, as soon as B = IB, 
the scope is not executed, and the computation proceeds from the statement 
immediately following statement (^ . Each time B = OB, the statements in the 
scope are executed, then V is incremented by E , and B is tested again. 
Thus, when the iteration is finished and B = IB, V has the value used during 
the last computation of the scope, incremented by E_. The scope will not 
have been executed for this value of V. (The value of V will be E, , of 
course, if B = IB before the scope is executed at all.) If, at any time, 
the ccsnputation transfers out of the iteration to another part of the 
program, the value of V will be the current value at the time the transfer 
was made. 

ixampj-ol ihe following program segrient illustrates this type of THR0UGH 
statement . 
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L = 1 
K = 1 
A = D(l,l) 

THR0UGH STl, F0R I = 1,1, I ,G. 10 
THR0UQH STl, F0R J = 1,1, J .G. 10 
WHENEVER A .LE. D(l, j), TRANSFES T0 STl 
K = I 
L = J 

A = D(I, J) 
STl GC^NTINUE 



THROUGH STATa^NT 

This program will locate the algebraically anallest element of the 
ten-by-ten, 100 elanent, matrix Q and the column and row subscripts of the 
first such element , if the smallest element occtirs more than once. 
The statements 

L « 1 and 
K = 1 

initialize the column and row subscript indicators to 1. 

The statement 

A = D(l,l) 
initializes the variable A, which is to hold the "smallest elanent so far", 
to the base element of the matrix* 
The statements 

THR0UGH STl, F^ I = 1,1, I .G. 10 
THS0UGH STl, F0R J = 1,1, J .G. 10 

define the iteration. The first will count from the first row through 
the tenth row and the second will count from the first colunm through the 
tenth colunn. Since the second THROUGH statement is nested inside the first 
it will count through the 10 columns for each of the 10 rows. 
The statement 

WHENEVER A .LE. D(l, J), TRANSFER 10 STl 

makes the canparison between A, the smallest so far, and the next untested 
element D(I, J). If A is less than or equal to D(I, J) it is still the 
smallest so far, so control is transferred to statement STl. 

If A is greater than D(I, j), transfer is not madej instead the state- 
ments 

K = I 

L = J 
are executed which update the row and column subscript indicators to the 
subscripts of the new "smallest element so far" and the statement 

A = D(I, J) 

updates the "anallest element so far" variable A, to the new "smallest 
elanent so far". 
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The statament 

STl C0NTINUE 

resxilts in no operation. It does, however, provide a statement inside the 
scope of the iteration to which the WHENEVER statement can transfer to 
continue the iteration but which will change nothing else. 

In all cases ever""" reference to an ex'^ression E. will involve its 
current value at the time of reference. Moreover, the variable V may have 
its value changed at Kiy time during the execution of the scope. In a 
statement of the form of section A.5'l> the value of V will be reset by 
the value of the next E. for the next computation of the scope. In a 
statement of the form of section 4* 5*2, the current value of V will be used 
for incrementing, testing, etc^ 
Examples: 

n n 1 

(i) To evaluate the polyncmial ex + c _,x ~ +...+c,x + c using the 

formula (.,.((cx+c ,)x+c ^)x +...+ c, )x + c (nested multiplication), 
n n-l n-2 1 o 

we may write the program: 

INTEGER J,N 

Y = 0. 

THR0UGH P0LY, F0R J = N, -1, J .L. 

P0LY Y = X ^- Y + C(J) 

^x' t^x uiic uLc;cuia.xiK '^■J* wiles oua-u^uciib ^m j.xj\jnjj.b v ^it^ *j^^ o^v.^« ^<^^« 

(ii) A Newton's Method solution i\^j = x, - jj-f — T ) °^ ^^^ equation 
f(x) = cos X - X = could be written as a single statement, using the 



criteria "jf(x)| < e and |x - x, , | 



(X,) 



fTZT 



iteration: ■ 



< e " for stopping the 



NEW THR0UGH NEW, F0R X = XO, (C0S.(X)-X)/(SIN. (X) +1.), 

1 .ABS.(G0S.(X) - X) .L. .EPSL0N .AND- 

2 .ABS.((C0SX=(X) - X)/(SIN.(X) + 1.)) -L. EPSL0N 

where XO is the initial guess. The digits "1" and "2" which. appear at the 
left of the second and third lines serve to show that the second and third 
lines appear on second and third cards; i.e., the statement is "continued" 
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from the first card and occupies three cards. See sec. 14. 1.2. 2, Continue 
ation Cards. 

(iii) If the transfomation of the iteration variable is to be per- 
formed within the scope of the iteration, one may use a zero increment. 
E.g., if J is an integer variable, and the scope of the iteration is to be 
performed for those multiples of 2 which are not multiples of 5 and which 
are less than the value of the integer K, one might write the iteration as 
follows: 

Hm^UQH END, F^R J = 2, 0, J .GE, K 



J = J + 2 

END WHENEVER J .E. (j/5) * 5, J = J + 2 

(iv) A table-look-up procedure using an iteration statement. Suppose 
that a string of alphabetic (or numeric) characters (i=e=j a "sentence"^ 
has been decomposed into single characters stored in 0(1), C(2), ..., C(K), 
where K is the length of the string. Then the first occurrence of a ccanma 
coxild be found as follows: 

Lm .THR0UGH L00K, F0R I = 1, 1, C(I) .E. $,$ .0R. I .G. K 

WHENSVSR I =G< K, TRhNSFER T(/( }^0C0iiA 
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4o6 Nested THR0UGH Statements 

As indicated in sec. 4»5» the "scope" of an iteration statement is the 
block of statements designated for repeated executions 



r 



THR0UGH END, F0R V = E^, E , B 



jcope 



END 



Seme of the statements within the scope of an iteration may, them- 
selves, be iteration statements. However, if iteration statement b is in 
the scope of iteration statement a, then the scope of b must be entirely 
within the scope of a. The following diagrams represent scane valid con- 
figurations: 

(1) THR0UGH K. F0R . . . (Iteration statement a) 



Scope 

of a 



Scope of 



THH0UGH M, F0R = . . (Iteration statement b) 



M 



(2) 



Scope of 
a 



Scope of 
b 



THR0UGH K, F0R ... (iteration a) 



THR0UGH K, F0R . . . (iteration b) 



K 



Here, although the scopes of a and b both end on statement K, iter- 
ation b is incremented and tested first. Therefore, iteration b is com- 
pleted before iteration a is incremented. 
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(3) 



Scope 
of a 



Scope of 



THR0UGH K, F0R « . . (Iteration a) 

« ft 

THR0UGH M. F0R »o, (Iteration b) 



M 



Scope of 



THROUGH M. F0R ... (iteration c) 



N 



K 



The 



(4) 

Scope 
of b 



following diagran represents an invalid configuration J 

Scope of THR0UGH K, F0R . . . (Iteration a) 



THR0UGH M^ F0R 1 1 1 (Iteration b) 



M 



When iteration statements occur in the scope of other iteration 
statements, they are said to be "nested." The "nesting depth" of an iter- 
ation statement is the number of iteration statements in whose scope it 
appears. The nesting depth of an iteration may not exceed 50. For example: 



(5) 



Scope of 
a 



THR0UGH K^ F0R ... (Iteration a) 



Scope of 
b 



THR0UGH M, F0R ... (iteration b) 



Scope 
of c 



THR0UGH M, FOR ... (iteration c) 



N 



M 



K 
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In example (5)s iteration a has a nesting depth 0^ iteration b has 
nesting depth 1, and iteration £ has nesting depth 2. In example (3)^ both 
b and _c have nesting depths of 1^ 

There are no restrictions on jumping into or out of the statements in 
the seope of an iteration. If the program Jumps out of the iteration and 
the iteration variables are not modified while outside the iteration and 
if control is rettirned to the stat«nent in. the iteration following the 
"jump out" statement, then the execution will continue as if it had not 
been interrupted « 
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4.7 PAUSE N0. Statement 

This statement has the form: 

PAUSE N0. n 

This statement indicates a breakpoint in the program, and it causes 
the computer to stop in such a way that the operator can manually start it 
and automatically go on to the next statement in the program. The number 
"n" is an octal integer containing up to 5 digits, vrtiich will be displayed 
on the computer console for the operator to note when the computer stops, 
thus indicating the point in the program at which the stop occurred. This 
statement is used only in very special circumstance sj the majority of users 
will have no need for it. 
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4'. 8 EXECUTE Statement 

The EXECUTE statement may have either of two forms i 

(a) EXECUTE C« (A,, A , .,,, A ) 

T«*iere C is the name of the fimctien and the A. are peimissible arguments 
for the particular function (see sec. 2»4, Functions) | and 

(b) EXECUTE C. 

where C. is the name of the fxmction and the particular function does not 
require airgvments. 

This statement is provided to pennit non-single valued functions to be 
executed. 

Examples To execute the function of section 2.4.2, Non Single Valued 
Functions, which sorts a list, write: 

EXECUTE S0RT. (LIST, N) 

It ie not meaningful to write 

EXECUTE SIN. (X) 
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4.9 END 0F PR0GRilM Statement 
This statement has the fonns 

END 0P PR0GRM 

This statement must be the physically last stateiaent in the program 
(i.e., the last carci)| it serves t& teminate coSipilation of the prograa. 
If it is executed it will terminiate the computation of the object program. 
Execution of this statement at execution time will transfer control to the 
operating system in which the ti^ansiatecl program is iMbedded. An Alternate 
way of terminating the execution of a program - 1*6* returning to the 
operating system - is to attempt to exeowte an input statement *^en the 
data has been exhausted (see s§c= 5»5, Basic Inpiit - Output StAteMents)^ 
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5 Input and Output 

In this section the various statements^ specif ieations, and lists which 
may be used to read in fresh data, output results, and save and read data 
onto and from auxiliary storage _, are described. 

5.1 Simple Input-Output Statements 

These statements are referred to as Simple Input-Output Statements be- 
cause their use does not involve the use of Format Specifications a.nd Input- 
Output Lists (sections 5-2 through 5'3)' 

Of these^ the input statements are used in order to read any data which 
may appear on cards following the user's program in which these statements 
appear. The output statements are used in order to cause results calculated 
when the program is executed to be printed on paper (listed) or punched on 
cards . 

5.1.1 READ DATA Statement 
The form of this statement is 

READ DATA 

This statement causes information to be read from cards . The values 
to be read and the variable names are punched in the data cards in a se- 
quence of fields of the form; 

V^ = n-^, V2 = ng, V^ = n^, ■•■> \ = \* 

The V,!- ..,, Y, are the variable names and n, j, ..., n, are the corresponding 
values . Reading is continued from card to card until the terminating mark 
* is encountered. Fields cannot be divided between cards, so that the last 
character in a card not terminated by an asterisk would normally be a 
comma . However, as a convenience, the end of the card is treated as an 
implied comma and hence this final comma may be omitted . The variable 
names may designate simple variables or elements of linear and two- 
dimensional arrays. The subscripts on the array variables must be integer 
constants . The values may be floating point, integer, octal. Boolean, or 
alphabetic with the forms described in section 2.1, Constants. 

5-0 
6/20/62 



SIMPLE INPUT-OUTPUT STATEMENTS 

For convenience in entering values of array elements, it is 
possible to designate only one variable name and have successive 
nximbers, written without names , interpreted as the consecutive 
values of the array, i.e., 

V(j) = n-^, ng, n ; .. ., n^^ 

would be the same as 

V(j) = n^, V(j + 1) = Hg, ..., V(j + k - 1) = n^ 

For two-dimensional arrays, successive numbers will be entered 
in succeeding columns of the designated row until the row, as de- 
termined from the current value of the dimension vector, is filled, 
and then the next row will be started . 

Zeros must be punched; adjacent commas (,,) are simply skipped. 
Blanks are ignored throughout except between dollar signs (which are 
used only to delimit a string of Hollerith characters). 

A simple integer variable may take on a value equal to a group 
of six or fewer Hollerith characters delimited by dollar signs, 
i.e., an alphabetic constant. 

Longer strings of Hollerith characters may be entered as elements 
of arrays. Such strings are divided into six character groups 
for storage . 

As an example illustrating many of the features described herein 
consider the data card sets 

XI = 1.2, Yl = -6.8, INDEX = k, A(U) = J.l, -10.95, 
12.6, MATRIX (2,1) = 25E-2, 1.8E-10, 3-l^E-8, 
STRING (1) = i:iEKDb0FbPRj6BLEM$ * 

A MAD program deck having in its body a READ DATA statement 
would be followed by the preceding data cards . At the time the 
READ DATA is executed (not compiled), these three cards would be 
read and the values in storage of the computer would be set as 
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XI wovild become 1.2 

Yl would become -6.8 

INDEX would become k 

A(k) -would become 3»1 

A (5) would become -10.93 

A(6) woxold become 12.6 

MATRIX (2,1) would become 25E-2 

MATRIX (2 ,2) would become 1.8E-10 

MATRIX (2, 5) would become 3.1^E-8 

STRIHG(l) would become EMJbjZjF 

STRING(2) would become ERjZ^BLE 

STRING (3) would become Mbbbbb 

Note that although the components of the data cards resemble sub- 
stitution statements (Xl=1.2, Yl=-6.8, etc.), these cards are not parts 
of the source language jsrogram and no attenrpt should be made to compile 
them. 

5.1.2 mWI EXULTS statement 

The form of this statement is 

PRINT RESULTS ^ 
where (_, is a list of the form 

1^ 2-' *''^ n 

and the T.'s are terms of the list. The T.'s may be 

(1) simple variables J 

(2) subscripted variables j 

(3) block designations of the form 

V(s^^ s^, ..., s^) ... Y(r^, r^, ..., r^^) 

where the .V(i, , i , ..., i ) are subscripted variables. 
Note that expressions n®y not appear as terms T pf the list but subscripts 
of subscripted variables in the list may be expressions (see section 2.2.^^ 
Subscripts of Subscripted Variables). 

5-2 
6/20/62 



SIMPLE INPUT-OUTPUT STATEMENTS 

The printed output is analogous to the data input in that the numbers 
printed are preceded by the appropriate variable name and an equal sign, 
e.g., X = -12 A, and only the initial elements of arrays are so labeled. 1 
Elements of three and higher dimensional arrays will be labeled with the 
equivalent linear subscript- If dummy variables (in a function definition) 
are included in the list the specific variables assigned to such variables 
during execution will not be labeled but simply preceded by ... 

Example statements are: 

PRINT RESULTS XI, Yl, Z(l), ..., Z(N+1), MTX(1,1) . . .MTX(HrN) 
PRINT RESULTS X (1,3,^,2), ..., X(N+l,i^,4,l) 

5.1.3 PRINT CjZ^MMENT Statement 

This statement has the form; 

PRINT C)Z)MMENT $S$ 

Here S designates a string of no more than 120 Hollerith cl^racters . 
These characters may not Include dollar signs and here blanks are not 
ignored. The string, delimited by dollar signs as indicated, will be 
printed. The first character will be interpreted as a carriage control 
code (see section 5.2 .1^.2 (b)) • 

An example statement is ; 

PRINT COMMENT $1 J0HN PUBLIC, MATH 373 PROBLEM 1$ 
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5.2 Format Specifications 

When information is read from (or punched into) a card into (or 
from) a computer^ it is necessary to know how this information has been 
allocated among the available columns of the card. Similarly, whenever 
information is to be printed by a printer (either on-line or off-line), 
it is necessary to know how this information has been allocated among the 
columns available on the printer. A description of each allocation is 
called a format specification. Usually, but not always, a list of variables 
(see section 5.3, The Input-Output List) whose values are to be printed, 
punched or read, is associated with a format specification. (Occasionally, 
the information is contained entirely in the format specification, so the 
list may be empj;y.) 

In substance, a format specification is a string of alphanumeric 
characters and a restricted subset of special characters, all terminated 
by the (very) special character "*". All characters which may be used 
are mentioned explicitly in the following subsections of this section* 

A format specification, in order to be used, must be stored in 
successive elements of a vector of integer modej therefore, it is stored 
in groups of six characters per element. A format specification should 
be stored in such a way that the first character in the specification is 
the left-most of the six in a vector element. If the specification is 
stored in such a way that the first character is not the left-most in 
an element, then those characters to the left of it mxist be blanks . 

A vector may be preset with a format specification by the use of 
a VECTOR VALUES declaration (see section 3.7). Also, the elements of 
the vector may be computed or read in as data . 

5 .2 .1 Single Line Format Specifications 
A single line format specification has the form: 

T T . . . T * 
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where the T^'s are terms of the format specification, the character,"," 
separates pairs '..of terms; and. is included in the format specification, and 
the character "*" follows the last term of: the.^format sp^ificatios and 
must appear in the format specif icationj it serves to terminate the 
specification. 

Any term T^ may be one of the following: 

(a) a Basic Field Description (see section 5.2.11, Basic Field 
Descriptions - A ResumS); 

(b) A multiple Basic Field Description (see section 5-2.12, Multijjle 
Basic Field Descriptions) j 

(c) a Scaled Field Description (see section 5.2 .13, Scaled Field 
Descriptions). 

5.2.2 Multiple Line Format Specifications 

o — — ^4^^:^ .^^v-v,j.j. j.^i^ oj.vjxi majr uc luauc ou X cx ci' ou uiurc unan one 

line or card at the time it is used. 

A multiple line format specification has the form 

where each X^ may have the form of a complete single line format specifica- 
tion without a terminating "*" but with carriage control if required (see 
section 5-2.14.2, Carriage Control) or any X.. may be blank to indicate a 
blank line or card for output or an ignored card for input; the character 
"/" is used to terminate a single line format specification X within the 
multiple line format specification, and the character "*" is used to 
terminate the last single line format specification of the multiple line 
format specification. 

5.2.3 Format Fields 
Each format specification describes successive fields across 
the available columns, starting from the left. If the specif ica-;' . 
tion describes fewer than the total number of available columns, 
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the remainder of the line or card will be filled in with blanks . 
If J at execution time , a format specification is used which de'-^- 
scribes more than the total number of available columns^ an error 
indication will be given and the problem will be terminated. 

The seven types of fields which nay be described in a format 
specification appear in the following list: 

S-Field; skip or blank information 

I -Field J integer 

P -Field; fixed point number 

E-Field; floating point number 

N-Field; octal number 

G-Field; BCD characters 

H-Field; Hollerith information 

The terminology "fixed point number" used in connection with 
F-fields and "floating point number" used in connection with 
E-fields is rather unfortunate but is difficult to avoid. It is 
meant solely to provide a means of distinguishing "the form of 
the information which occupies an F field" from "the form of the 
information which occupies an E field" and throughout this section 
is used solely for that purpose. MAD does not handle "fixed point 
numbers" internally. 

"Format Field" is a name for two concepts taken together. 

First, there is the "Basic Field Description". This is an 
item which appears in a format specification. 

Second, there is the "Field Infonnation" . This refers to the 
information and the form of the information which may appear on 
a card or a printed line in the field described by the corresponding 
Basic Field Description. 

Both the Basic Field Descripton and the Field Information for 
each of the seven types of format fields are described in the 
following seven sections (section 5.2 .it- through 5.2. lO). 
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5.2.4 S-Fields 

5.2.-4.1 S-Field Basic Field Description 
The "basic field description for an S-field has the form: 

Sn 

where the character "S" indicates an S-field and n is a decimal integer 
equal to the niMiber of columns in the field. 

5.2 .i<-. 2 S-Field Information 

If an S-field basic field description appears in a format specifipation 
which is used for input, any information which appears; in the corresponding 
columns will he ignored. 

If an S-field is used for output the corresponding columns will be 
blank . 

Example: The basic field description 

S51 
indicates that thirty-one columns are to be skipped. 

5.2.5 I-Fields 

A number in an I-field has the form of an Integer internally. 

5.2.5.1 I-Field Basic Field Description 

The basic field description for an I-field has the form 

In 

where the character "I" indicates an I-field and n is a decimal integer 
eqiial to the number of columns in the field. 

5.2.5.2 I-Field Information 

The information in an I-field may have any of the following forms: 

+m 

-m 

m 
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where m is any integer satisfying 

< m < 9999999999 

If^ as in the third form, no sign is punched, the integer 
will be assumed to be positive. Any possible sign must be counted 

If an I-field is used for input: 

All blanks in the field are ignored; 

If the entire field is blank, the value will be set 
equal to minus zero; 

Leading zeros need not be punched j 

Trailing zeros miist be punched) 

If the integer is negative the "-" sign must be punched. 

If an I-field is used for output: 

For positive integers, "+" signs are not printed or 
punched; 

For negative integers ^ "-" signs are printed or punched; 

If the integer internally contains fewer digits than 
provided for by the field size, it will print right 
justified in the field with blanks in the remainder 
of the field; 

If the integer internally contains more digits than 
provided for by the field size, the least significant 
digits will be printed and the sign and the remainder 
of the most significant digits will not appear; 

Wo decimal point, ".", is printed or punched in an 
I-field . 

Example: The basic field description 

13 
will cause three columns of an input card to be treated as a 
decimal integer. The same basic field description will cause 
information printed or punched into some partic\ilar three columns 
to have integer representation. 
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5.2.6 F-Fields 

A number in an F-field has the form of a floating point number 
internally . 

5.2.6.1 F-Field Basic Field Description 

The basio field description for an F-field has the form 

Fn.k 

where the character "F" denotes an F-field, n is a decimal integer 
equal to the number of columns in the field, k is a decimal integer 
equal to the number of digits to the right of a missing decimal point, 
and the character "." (period) separates n and k. The integer k is 
always interpreted mod 10, e.g., k = 13 is equivalent to k = 3. 

5.2.6.2 F-Field Information 

The information may have any one of the following forms; 

+ra-A +m.^E+e +mjlE+e 

•^ m./jrfe mJ^E+e 

X. +ra.XEe +m/Ee 

ni -<-' m.^Ee mJ^Ee 

+m.j2.+e ±ini+e 

ra-H+e r&jL+e 

where ra is a decimal integer, jl is a decimal fraction, and 
e is a decimal integer equal to the e^onent of the power 
of 10 by which the number + m.|or +mjf is to be multiplied! 
e may contain one or two digits, no more. 

The character "E", if punched^ indicates that an exponent 
follows J 

If a sign ("+" or "-") does not appear as the left-most 
character, the number is assumed to be positive^ 

The sign which follows the character "E" is the sign of the 

exponent, ej 
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If no sign character follows "E" the exponent e will be ass-umed 
to be positive; 

Jlotice that the "E" is necessary only if the sign of the ex- 
ponent is not punched and that the sign of the exponent is 
not necessary only if the "E" is punched and the exponent is 
positive; 

All blanks are ignored; 

If an entire F-field is blank the value will be set equal to 
minus zero; 

If the F-field information is in one of the forms with a decimal 
point ("•")j this "." in the field information will override the 
effect of the number k in the basic field description. In that 
case, trailing zeros in X need not be punched; 

If the form of the field information is one for which the decimal 
point is not punched, then trailing zeros must hot be omitted; 

Any number of digits may be used in the field but only eight 
digits of precision are retained. 

In determining the field size, n, in. the basic field description 
the count must include any possible occurrences of the sign 
of the number, "+" or "-", a decimal point, ".", an "E", and 
an exponent sign "+" or "-" as well as maximum number of digits 
in m, jL, and e, combined . 

If an P-field is used for output: 

The printed or punched information will have one of the following 
forms : 



m 



.1 



-m. ,x^ 



m 
-m 



where m is a decimal integer and X is a decimal fraction rounded 
to k digits; 

The character "+" is not printed, the form' without sign represents 
a positive number; 

One of the two forms m or -m will occur when k=0 (the "." is 
not printed); 

The nimber will be right justified in the field in all cases; 

If the field size, n, is larger than required, for information, 
blanks will be printed or punched in the remaining columns to 
the left; 
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If the field size is smaller than required, information will 
be output from right to left until the field is exhaiistedj 
notice that a sign, "-", and/or a decimal point, ".", which 
otherwise would print or punch may he lost in this event; 

In determining the field size n the count must include possible 
occurrences of the characters "-" and "." as well as the 
maximum total number of digits in m and £ . 

Example I The basic field description 

F9.5 

will cause the following data all to be given the same machine representa- 
tion, namely that of 3.9962x10 ; 

+39962 
39962 
399.62 
3.9962E+2 
3-9962E2 
3996 .2E-1 
etc. 

The basic field description 

F9-3 

will cause the machine representation of 3.9962x10^ to be printed or 
punched as 

399-620 
right justified in the 9 columns. 

5.2.7 E-Fields 

A number in an E-field has the form of a: floating point, number ; 
internally. 

5.2.7.1 E-Field Basic Field Description 
The basic field description for an E-field has the form: 

En.k 

where the character "E" denotes an E-field^ n is a decimal integer equal 
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to the niJinber of columns in the field, k is a decimal integer equal to 
the number of digits to the right of a missing decimal point and the 
character "." separates n and k. The integer k is always interpreted 
mod 10^ e.g.j k = 22 is eq,uivalent to k = 2 . 

If ah E-field is used for input: 



The information 


may- 


have any one 


of 


the 


following forms: 


+m. .£. 




+m.£E+e 






+m£E+e 


m.X- 




m.J6E+e 






m I E+ e 


±mt 




+m . X'Ee 






+m X Ee 


m.-C- 




m .k Ee 

+m . }L±e 

m,X+e 






m.ilEe 

+m)L+e 

mj^+e 



where m is a decimal integer, JC is a decimal fraction, and 
e is a decimal integer eq.iial to the exponent of the power of 
10 which the number +m.)L or +m X is to be multiplied; e 
may contain one or two digits, no more; 

The character "E", if punched, indicates that an exponent 
follows; 

If a sign ("+" or "-") does not appear as the left-most 
character, the number is assumed to be positivej 

The sign which follows the character "E" is the sign of the 
exponent, e| 

If no sign character follows "E" the exponent e will be 
assimed to be positive; 

Notice that the "E" is necessary only if the sign of the 
exponent is not punched and that the sign of the exponent 
is not necessary only if the "E" is punched and the exponent 
is positive; 

All blanks are ignored; 

If an entire E-field is blank the value will be set equal to 
minus zero; 

If the E-field information is in one of the forms with a 
decimal point ("."), this "." in the field information will 
override the effect of the number k in the basic field 
description. In that case, trailing zeros in J^ need not be 
punched; 
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If the form of the field information is one for which the dee- 
imal' point' is not' punched, then trailing zeros- must not be 
omitted; 

Any number of digits may be used in the field but only eight 
digits of precision are retained. 

In determining the field size, n, in the basic field description 
the count must include any possible occurrences of the sign of 
the number, "+" or "-", a decimal point, ".", an "E", and an 
exponent sign "+" or "-" as well as maximum number of digits 
in ra, £ and e, combined. 

If an E-field is used for output: 

The printed or punched information will have one of the following 
forms : 

m .X. Eb'e, Bp mEbe, e 

-m .1 E-e-|_e2 mE-e, e 

-m.XEbe^eg -mEbe,e 

-lu. /v. xj-v_, Cp -ujii-e, e^ 

where m is a decimal integer, -^ is a decimal fraction rounded 
to k digits, "b" represents a blank and e-i_ and b are the two 
digits of the exponent or the power of 10 by which the m.£, or -m-X- 
is to be multiplied; e^^ is always punched or printed, even when 
equal to zeroj 

The four forms in the right hand column result when k = Oj 

The number will be right justified in the field in all cases; 

If the field size, n, is larger than required for the information, 
the remaining columns to the left will be: blank; 

If the field size is smaller than required, information will be 
output from right to left until the field is exhausted; notice 
that signs, "-", and/or decimal points "." which otherwise 
would print or punch may be lost in this event; 

In determining the field size, n, the count must include possible 
occurrences of the characters "-" and ".** and the two exponent 
digits as well as the maximum total number of digits in m and jl. 

Example: The basic field description 

E10,2 

will cause the following data all to have the same machine representation, 
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namely that of 2^^.629 

21^.629 
24629 -2-5 
24629-5 

The baaid field description 

E10.2 
will cause the machine representation of 24.629 to be printed or punched as 

246 .29-01 
right justified in the 10 columns. 

5.2.8 K-Fields 

^ere is a one-to-one correspondence between the 36 bits internally 
and the bits represented by the number in a K*field. Any octal digit 
represents 5 binary bits as follows: 

Octal Binary 

000 

1 001 

2 010 

3 Oil 

4 100 

5 101 

6 110 

7 111 

5.2.8.1 p:-Fleld Basic Field Description 
The basic field description for a K-field has the form 

Kn 

where the character "K" denotes a K-field and n is a decijnal integer 
equal to the number of columns in the field . 
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5.2.8.2 K-Field Information 
If a K-field is used for input: 

The information in a K-field may have any one of the following 
forms : 

+ P 

- P 

P 

where p is an octal integer satisfying 

< P < 7777777T7777J 

The octal integer p either contains or Implies a value for each 
of the 56 binary hits in a computer wordj the left-most of these 
36 bits coincides with the sign bit; -the sign of the number is 
determined by the .logical "or" of the: sign/ if any^: and the 
left-most of the 56 bits as given by the following list f"+" is 
0, "-" is l)t ■ 



Sign- in Field 


Left-Most Bit 


Sign of Number 


+ 





+ 


- 


1 


_ 


b 





+ 


+ 


1 


- 


- 





- 


b 


1 





All blanks in the field are Ignoredj 

If the. entire KAfield"' is" blank . the value will be set equal to 
plus zero J 

Leading zeros need not be punched; 

Trailing zeros must be punched; 

The field size must include possible occurrences of a sign 
chaDracter in the count. 

If a K-field is used for output: 

The information punched or printed will have the following form: 

P 
where p is an octal integer; 

5.15 
6/20/62 



- FORMAT SPECIFICATIONS 

No sign is punched or printed] the sign of the number is indicated 
by whether the left-most of 12 octal digits is less than k (+) or 
greater than or equal to h (■-)| 

The number is right justified in the fleldj 

If the field size is larger than required for the information the 
remaining columns to the left are filled with blanks; 

If the field size is smaller than required for the information 
the information is:, output from right to left until the field is 
exhausted and the most significant part does not appear. 

Example: The basic field description 

K5 

will cause the appropriate 5 columns of data card to be interpreted as 

containing three octal digits . For example the configuration 

will become the machine word 

000 000 000 000 000 000 000 000 000 000 100 Oil 

5.2.9 C-Fields 
The Hollerith card code for each character corresponds to the internal 
6-bit BCD code for the same character. 

5.2.9-1 C-Field Basic Field Description 
The basic field description for the C-field has the form: 

Cn 

where the character "C" denotes a C-field and n is a decimal integer equal 
to the number of columns in the field. 

5.2.9-2 C-Field Information 
The information in a Cwfieid tea the fowni 



where c is any string of Hollerith characters available on the computer 
(see section 1.1, Tbe Hollerith Character Set). Blanks are not ignored. 
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If a C-field is used for input: 

Characters are taken from the field from left to right either 
until six columns have heen read or until the field has been 
exhausted whichever occurs first j 

The information is left justified internally; 

If the field size, n, is less than six, the right-most 6-n 
characters are filled with blanks; 

If the field size n is greater than six, the right-most n-6 
characters are lost; 

If a C-field is used for output: 

Characters are output from left to right either until six have 
been transmitted or until the field has been exhausted, which- 
ever comes first; 

The information is left justified in the field; 

If the field size, n, is greater than six, the right-most n-6 
characters are filled in with blanks; 

If the field size, n, is less than six the right-most 6-n 
characters are lost . 

Examples : 

(1) If the format specification 

2C5* 

were used to read a card punched with 

ABCDEFGHIJK 

in columns 1 through 11, the two computer words involved would 
contain 

ABCbbb 

DEFbbb 

(2) If the format specification 

C6* 

were used to read the card of example (l), the single computer 
word involved woiild contain 



ABCDEF 
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(5) If the format specification 

were used to read the card of example (l), the two computer 
words involved would contain 

ABCDEF 

HIJbbb 

5.2.10 H-Fields 
Ihe purpose of the H-Field is to permit the inclusion of strings of 
Hollerith characters (see section 1.1, The Hollerith Character Set) directly 
in the foriDat specification itself. 

5.2.10.1 H-Field Basic Field Description 

The basic field description for an H-field has the form 

nHc 

where the character "H" denotes an H-field, c is any string of Hollerith 
characters and n is a decimal integer equal to the number of characters 
in the. string c counting from the character immediately following the "H" . 
Although, for every other type of fieW> the basic field . description.' 
which appears in a format specification must be separated from a following 
basic field description by a character ",", for the H-fielda, since the 
number n explicitly defines the nimiber of characters in the basic field 
description, the character "|" after the last character in the string c 
is optional . 

5.2.10.2 H-Field Information 

If an H-field appears in a format specification which i$ used for 
input, the Hollerith characters which appear in the corresponding card 
colianns will replace the n characters in the string c in the format 
specification itself. 
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If an H-field is used for output, the n characters of the string c 
will be printed or punched in the corresponding columns . The characters 
of the string c in the format specification are not disturbed. 

Example: The basic Field Description 

3HABC 

or alternately, 

3H 

vill cause, upon reading of a data card, the characters "ABC" in the field 
description to be replaced by the Hollerith characters from the appropriate 
three columns of the card • 

The basic field description 

5HABG 

will cause printing, or punching, of the characters "ABC" on the output 
in the appropriate coli:imns . 

5.2.11 Basic Field Descriptions - A Resume 

The purpose of this section is merely to gather together all in one 
place a recapitulation of the possible seven types.of basic field descriptions. 
They are: 

Sn (see section 5'2.^.l) 

In (see section 5'2-5-l) 

Fn.k (see section 5'2.6.l) 

En.k (see section 5-2.7'l) 

Kn (see section 5 -2 .8.1) 

Cn (see section 5-2 .9-1) 

nHc (see section 5 •2.10.1) 

5.2.12 Multiple Basic Field Descriptions 

If several consecutive fields can be described by the same basic field 

description, repetition may be avoided by using a multiple basic field 

description. 

5.19 
6/20/62 



FOEMAT SPECIFICATIOKS 

(a) a multiple basic field description may have the form 

ID 

where D is any one of the basic field description forms of 
section 5.2 .10^ Basic Field Descriptions - A Resuafe, and where 
i is a decimal integer eqixal to the nianber of consecutive fields 
with form Dj 

(b) if M is a basic field description of the form 

D 

or a multiple basic field description of the form 

iD 

or a scaled field description Sj where S may have the form of 
either a scaled basic field description (see section 5.2.15.I) 
or a scaled multiple basic field description (see section 5. 2. 15-2 ), 
then a multiple basic field may have the form 

j(Mj^, ^^, ..., M^) 

where j is a decimal integer equal to the niaraber of times the 
group of q field descriptions NL, VL, ..., U is to be repeated 

aiiu. niiSi c oiic >-j.ja.4.av-L.cj.o { iaxiU J UCXXlUjL L> bXlti gTOUp Oi I leiCL 

descriptions on the left and right respectively. 

(c) a multiple basic field description may have no other form than 
those given by (a) and (b) above. 

Notice in particular that parenthesized groups may not be nested. 

Examples ; 

(1) the format specification 

5F10.5,El8.4,2E9 .1,512* 
is a short way of writing the format specification 

F10.5,F10.5,F10.3,El8.i)-,E9.1,E9.1jI2, 12^12* 

They are logically equivalent. 
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(2 ) The format specification 
3E10.3,2(I2,3F10.1),2C5* 
is logically equivalent to the format specification 

ElO .5,E10 .5,E10 .3,I2,F10 .I.FIO .l,F10;.l, I2,F10 .1,F10 .1,F10 .1,C5,C5* 

5.2.13 Scaled Field Descriptions 
It is possible to write an F-field or an E-field description with a 
scale factor which will be applied to every number to which the description 
is applied. 

5.2.13.1 Scaled Basic Field Descriptions 
A scale factor may be applied to a basic field description. Such a 
scaled basic field description may have any one of the following forms: 

sPD 
+ sH) 
- sPD 

where D is either an F-field basic field description (see section 5.2 .6.1) 
or an E-field basic field description (see section 5.2.7 -1)^ the character 
"P" (for "power") denotes a scale factor and s is a decimal integer equal 
to the exponent of the power of 10 times which the number is to be multiplied . 
If the exponent is negative the sign "-" must be punched; if the exponent 
is positive J the sign "+" is optional. 

5.2.13*2 Scaled Multiple Basic Field Descriptions 
A scale factor may be applied to a restricted form of a multiple basic 
field description. Such a scaled multiple basic field description may 
have any one of the following forms 

sPiD 
+ sPiD 

- sPiD 

where D is either an F-field basic field description (see section 5 -2 .6.1) 
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or an E-fleld basic field description (see section ^.2.^.1), i is a decimal 
integer equal to the number of consecutive fields described by D (so that 

iD 
alone has the form of the primitive multiple basic field description of 

a decimal integer equal to the exponent of the power of 10 times which 
the number in each of the consecutive i fields is to be multiplied. If 
the exponent is negative, the sign "*•" must be punched] if the exponent 
is positive, the sigii ''4-" is optional. 

5.2.15.3 Scaled F-Fields 
When a scale factor is applied to an F-field the following formula 
is true: 

EN = IREN * 10 ^"^^^ ^^^*°^ 

where EN stands for "External Number" and IREN stands for ^Internal Re- 
presentation of External , Number " . The scaling (multiplication by 
10 ) is done before tJie. number is converted for output and 

after the number has been read and converted for input. Notice that 
scaling actually changes the value of the number in an F-field. 

Example : 

Suppose the format specification 

3F7.3* 

is used and as a result the three numbers which follow are printed s 

^ bbO.522b-l.567b93.67i 

If, instead, the format specification 

-2P2F7.3,F7.3* 

were used, the same numbers would print as: 

bbO .005b-0 .Ol6b93 .67I 
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5.2 .15 A Scaled E-Fields 
When a scale factor is applied to an E-field the basic field de- • 
scription may he used only for output. Although the number is modified 
the exponent is also modified so that the value is unchanged^ only the 
form in which the number is printed is changed. 

Example : 

Suppose that the format specification 

El8.i+* 

is used and as a result the niimber 

bbbbbbb0.9321E-03 

is printed. 

If the format specification 

2PE18.4* 

were used instead, the same number would print as 

bbbbbbb95 .2100E-05 

5 .2 .ik Format Specifications and Reading, Punching and Printing 
With the following two exceptions, specifications for reading or 
punching cards and printing lines are identical. 

5. 2. 1^1-. 1 Available Colianns 
(a ) Cards 

Whether reading or punching, the maximum number of card columns 
is 80. A format specification used for cards may describe 80 
columns. If more than 80 columns are described at execution 
time indication of an error is given and the job is terminated. 
It is an often useful convention that the card columns used 
for data be limited to 72, leaving the remaining columns for 
purposes of identification. Machine configurations at some 
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(b) Priirt Lines 

As in the case of cards j, the niimber of characters available for 
a print line depends on the equipment beijag used. When printing 
on-line 119 columns are available. Otherwise (off-line) 132 
columns are available. On-line printing is to be used extremely 
rarely. If more than 132 columns are described at execution 
time, indication of an error is given and the problem is 
terminated . 

5.2.lil--2 Carriage Control 

(a) The character destined for the first column of a card, when 
pimchingji has no special significance. It is regarded in the 
same light as the remaining 79 columns . It is punched in coliamn 
1 and the successive characters go to successive coliomns . 

(b) The left*most character output for a print line is the "carriage 
control character . " It controls the printer carriage "vertical 
motion," immediately before the line is printed, i.e., "controls 
the preprint skiiTt" The code for the carriage control character 
is given by the following table; 

If the left-most The preprint 



blank single space 

d double space 

•^ no space 

'" triple space 

1 sheet eject (skip to next page) 

^ skip to next half page 

k skip to next quarter page 

6 skip to next sixth page 

This control character is specified in a format specification in 

addition to the 119 or 132 information positions . 
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In any format specification which is used to print a line, the left- 
most character output is detached to be used as carriage control and is 
not printed . The appearance of the line which is printed is as if the 
carriage control character had gone to "print position zero" which is non- 
existent. The second character is printed in print position 1, the third 
in print position 2, etc. 

Note that "blanks count as characters . 
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5.3 The Input-Output List 

The purpose of the input-output list is to specify sources or des- 
tinations of the information transmitted. An input-output list, which 
will "be referred to as^, has the form 

where the T. 's are the elements of the list. The elements of an l/o list 
may be: 

(1) simple variables or subscripted array variables j 

(2) blocks, which may be written in one of the two forms: 

(a) V(s^, Sg, ..., s^) ... V(r-|^, r^, ..., r^) 

where V(s^, s^, ..., s^) and V(r^, r^,, ..,, r^) are 
subscripted array variables with the s.'s and r. 's any 
integer valued arithaiBtic expressions and where the 
characters "..." indicate that all consecutive elements 
of the n-dimejisional array from V(s^, s , ..., s ) through 
V(r,, r , ..., r ) inclusive are to be transmitted; 

(b) V(st^, Sg, ..., s^^), ..., V(r., r„, ..., r ) where commas 
(",") appear between the end elements in accordance with 
usual mathematical notation. The appearance of the 
commas is the only difference between (a) and (b) and the 
meafeings are identical . 

If the block specification is used in an Input-Output list for reading 
or writing binary tape , the linear subscript corresponding to s, ^ s , ..., s 
may not be greater than the linear subscript corresponding to r, , r„, .,., r . 
The elements of the list are transmitted from left to right, but, ^or binary 
tapes , within a block the^ order- 'is ;V(r, . r., ..., r ) to V(s, , s^, ..., s ), 
i.e., right to left. 

In the particular case when a list is used for output only and not 
input, additional elements of the list may be: 

(5) any constants j 

(k) any expressions. 
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Example of a list which may be used either for input or output: 

A,B,K(5),J(25*I-L),A(K+1) . . .A(L*2) 
Example of a list which may he used for. output only: 
AB,r)2.5,MTY(l), . . .,MTX(n),P(i1+),J(I,K) 
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5 A Relationship Between the 10 List and the Format Specification 
The "list" consists of a set of names of variables to or from which 
information is to flow. -• Except for Hollerith strings imbedded in the 
format specification itself (see section 5-2»10) and the S fields, each 
field in the specification refers to one item on the list . For this 
purpose, a regional entry on the list, such as A(6)...A(8), counts as 
several names of variables (in this case, the three variables A(6),A(7)j 
and X(8^). During the transmission of information, the input or putput 
subroutine scans both the list and the specification simultaneously, 
correlating corresponding entries, and associating a field size, a type 
of conversion, etc., to each variable. If a Hollerith string is encountered 
in the specification, it is immediately transmitted, and it is not associated 
with any item on the "list " . 

For example, if the list consisted of?. 

A, B(l,l), I, K 

where I and K were integers, and the others floating point, and the 
specification were 

IHl, F11.5, -2PEl4»i^, S15, aHM=I5^ S9, 2HJ=I3* 

we might find a printed line like the following (at the top of the next 
page because of the lEl); 

k^6 .010bb-l6 .1251E+02bbbbbbbbbbbbbM=50bbbbbbbbbJ=17 

The same list would look as follows with the format specification 

IHl, 2F11.3, Sl6, 2HM=I3, S9, 2HJ=T3*- 

i|-56 .010bb-l6l2 .510bbbbbbbbbbbbbbbbM=50bbbbbbbbbJ=17 

As stated above, a specification may not account for more than 80 
columns on a card. It may happen^ however, that a list calls for more 
information than can appear on a single card.. Or perhaps only a certain 
part of each card is to be read. The determining factor in every case 
is whether or not the entire list has been accounted for. After each 
card is read according to the format specification, the list is consulted^ 
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if it is not yet satisfied, another card is read, and so on. It is important 
to realize that the specification is not necessarily scanned from the 
beginning when a new card is read. In fact, the specification scanner moves 
to the left from the end of the specification (the *) until it hits a left 
parenthesis not in an H field, (if there is no left parenthesis, it will 
move to the beginning of the specification). It then examines the characters 
just to the left of this left parenthesis to see if they are a multiplicity 
indication (see section 5.2.2). The information from this left parenthesis 
(together with the multiplicity, if any) to the end of the specification now 
becomes the format specification until the list is satisfied. A similar state- 
ment may be made for printed or punched output . 
Thus, if the specification which follows 

5F10.5 A(F10.3, 6HBETA = 12)* 

is used with an 1/0 list which contains more than 11 elements (say 19 or 
27), then the first line printed (or read) would have three fixed point 
numbers. Subsequent lines would all be printed (or read) according to 
the specification 

ij-(F10.3, 6HBETA = 12)* 

until the last element of the l/0 list were transmitted . 

As another example one might have an integer equal to a count of sub- 
sequent cards on a first data card, followed by many cards, each with six 
floating point numbers. The specification might then be: 

i6/(6eio.5)*. 

Only the first six columns would be read on the first card, and only 60 

columns would be read on subsequent cards . The remaining columns are 

ignored and may contain any legitimate H9llerith characters . 

If a specification contains a Hollerith string of the form nHa,a„...a . 

1 2 n' 

certain conventions are observed: 

(1) If the list is satisfied, but the next field specification is 

a Hollerith strixig, the sLrixig is transmitted anywayj 
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(2) On input, i.e., reading from cards, when a Hollerith string is 
enconntered in the specification, the information in the 
corresponding columns of the input card will be "brought in and 
will replace the Hollerith s.tring itself within the format 
specification. This can then he used as a specification for out- 
put. ■For-eKam^le, this is useful for labeling a set of data 
and causing the label to appear on the output along with a date, etc. 

Thus J a card punched as follows; 

1 DATA SET NO. 5~A JULY 51, 1959 J DOE 

might be read in with a format specification 

72H (72 blank spaces) *. 

Later, this specification co\ild be used to print the same information as a 
heading for the results. Note that "1" provided for carriage control for 
the printing. 

WARNING; The specif icationseST2* and 72 (1H)*, while indicating 72 blank 
spaces, do not allow the reading in of an entire card, as indicated above, 
since they do not provide a storage region of 72 characters in length into 
which the information on the card may be read and stored until needed . 
See the example in section 3.7, VECTj^R VALUES Declaration. 

5.^.1 Input and Output of Boolean Values 
For the piorposes of input and output the Boolean values OB and IB are 
considered as integers and an I-field should be used to transmit them. 
Externally these values have the form of the digits and 1 respectively. 

5-^ -2 Input and Output of Statement Label Values 
Output of statement label constants may be accomplished by means of 
a C-field or an H-field. Input and outputof variable values may be 
accomplished by means of K-fields but involves concepts beyond the scope 
of this manual. This will be of no use to most users and should not be 
attempted without a deeper investigation of MA.D. 
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5.^.5 Input and Output of Function Name Values 
Output of fimction name constants may be accomplished by means of a 
C-field or an H-field. Input and output of variable values may be accomplished 
by means of K-fields but involvei-concepts beyond the scope of this manual. 
This will be of no use to most users and should not be attempted without 
a deeper investigation of MAD. 
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5-5 Basic Input-Output Statements 

Of these, the input statements are used in order to read any data 
which may appear on cards following the user's program in which these 
statements appear. The output statements are used in order to cause 
results calculated when the object program is executed to he printed on 
paper (listed ) or punched on cards . 

Throughout the following sections the following notation will be used: 
F: F is the name of the first element of a format specification "Vector 
which must be of integer mode. F may be written as 

(a) an array name, 

(b) a subscripted variable^ 

(c) a simple variable in case the format specification is less than 
or equal to six characters. 

See section ^.2, Format Specifications. 

<^iC is an input-output list, see section 5-5, 'The Input-Output List. 
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5.5.1 READ FORMAT Statement 
This statement has the form 

READ Fj^FvMAT F, X 

This statement causes data in the users data deck to be read into 
list <>Cin accordance with the format a't F. If the purpose is merely to 
read characters into a format specification comprised entirely of an H-field, 
in which case the list will be empty, the statement may be written 

READ F0RMAT F 

5.5-2 PUNCH F|6rmAT Statement 
This statement has the form 



PUNCH F0RMAT F,^ 



This statement causes data in the listX^.'^o '^s punched onto cards 

in accordance with the format at F. If the purpose is to punch using a 

format specification comprised entirely of an H-field, in which case the 

list will be empty, the statement may be written 

PUNCH Y0Tm.T F 
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5.5.3 PRINT f^RMT Statement 
This statement has the form 

PRIFT F0RMT F, ^ 

This statement caxises data in the list<5Lto be printed in accordance 
with the format at F. If the purpose is to print a comment or heading 
line using a format specification comprised entire of an H-field, in 
which case the list will be empty, the statement may be written 

PRINT FORMAT F 
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5.6 FEINT ^N LINE F0EMA.T Statement 
This statement has the form 

PRINT 0N LINE FORMAT ^ ,^ 

The listoC'^ls printed on-line in accordance with the format at-F. .If 






tne purpose xs t,0 prino a comment uaxiiy a j-ui-iijao B^c;i;xj.xt;aijj.uii ijutuyixocu 
entirely of an H-field, in which case the list will be empty, the statement 
may be written 

PRINT 0N LIKE FORMAT F 

After <i has been printed, a skip of one-sixth page is produced to allow 
the operator to read the comment . 

This statement provides the facility for printing comments to the 
computer operator on the on-line printer at execution time. The usage is 
quite restricted! and th^ statement will not bS'"used by most users- 
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5.7 Auxiliary Tape Storage Statements 

In this section J the following notation will be used: 

: N is an integer valued expression. The value will be used as a tape 
number and must be in the range [l, 10]'. To determine which numbers 
in particular may be used, see the manual describing the operating 
system in which MAD is imbedded. 

5.7.1 READ BCD TAPE Statement 
This statement has the form: 

READ BCD TAPE N, F, <^ 

This causes the listtjCto be read from tape N in BCD mode in accordance 
with the format at F. 

5.7-2 WRITE BCD TAPE Statement 
This statement has the form: 

WRITE BCD TAPE N, F, ^ 

This causes the list oC to be written onto tape N in BCD mode in 
accordance with the format at P. 

5.7.3 READ BIHARY TAPE Statement 
This statement has the form: 

READ BIIARY TAPE H, cS-^ 

This causes consecutive words from the next record on tape N to be 
read into the list ©C in binary mode. Transmission stops when the list 
is exhausted or when the end of record is reached, whichever occurs first. 

5.7.^ WRITE BINARY TAPE Statement 
This statement has the form: 

WRITE BINARY TAPE N, £ 

This causes the list oC to be transmitted as consecutive words com- 
prising the next (l) record on tape N, in binary mode. 
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5.7.5 REtflMD TAEE Statement 
This statement has the form: 

REWIiro TAPE N 

This causes tape N to he rewound, i.e., positioned at load point. 

5.7.6 EKD 0F FILE TAPE Statement 
This statement has the forms 

END 0F FILE TAPE N 

This causes an end-of-file mark to be written on tape N. 

5.7.7 BACKSPACE RECORD pF TAPE Statement 
This statement has the form: 

BACKSPACE BECKED f>Y TAPE N 

This causes tape N to be moved backwards to the beginning of the 
preceding record . 

5.7.8 BACKSPACE FILE Statement 
This statement has two -forms: 

(1 ) BACKSPACE mm: 0F mPE : If ^'I": ■ 

This causes tape N to be moved backwards until an end-of-file mark, 
the load point gap or the load point is encountered. If it is an end-of-file 
mark which is encountered, this statement will cause the tape to be 
positioned Immediately preceding the end-of-file mark. 

(2) BACKSPACE FILE ^F TAPE K, IF. Lj^AD P0INT TRANSFE?^ T^ S 

The action for this statement is exactly as the preceding one with 
the additional facility that, if the tape is already at load point, the 
program transfers to the next executable statement given by the; 'Statement 
label expression S. 
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5 •7-9 Action of End-of-File and End-of-Tape 
If an error (improperly formed format specification, invalid data, 
a tape check, etc.) occurs during any input-output statement, the system 
subroutine ERE0R. is automatically entered. The subroutine ERR0R. sets an 
error flag and returns control to the system in which the translated program 
is imbedded . 

If an end-of-file is encountered while executing a READ F^EMAI, READ 
BCD TAPE, or READ BIMARY TAPE statement the subroutine SYSTEM, (a subroutine 
which returns control to the system in which the translated program is 
imbedded) is automatically entered. This action can be changed by executing 
the statement 

EKECUTE SETEJ^F. ^(S) 

The subroutine SETEJ^Jf. sets the read routines to transfer to the executable 
statement labeled S when an end-of-file is encountered. If the statement 
EXECUTE SETEjz5f. (o) is executed the read routines will be reset to enter 
SYSTEM, when an end-of-file is encoimtered. 

If an end-of-tape is encountered while executing a WRITE BCD TAPE or 
WRITE BINARY TAPE statement, no special action is taken. If writing continues, 
the tape may run off the reel . This can be changed by the executing the 
statement 

EXECUTE SETETT. (S) 

The subroutine SETETT. sets the write tape routines to transfer to the 
executable statement labeled S when an end-of-tape is encountered. Executing 
EXECUTE SETETT. (o) will reset the write tape routines to the normal situation. 

The subroutines SETE0F. and SETETT. may be executed as many times as 
desired. Only one setting is in effect for end-of-file (that specified by 
the latest execution of SETE^F.) and end-of-tape (that specified by the 
latest execution of SETETT.), i.e., each setting cancels the previous one. 
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6 Function Definitions 

In this section the executable statements and declarations necessary 
to define function subprograms are described; for use of defined functions 
see secsc 2. 4»1> 'Single Valued Functions, 2. 4« 2, Non-fSingle Valued Functions, 
and Z'k'5} Arguments of Functions* 

A "function definition" is a sequence of executable statenents from 
section 4 and declarations froa section 3 qualified by and delimited by 
the declarations and executable statements of sec. 6.3, Function Definition 
Statements, below. Further restrictions are imposed upon the sequence of 
statements by the considerations given in sec. 6.4, Internal Function 
Definitions, and sec. 6.5, External Fimction Definitions. 

As has already been alluded to, there are two main types of programmer 
defined functions: the internal function and the external function. 
Either of those types may include single valued functions or non-single 
valued functions (non-single valued functions are sometimes referred to as 
"procedures"). 

A "recursive function" of either type may also be defined. For the 
purposes of this manual a recursive function is a function whose definition 
calls for the function being defined or calls for a fvmction which 
ultimately calls for the function being defined. 

The general structure of the definition of a recursive function is the 
same as for any other function but it will include seme of the statements 
described in section 7, List Manipulation Statements. 

When defining a recursive function, it is necessary to keep in mind 
that names are used as function arguments and not values (see section 15, 
Examples, for fiirther considerations with regard to recursive functions). 

The name of a defined function must be distinct from the nanes of 
functions already available to the translator. For a list of these 
names see the manual of library subroutines available. 

Each function definition (except Single Statement INTERNAL FUNCTI0N 
Definitions, sec. 3«8) may define any number of functions and/or any 
number of procedures. 
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60 1 Dummy Argiaaents 

Scane functions may be defined which require no arguments-. Often, 
however, it is convenient to be able to specify which values are to be used 
during the evaluation of a fianction, i^e^j to specify arguments, A dis- 
cussion of the form of the arguments when a function is used appears in 
sees. 2' 4:1s. Single Valued Functionsj 2 = 4, 2, Non-Single Valued Functions j 
and 2.4=5 J Argiments of Functions » 

We are concerned here with the form of the argisnents when the function 
is defined s these are called "dummy arguments". 

A dummy argument is used in the statements of the function definition 
to indicate a correspondence between the variable name in the statement 
and a position in the list of arguments. When the function is used every 
occurrence of the "dimaiiy variable"^ say A. 5 in the statements of the 
definition will be replaced by the value of the expression in the position 
of A. in the list of dismii'"" ar^'umentss 

The fomi of a dtrnmy argument may bet 
(i) a name of a simple variable, 
(2^ a name of an array variable, 
(3' a name of a function.-, 

riotice that neither constants nor subscripted variables may appear as 
dummy arguments and that arguments of a function do not appear following 
the function name in the dummy argonent list. 

Names which appear as dvanmy arguments in any function declaration may 
not appear in a PR0GRAM C0»M0N, ERASABLE or EQUIVALENCE declaration- 

The modes of dummy arguments must be declared as for other variables 
if of other than normal modeo 

Dianmy arguments which are array names need not be dimensionedo 
Examples The list 

( SIMPLE, FCNo, ARRAY) 

is a valid dianmy argument listo The following list is nots 

(FCNo(X,Z+2o,VARl), ARRA.Y(l,9)) 
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6.. 2 The Use of a Defined Function 

In the use of a function, either in an expression or in an EXECUTE 
statement, the arguments may be constants, variables or expressions. 
However, if one of the dummy arguments in the definition appears to the 

statanents, then it is not meaningful to use a constant or an expression 
in the position of that argument in the call for the function. The 
arguments in t]ie call for a function cannot be checked by the translator 
for correspondence in mode and number to the arguments in the definition. 
It is left to the user to see that they agree. 
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63 Function Definition Statements 

631 INTERNAL FUNCTI0W Declaration 
This declaration indicates that an internal fvinction is to be defined 

by,: feh« statements which follow.. This declaration has the form; 

INTERNAL FUNCTI0N. (A^, A^s - - s A^J 

where the A. are the disBmy argianents for the function ■-• Note that n© 
function name appears in the declaration= This statement must necessarily 
appear as the first statement of an internal function definition. Section 
6»45i Internal Function Definitions, contains further discussion regarding 
the definition of internal function&o 

6.3=2 EXTERNAL FUNCTI0N Declaration 
This declaration indicates that an external function is to be defined 
by the statements which follow. This declaration has the form: 

EXTERNAL FUNCTI0N (A, , A_, » o » , A ) 

X < n 

where the A, are the dummy arguments for the function. Note th*t no 

function name appears in the declaration^ This statement must necessarily 

appear as the first statement of an external function definition Section 

6:5j External Function Definitions, contains further discussion regarding 

the definition of external functions = 

6o3»3 ENTRY T0 Declaration 
This declaration has the form: 

ENTRY T0 F= 

where F.. is a name of a function which is to be defined by this function 
definitioHo This declaration marks the entry to the function F. = Execution 
of the function F= begins with the first executable statement following 
the ENTRY T0 F- declaration: 

In a single program a given function name may appear in only one ENTRY 
T0 decli rationi also a given function name may appear in only one ENTRY T0 
declaration within a given function definition. However, a single function 
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definition may contain many ENTRY T0 declarations each vdth a distinct name 
Fo , to d esignate entry pointg to several function naijies within the function 
definition « 

6.3 » 4 FUNCTI0N RETURN Statement 

This is an executable statement which terminates the execution of the 
function and returns to the place in the program where the function name 
was used in such a way that computation may precede using the result or 
results of the function^ 

There are two forms of the FUNCTI0N RETUR.M statement: 

(1) FUNCTIjZ&i RETURN E 

where E is any expression, the mode of which is consistent with 
the mode of the function being evaluated. At the time of return 
the value of the function will be made the value of E= The mode 
of the expression E is not checked by the translator^ If a 
fimction name is to be used in any expression, i=ec, if it is to 
be a "single valued function" (see sec. 2o4.1j Single Valued 
Functions), then the FUNCTION RETURN statanent which teminates 
the execution of the function must have this form» 
Example: The statement 

causes the value of the function in which it appears to be set 
equal to the expression 

BETA/K5-4.*D 

(2) FUNCTI0N RETURN 

This form of the FUNCTION RETURN statement, which does not include 
an expression, is used to terminate the execution of siny function 
which is not a "single valued function" in the sense of this 
manual, i.e., any function which is not used in an expresbion.. 
It should be noted that a f iinction may have a single result » such as 

setting a single variable equal to a calculated value, and still not be a 

"single valued function" in the sense of this manualo 
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On the other hand, a function may be defined which yields several 
results, such as setting several variables equal to several calculated 
values J but which may nevertheless be a single valued function if it 
terminates execution by a FUNCTI^ EETURN statement of the fom under (1) 
above : 

A single FUNCTION RETURN statement may be used by many functions 
within the function definition » 
Examples? , 

(1) The internal function definition 

- ' INTERNAL FUNCTI0N (X,Y,Z) 
ENTRY T0 FCN» 

X = SINo(Y) oP. 2 + SIN.(Z) oP. 2 
FUNCTION RETURN 
END 0F FUNCTI0N 

has a single result? the value of X is set equal to the value of the 
expression 

sin^(Y) + sin^(Z) 

Nevertheless, this is not a "single valued function" since the FUNCTI0N 
RETURN statement used contains no expression . 

(2) The internal function definition 

INT131NAL FUNCTI0N (X,Y,Z) 

ENTRY T0 FCN2. 

X = SINo(Y) 

Y = C|iS.(Z) 

Z = ATANo(X) 

FUNCTI0N RETURN X + Y + Z 

END 0F FUNCTI0N 

is a "single valued f-unction" whose value is equal to the value of the 
expression 

X + Y + Z 

even tho.ighj in addiction, new values for all of X^ Y and Z are calculated:- 
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6o3o5 ERR0R RETURN Statsnent 
This statanent has the form: 

ERR0R RETURN 

It is an executable statement which may be used in a function defin- 
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If the right-most argument in the call of the function in the defin- 
ition of which the ERE0R RETURN statement appears is an expression of 
statement label mode (see seco 2»10, Statement Label Expressions), and if 
the ERR0R RETURN statement is executed, then execution of the function is 
terminated and control is transferred to the statanent the label of which 
is equal to the value of the aforementioned statement label expressiono 

If the ERR0R RETURN is used in a function definition and a call for 
the function is made such that the right-most argument is not of statement 
label mode, and if the ERR0R RETURN statement is executed, control will be 
returned to the system and the execution of the program will be terminatedc 
An error indication will be printed= 

6o3o6 END 0F FUNCTI0J Declaration 
This declaration has the form: 

END 0F FUNCTI0J 

It delimits the extent of any function definition and hence must 
appear as the physically last statanent in any function definition (with 
the exception of the Single Statement INIEHKAL FUNCTION Definition of 
section 3 08)0 It is not an executable statement and the execution of the 
function may not be tenninated by attempting to execute ito If it is 
executed an error indication will be given and the problon will be termi- 
nated at execution times 
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60 4 Internal Function Definitions 

An internal function definition is a function definition which is to 
be translated as part of the main program c 

Names used as dummy arguments must be distinct from all other names in 
the pro gran o 

The name of the defined function must be distinct from all other nanes 
used in the pro gram =. 

Names of variables or ftmctions which do not appear as dummy arguments 
but which are defined elsewhere in the program may also appear in the 
internal function definition statements = The value of such a variable or 
function is the current value at the time the internal function is used. 

Internal function definitions of all kinds (including single stat«nent 
definitions)^ may occur anywhere in the program, except within another 
internal function definition Ihtei^nal function definitions may occur 
within external function definitions : 

Examples The following is an internal function definition of the function 
whose name is NAMEFo J 

INTERNAL FUNCTI0N (P,DjFc) 

ENTRY T0 NAMEFo 

TMP = F 

P = F.(D) 

D = Fo(TlMP) 

FUNCTI0N RETURN 

END 0F FUNCTI0N 

This is a single valued function whose value is a constant integer zeroc 
In addition the function sets P equal to F-(D) and D = F (P) where F= is 
any defined function name (other than NAMEF. ) used in a call for NMEFc = 
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60 5 External Function Definitions 

An external function definition is a f-onction definition which is to 
be translated completely independently from the main program in which it is 
to be usedo The external function definition as a whole appears completely 
outside any other programo Because of this an external function definition 
is an entirely separate^ complete program.. Thus it must contain its own 
DIMENSI0N declarations, mode declarations, FI^RAtf C0M0U declarations, 
ERASABLE declarations and any others which it itself requires. However,, 
it is teiminated by an END 0F FUNCTI0N statement and not by an END 0F 
PR0GRAM statement as are other programs o 

Declarations required on the dvsm.j arguments appear within the external 
function definition.. 

Since an external function definition is a complete program there 3,3 
no conflict between names used within it and names used within a main 
program using it or names used within another external function used by the 
same main programo 

Examples The following is an example of an external function vAiose value 

2 
is 1/x if < X < 1 and 1/x if x > lo If x < 0. one obtains an error 

return (see section 6c 3=5, ERE0R RETURN statement )o 

A EXTERNAL FUNCTI0N (X) 



J 


ENTRY T^ irvSP- 




G 


WHENEVER X=GoO. oANDo 


X oLE 


C 


FUNCTION RETURN X 


e.P= -1 


H 


0R WHENEVER X .G: 1. 




D 


FUNCTI(3^ RETURN X 


.P. -2 


T 


0THERWISE 




E 


ERR0R RETURN 




K 


END 0F C0NDITI0NAL 




B 


END 0F FUNCTIC^ 





Here the statements are all labelled only for reference in what follows =. 

This definition program defines a single-valued function of X, called 
INVSF=o Since no mode declaration is given it is assumed by the translator 
that X is floating point even though a different normal code were declared 
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in the main progranic. The value of INVSF: (X) is computed by the use of a 
compound conditionals If < x < Ij (statement G) then statement G is 
executed J causing a return to the calling program with th-? value ~ = If ^ 
the condition < x < 1 is not true^ then the condition jt > 1 is tested 

(statement H)» If ^ > 1, statement D is executed: Finally, if neither of 
the conditions 0<x<lorx>lis true^ then statement I finds that 
X < and stat«nent E (the error return) is executedo If 

A = B - D 

Z = T(I) + INVSFo(Y) * T(l-l) 

Y(I) = Z + R(J) « 2o5 

is part of a program and the error return statement is executed dicing the 
evaluation of INVSFo(Y) (io8„, Y > O), then control is returned to the 
system in which the translated program is imbedded, with an error flag set. 
If 

A = B - D 
F Z = T(I) + INVSF=(YjE|l) * T(l-l) 
S Y(I) = Z + R(J} * 2.5 



ER Z = 
L Y(I) = 1. 

o o o 

is part of a program and Y < Oj then when the EHR0R RETURN statement is 
executed control transfers to statement ER (then goes on to L)^ instead of 
finishing the execution of statement F (and then going to S)o Note that 
the END 0F FUNCTION statement will never be executed, but must be present in 
the definition. 
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7 jLiXSw MampUj-SuXon D'uau@ni6nuS 

These executable statements facilitate the writing of recursive 
internal and external functions (see sec. 6, Function Definitions). They 
cause the designation and use of a vector for the temporary storage of 
data and function returns o 

7ol SET LIST T0 Statement 

This statement has the forms 

SET LIST T0 V 

where V is the name of an array element, i=eo, either an array variable 
name or a subscripted variable = The name V specifies the initial element 
in a vector to be used as temporary storage = Consecutive elements will be 
used as required by executions of the statements iiriaich follow in this 
section The number of elements in the list-vector V is deteimined by 
the amount of data the user specifies to be stored there o 

In the following the teraiinology "SAVE statement" will refer to any 
SAVE DATA or SAVE RETURN statement (see secsc 7o2 and 7^3)= The terminology 
"REST0RE statement" will refer to any REST0RE DATA or REST0HS RETURN state- 
ment (see secso 7-^4 and 7»5)» 

A SET LIST T0 statwent must be executed before any SAVE or REST0RE 
statement is exeeutedo 

A SET LIST T0 statement defines the list-vector to be used by all SAVE 
and REST0RE statements until another SET LIST T0 statement is executed^ 
Examples The statement 

SET LIST T0 VECT0R (12) 

results in the assignment of the first element of the current list-vector 
to the 12th element (where is the 0th element), of the vector VECT0Ro 



7oO 
6/20/62 



SAVE DATA STATMENT 

7o2 SAVE DATA Statement 

This statement has the forms 

SAVE DATA ^ 

where #Chas the most general form of an input-output list (see sec 5^ 3s 
The Input-Output List). This statement causes the current values of the 
elements of the list ^ to be stored^ in the order of their appearance from 
left to right in the list, in consecutive elements of the current list- 
vector as specified by the most recently executed SET LIST T0 statement^ 
starting vdth the first available element of the list-vector. I'he element 
which is the first available element of the list-vector is either 

(1) the initial element of the list-vector V as specified by the most 
recent SET LIST T0 statement if no other -SAVE statement has been 
executed since the SET LIST T0 statement ^ 

or 

(2) the element one past the last element used by the most recently 
executed SAVE statement: 

Examples; 

If the following segnent of program is executed; 



SET LIST T0 LIST3 



SAVE DATA MAX, MIN 



SAVE DATA LISTl(O), .,=, LIST1(4) 
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the statement 

SET LIST T0 LIST3 

assigns the first element of the ciirrent list-vector to LIST3(0)= 
The statement 

SAVE DATA MAX^MIN 

stores the values of MAX and MIN as if the following two substitution 
statements were executed; 

LIST3(0) = MAX 
LIST3(1) = MIN 

During the execution of the statements a the next available element in the 

current list-vector is LIST3(2)« 
The statonent 

SAVE DATA LISTl(O), 0=0, LIST1(4) 

stores values of the el^nents of the block as if the following substitution 

statements were executed; 

LIST3(2) = LISTl(O) 

LIST3(3) = LISTl(l) 

LIST3(4) = LIST1(2) 

LIST3(5) = LIST1(3) 

LIST3(6) = LIST1(4) 

During the execution of the statements 3 the next available element in the 
current list-vector is LIST3(7)o 
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70 SA¥E RETURN Statement 

■ This statement has the forms 

SAVE RETURN 

This statement is meaningful ©nly when it appears in a ftinction 
definitions It causes the location in the program which called the 
fxmction to which the function is to return upon ccanpletion to b« stored 
as the next available element in the current list-vectoro The determi- 
nation of the next available element is made as under SAVE DATA Statement 
above = 
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7=4 REST0RE DaTA Statement 
This statement has the form 



REST0RE DATA ^ 
where^^ has the restricted form of an input list (see seco 5.3, The Input- 



/ (■•■•■^ -r^i-A- I -i ^4- 1 



vyuuuui. i-ij-auf^ ij. a. ouucij. ui ii iicuucB iM, j m , oooj H ^n uata exemenx-sj are 
designated by the list in oraer from left to right, and if elements throiagh 
the kth of the list-vector have been used bj SAVE statements, then the 
variable N. is replaced by the value of the k-i+l«st element of the list- 
vector, i=e=, N^ is replaced by the kth element, N by the k-l'st element, 
000 J N^ by the k-n+l'st element Notice that the form of the list implies 
that n names may be designated vd.thout the appearance of n terms in the 
list| a term which is a block designation must be counted to represent as 
many names as elements in the blocks 

After this operation is completed^ these n elements of the list-vector 
are made available in the list-vector. This means that the value of the 
first name in a SAVE statement list executed immediately after a RaST0RE 
DATA statement will be saved in the k-n+l'st element of the list-vectoro 
Example : 

Suppose that the following segnent of program were e^tecuted after j in 
time, the execution of the program segment in the example of section 7o2, 
SAVE DATA Statement! 



REST0FE DATA LIST1(4), -o,, LISTl(O) 



REST0RE DATA MIf.,MAX 



the first REST0RE DATA statement has an effect exactly the reverse of the 

second SAVE DATA statonent in the earlier example and the second REST0RE 
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DATA statement has an effect exactly the reverse of the first SAVE DATA 
statements During the execution of the statements X, the next available 
element in the current list-vector is LIST3(2)c During the execution of 
the statements \i, the next available element in the current list-vector 
is UST3(0). 
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7o5 REST0RE RETURN Statement 
This statement has the form 

REST0RE IffiTURN 

This statement is meaningful only vdien it appears in a function 
definitiono It causes the current last elemeni in the current list-vector 
to be used as the location to which control is returned upon completion of 
the function program, ice., when the next FUNCTI0N RETURN statement is 
executed <> The last element of the list-vector is then made available for 
use by the next SAVE statements 
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12 Advanced Features 

Tnis section is devoted to a discussion of features which are exten- 
sions of the basic features described in preceding sections o An vtnder- 
standing of the preceding sections is presumed hereo 

12»1 Subscription Redefinition 

In order to conserve storage it is sometimes desirable to store less 
than the total number of elements of an arrajo For example, it is often 
desirable to store only half of a synmetric matrix or only the non-zero 
elements of a sparse matrixo 

This can be done using the features of MAD already described if the 
user is content to use linear subscripting exclusively for the matrix in 
question o 

Since the standard array subscripting process used by MAD presumes 
the presence of an entire array stored in the order described in section 
2=2»2s Arrays^ however^ it is not possible to use the array subscripting 
notation for an abbreviated array with the standard array subscripting 
procedureo 

A special subscription routine may be written by the user as a single 
valued internal function, external function or assembly language subrputine 
+ Q ^^^Tni°^ t^e use**^ t^c u°e the arra°''' subscripting notation for an abbrevia"^ **d 
arrayo 
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Hequirements Imp'^sed on th» Subscription Routine 



(1) The arguments for the routine must have the follovdng forms 

(KAME^ S^s, 32? « = = ? S^) 

where NME is the name of the abbreviated array and the S. 's have the values 
of t.he n subscripts resulting from the conventional notation;^ i-rea. the use of 
the subscripted variable 

NAME(S^, S^, =-, S^) 

(2) The routine must be a single valued function whose value is equal 

to the l inear subscript of the array NAME which corresponds to the array 

subscript (S- , S^, oooj S )» 
i -i n 
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A subscript routine written in assembly language may use the fact that 
the location of the dimension vector^ ADIM(0)s appears in the decrement 
part of the calling sequence parameter containing in the address part the 
first l©cation of the array^ A(0) (see sec= 14= 3 > Structure of Subroutines) = 
Examples Suppose that the matrix A ef order n is an upper triangular 
matrixj with all elements below the main diagonal equal to zero» Instead 
of storing the entire matrix^ it is desired to store only the elements above 
and on the main diagonal, i:e= j only the elements a. . for which i < j where 
i is the row index and j is the column indexc The elements are to be stored 
by row^ i=e« so that 3 

(1) if the element a. .. where j < n, is in matrix location A(k) then 
the element a. is in matrix location A(k+l)| and 

(2) if the element a. is in matrix location A(^j then the element in 

%n -^ 

matrix location A(i'^l) is the element a ., . , » 

*^ 1+1; 1+1 

Theftj assume further that the matrix A is stored in the MAD array A 
with the base element A(lsl) equivalent to the element A(B), and the array 
element A(B-=1) is set equal to zero whenever ccmputation is performed using 
the matrix A. 

If SUBFo is the name of a function written to calculate the correct 
linear subscript given the conventional matrix subscripts i and j then the 
desired value of the function is given by 

(B-l<> if i > j 

SUBF» (A,i,j) -l^(^^,^ ,w, ,s 

H2(n+lj'=-ij(i-l) ^ . o ^r, . f> . ^ , 

^ ^' ^2 +J-i+Bs if 1 < J 

where i and J may be^ as usual, any integer expressions^ but whose values 
should satisfy l<ij(j<no 

The value of SUBF= (A»ijj)j then^ is the linear subscript of the 
element A(ijj) if A(i,j) is one of the elements which is actually stored^ 
and is the linear subscript (B-1) of an element whose value is equal to 
zero J which is all that is necessary^ otherwise » 

The 'subscripting function must have access to the values of n and Br 
If it is an internal function they can be gotten from the dimension vestoj- 
of A by using the name for it, say ABIMj y^idh the main program used. 
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OtherKisej for external and assembly language functions j the values of n 
and B can be made available through PRgiGRAM C(MM0N. 

12rl-. 2 Use of the Subscription Routine 

Once a routine is written to perform subscripting in other than the 
conventional manner, it is still necessary to indicate to MAD that this 
new roui^ine is to be used instead of the standard one: 

Aa always the array,. A,, in question must appear in a DIMENSION 
declaration in the following form 

DIMEKSIldN -" A(PjADIM), ... 

where ADIM is the name of the dimension .-ector of A = 

In this case ( that is when redefining subscription) a' statement of 
one of the following two foims must be used to se^ at least the first 
slen„er.t cf the diicension vectors 

VECT0R VALUE.S ADIM = SUBF« , d 

\'ECT0R VAIUES ADIM = SUBF. , d^, d^, o ,. ,-, , d^ 

where SUBFc is the nam® of the subscription function to be used with the 
array A and the d^s are integer constants-. 

As discussed in section 3 = 6c 2c, Matrix DIMENSION Declaration, it, was 
possible to compute or read in as data the entries in the dimension vector^ 
This still is possible when redefining subscription but each of the entries 
aeseribed m section 3~6o2 must be stored in the dimension vector in the 
element whose linear subseriprc is eqiial to 1 plus the Ziinear subscript, of 
the element a.s described in section 3-6o2o 

The name of the subscription! function must appear in the first element 
cf the dimension vector and may be put there in no other way than by 
presetting with ¥ECT0R VALUES. 

The first form of VECT0R VALUES declaration listed above may be used 

when the dimension entries are to be computed or read in| even if the array 

dimensionj in the element of d, is tc be computed or read in^ an integer d^ 

possibly a dianmy.; must appear. The second form is that used when elements 

can be preset, where d-, is the array dimension, d is the linear subscript 

of the base element,, etc^ 
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12=1.3 Available Subscription Routines 
The following two special siabscription routines are available in the 
subroutine library? 

(a) SIMMc 

The subroutine SIMM^ calculates the correct linear subscript for a two 
dimensional s symmetric matrix for whish only the elements above and on the 
main diagonal are stored, by rowso 

SIMMo has access to the dimension vector of the matrix in question so 
no information other than the arguments is necessary. 

(b) TRANSPc. 

The subroutine TRANSP: calculates the correct linear subscript for a 
complete two dimensional matrix^ A^ stored as the transpose of A| ioe. the 
transpose of A is stored by rows or A is stored by columns= 

For exaiiple^ a reference to the element A(6j2) using subroutine TRANSP. 
will result in the linear subscript for the element A(2^6)c 

TRANSP -. has a-ccess to the dimension vector of the matrix in question 
so no information other than the arguments is necessary: 
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13 Recommended Programming Practices 

This is an expandable section designed to have added to it, from time 

to time, recommendations for programming procedures which -Mill result in 

increased efficiency in some manner, eog= speed of canpilation, speed of 

execution, conservation of storage, etc= 

For this reason this section contains a separate table of contents 
By definition, this section is incomplete at any given time and there 

may exist methods .;','-. . :r!.'-n these which appear here which could be used 

to accomplish the same ends. 
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13.1 Table of Contents, Section 13 
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13.2 Efficient Programming of Boolean Expressions 

13 '3 Use of Parameters in Function Definitions 
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13-»2 Efficient Programming of Boolean Exprsssions 

The object program produced by the MAD translator evaluates the terms 
of a Boolean expression frcan right to left a 

It is clear thatj in a Boolean expression of the fonn 

•Tl ojSfeo T2 »0ao — =0R. Til 

which may be called an "o0R„ expression" ^ where the Ti«s are any pennissible 
terms and where the Boolean operations 

=EX0Ro 
= AND. 

do not appear in the expression, the value of the .0Ro expression is 
certain to be IBj true, as soon as, in the evaluation process, the value 
of some Tl has been found to be IB; The object program produced by the 
MAD translator evaluates the terms of an c0R= expression until the value 
of some term is IB or until all the terms have been evaluated, whichever 
happens firstc Thus possible remaining terms j, after the evaluation of some 
term whose value is IB, are not evaluated needlessly-. 

Thiss together with the fact that the tenns are evaluated from right 
to left as mentioned above ;, implies that^ if the programmer has any way of 
judging, he should write the terms of an o0R, expression from left to right 
in order of ascending likelyhood of truthfulness, ThuSj in an expression 

Tl .0Ro T2 

if T2 is likely to be "true" (have value IB) more often than Tl, the 
expression should be written as it appears above^ otherwise the order 
should be reversed;, i^e.- 

T2 o0Rc Tl 

It is clear that^ in a Boolean expression of the form 

Tl oANDo T2 cANDo — »ANDo Tn 

which may be called an "<.AND» expression", where the Ti's are any permissi- 
ble terms and where the Boolean operations 

=1X0R= 
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do not appear, the value of the oAND» expression is certain to be OB, 
false, as soon as^ in the evaluation processj the value of some Ti has been 
found to be OB^ The object program produced by the M4D translator evaluates 
the terms of an cAND=. expression until the value of some term is OB or 
until all the terms have been evaluated j whichever happens first « Thus 
possible remaining terms, after the evaluation of some term vdiose value is 
OB, are not evaluate,d needlessly. 

This, together with the fact that the terms are evaluated from right 
to left as mentioned above, implies thatj, if the programmer has any way of 
judging j, he should write the terms of an »AND= expression from left to 
right in order of descending likelyhood of truthfiilnesso Thus, in the 
expression 

Tl »AND» T2 

if T2 is likely to be "false" (have value OB) more often than Tl, the 
expression should be written as it appears above, otherwise the order 
should be reversed, i«eo 

T2 .AND. Tl 
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13 3 Use of Parameters in Function Definitions 

During compilation of a function definition, MAD records in a table, 
called the Parameter Use Table, all references in the function definition 
to the dianmy arguments of the functions which will be references to the 
actual argTjnents when the function is used = 

These entries are used to compile into the function subroutine, 
instructions to initialize the references upon each call of the functiono 

If too many references are made the Parameter Use Table may not be 
able to contain them all and compilation cannot be successful (in which 
case the coinment PARAMETER USE TABLE EXCEEDED is printed) o 

Thus it may be necessary to minimize the number of entries in order 
to make compilation possible and it is desirable to minimize the number 
of entries both to shorten the object program and speed up the execution 
of the subroutine^ This may be done in any of the following wayss 

(1) If the dummy variable X supplies data to the subroutine^ i=eo is 
input to the subroutine^ and is referred to several times, one may use 
the substitution statement 

Y = X 

where Y is not a dummy variable j immediately after entry to the subroutine 
and use Y instead of X thereafter^ Then initialization of references to 
X is done only in the substitution statement and there is only one entry 
for X in the Parameter Use Table =. 

(2) If the address of a variable is needed, as in the case of an output 
argwnent or an argnnent which Is an array name^ it is not possible to use 
the method of (l) above = Instead it may be possible to put the variable 
or array in PRQ^GRAM C01M0N5 by means of identical declarations in both the 
main program and subroutine if the subroutine happens to be either an 
external function or an assembly language subroutineo Then the variable 
or array will not be an argument at alio 
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14 Mechanics of Using MAD 

14ol Card Format 

In order to compile a source language program using the MAD translator 
the program must be punched on cards » These cards have the following 

format : 



Columns (inclusive] 


Content 


1-10 


Statement Label 


11 


Remark or Continuation Designation 


12 - 72 


Statement 


73-80 


Identification 



14»1»1 Statement Label Field 
Statanent labels may be punched anywhere in columns 1 - lOo Spaces 

are ignored* 

14c lo 2 Column 11 

Except for a "Hemark"', a statement must begin on a card which has a 
blank in column llo 

14»1'2=1 Remarks Cards 
If an "R" is punched in column 11, the card is a Remark Declaration 
(see seco 3=1) » 

140 1" 2c 2 Continuation Cards 

A card which has a decimal digit Oj 1$ — , 9, punched in column 11 
is a continuation cardo Continuation cards may be used in the event that 
a statement cannot be put on one cardo The order of the digits appearing 
on consecutive continuation cards is not significant,* the statement is 
ordered by the physical order of the cards comprising it. 

A maximum of nine continuation cards may be used, so that a statement 
may be ccmpri' ed of no more than ten cards« 
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14=1=3 The Statement Field 

MAD language statements may be punched anywhere in colunms 12 - 72 and 
only in columns 12 - 72; Spaces ar® ignored except when they appear bstween 
pairs of dollar signs (*"$"' s).. 

14»lo4 The Identification Field 

The information in the identification field is not translated by MAD= 
It is transcribed onto the source language listing produced during compil- 
atlon» The user may punch any legitimate characters in this fields 

It is good practice to include in this field identification information 
and sequence ntmbers whieh define the order of the cards in the prograjn. 
deck» 
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1a. c, 2 Diagnostics 

During the prosess of translation many kinds of errors in the formation 
of statements and the allocation cf storage can be detected-. To understand 
this error detection and the subsequent printing of diagnostic soraments seme 
knowledge of the structure of the translator is helpful ^ The translation 
from statements to machine code is accomplished in three major sections? 

(1) The deccxaposition of the original statements into arrays of binary 
operations and pseudo-operations s 

(2) The analysis of all of the declarative information in order to 
allocate variable storage and identify the arithmetic types (i<,e=. 
modes) of variables. 

(3) The combination of the information produced from (1) and (2) to 
translate the arrays to relocatable binary programs » 

When an error is encountered in one of these sections the translation does 
not proceed to the next sectiono However 5 insofar as possible, the entire 
spt cf ^rgtements is processed through the section in which the error is 
detected and therefore more than one error may be detectedo It should be 
understood thenj that not all detectable errors may be found because; 

(a) They are detectable only in a later stage of the translationo 

(b) Scaae types of errors make it impossible to attemp* ■^'urther 
detection within the section in which it occurs I 

(c) One error may actually obscure another errors 
Occasionallys an error in one statement may be such that it causes the 
translator to misinterpret a second statement; thus giving an error indica= 
lion even though no error exist? i.n tne later statements 

The printed diagnostic comment may very often have an alternative or 
ambiguous fonrio This results from the fact that it is frequently not 
possible to determine what form was intended, ■nerelj' that the present 
structure is not admissible? and therefore some of the alternative possibil= 
ities are suggested by the ccmmentc 

After MAD has ccanpieted translation of a source language programs a 

list is printed of all variable names which appeared in the program, but 

which appeared only onceo This list does not include names appearing in 

any of the following declarations: 
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PR0GRAM C0K0n 

ERASABLE 
DIMENSI0N 
VECT0R VALUES 
EQUIVALENCE 
Variable names which appear in the list are all assigned to the same 
location, under the assumption that they are not purposely used for anything 
except perhaps redundant labelling of statements.- 

The list is a valuable debugging aid due to the fact that it is very 
likely that in it misspelled names will appear = 
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i4-3 Structure of Subroutines 

The information in this and the following sections is to be found in 
much greater detail in other installation write-ups^ However^ the following 
sections should be sufficient for the general use of MAD= 

Subix>utines which are written for use by MAD programs, vdiether written 
in MAD as functions or in assembly language code, must- be relocatable and 
must operate from the calling sequences the translator produces ^ Consider, 
for example the function cal.'L 

FN. (A,B,C) 
which might appear in the body of a statement o Assume that B is an array 
which has an associated dimension vector BDIM= Using assembly language 
notauXon xor ij^-tustrative purposes^ '^^he caxxxng sequence produced wcuj^d bes 

TSX FN, 4 

TXH A 

TXH BjCjBDIM 

TXH C 
Input-output routines utilize two types of parameters, the regional 
and single variable types : In addition an error return is given as well as 
a format specification location. The paraneter operation code used is STR 
and the end of the parameter list is indicated by an STR operation with a 
blank address o Thus the statement 

READ FORMAT MT, BETA, X(l) = = .X(100),, K 

vould produce the calling sequence 

TSX READ, 4 

STR ERR0R 

STR FMT 

STR BETA 

STR X-ljOjX-lOO 

STR K 

STR 

On occasion it is useful to use the regional notation in subroutines 
which are not in the input-output category^ for examples G. (GAMMA, DELTA^ 
Z(10). =.Z(20))c The calling sequence would be 
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TSX G,4 

TXH GAMMA 

TXH DELTA 

TIX Z-10,0sZ-20 

It is important to notice that in this example, as v/ell as in the first, 
the parameters s if executed as instructions, would produce no operation. 

It is beyond the scope of this manual to discuss the structure of 
relocatable programs- It is sufficient to say that a relocatable program 
must containj in addition to the actual instructions in the program, 
infoimation as to which addresses must be relocated at the time of loading 
for execution and which addresses must not. In addition, the first card 
(or record) of such programs must contain information about the size of 
the program, the number of subroutines it calls on, the amount of storage 
it will share with other subroutines, the location of the list of sub= 
routines it calls on, and the names by which the routine itself is referred 
to. The symbolic names of the subroutines called on must appear as the 
first words after this infonnation= 

The execution of MAD programs requires the use of a loading routine to 
relocate and store the program and subroutines. A slightly modified BSS 
FORTRAN loader is automatically produced by the operating system in which 
MAD is imbedded ^ Also there are certain subroutines which may be auto- 
matically called for by a MAD program without an explicit reference to them 
in the source program. 
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SYSTMS SUBROUTINES 

14« 4 Systems Subroutines 

The use of the following nam^ f or functions (subroutines) should be 
avoided except where the operation is the one indicated here^ 

SYSTEM - Entry to this routine causes a return to the operating systouc 
The END 0F PR0GRAM statement produces a call for this routine. 
ERE0E - Entry to this routine also causes a return to the operating 
systenic Howeverj if a dump of storage was requested of the 
operating system such a print of storage will be produced before 
the return to the systanc The ERR0R RETURN statement may 
produce a call for this routine < 
TAPEWE - Entry to this subroutine causes BCD information to be written 
on tape= The argiments are; (1) location of the format 
specification, (2) tape number, and (3) a list of variables to 
be written. The WRITE BCD TAPE statement produces a call for 
this subroutine^ 
PRINT - Similar to TAPEWR above except that the use of the peripheral 
output tape is implied- The PRINT F0RMAT statement produces 
a call for this subroutine^ 
TAPERD - Entry to this subroutine causes BCD information to be read 
from tapeo The arguments are? (1) location of the format 
specification J (2) tape number, and (3) a list of variables 
to be readc The READ BCD TAPE statement produces a call for 
this routine: 
READ - Similar to TAPERD above except that the use of the peripheral 
input tape is implied The READ F0RMAT statement produces a 
call for this subroutine. 
PUNCH - Entry to this subroutine causes BCD information to be written 
on the peripheral output tape. The arguments are the same as 
those of PRINT and the PUNCH F0MAT statement produces the call 
C01MNT- Entry to this subroutine causes BCD information to be printed 
on the attached (on-line) printer c The arguments are the same 
as those of PRINT except that the line spacing is not given in 
the format specification since an automatic 1/6 page skip is 
produced. The statement producing the call is PRINT 0N LINEo 

14o6 
6720/62 



SYSTEM SUBROUTINES 

SETE0F - This name may appear explicitly in the form S£TE0F.(S)s mere S 
is an expression of statement label mode designating the point 
of rettirn when an end-of-file is encountered during the reading 
or writing of magnetic tape, 

SETETT = This name may appear explicitly in the form SSrETTc(S), vihere S 
is an expression in statement label mode designating the point 
of return when an end-of-tape is encountered during the reading 
or writing of magnetic tapso 

SETERR = This nans may appear explicitly in the fonn SETERR=(S)s, x*here S 
is an expression of statement label mode designating the point 
of return after an illegal data character or illegal fonnat 
specification has been encountered by an input-output sub- 
routine = 

It should be understood that the subroutine described above may be 
called by other subroutines as well as statements^ Thusj for example^ 
SETE0F is called by TAPERD^ TAPEWR calls PRINT, and all of the input- 
output routines call ERR0Rc 
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15 Exan^les of MAD Programs 

The following examples illustrate how some programs may 'be written 
in the MAD langiage. Since they were written to illustrate as many features 
of the language as possible^ they are not necessarily the most efficient 
or elegant programs which could have been written. 
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SCIENTIFIC 

15-1 Sdientific Exam|3les 

Example 1 

2 
Problem s To solve the quadratic equation ax + bx + c = for various sets 

of coefficients a^ h, and c. 

Analysis ; Let x^ and x be the two roots of the equation. Then their values 

are found by the formulas. 



-b +■ Vb^ - ij-ac - b - Vb^ - it-ac 



X-, = ' — X, 



1 ~ 2a 2 ~ 2a 

whenever a / 0. The single root X-, of the equation when a = is x-, = -c/b . 
The input values of a, b, and c are printed immediately after they are brought 
in to help in finding trouble spots during the development of the program 
(not as necessary here as in longer problems^ but a good idea ' ) • 
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Print 



■~C__l£. 



^ 






a 



Print 
"Linear 
Equation" 






< 




/7^ 
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Print "FEhI 
SOLUTIONS" 
("b +1fd)/2a 

"-b «1fd)/2a 



h^ 




Print "COMPLY 

SOLITIOKS" 

R(X^)=-b/2a 

R(X^)= -b/2a 
l(X^)-T^/2a 





T 



?'Iote- 



R(X^) ana 



(X, ) are zhe 



real and imaginary par^s 
of X-, , and J sitnlli,rlj>^ 
RCXg) and iCx^,) are the 
real and imaginary parts 
of Xg. 



SCIEHTIFIC 



R 

RMAIN PROGRAM 

R 

GAMMA READ FORMAT INPUT»A,B,C 

PRINT FORMAT GHECK,A,B9C 

WHENEVER A .NE. O9TRANSFER TO ALPHA2 
ALPHAl PRINT FORMAT LINEAR, -C/8 

TRANSFER TO GAMMA 
ALPHA2 D = B .P. 2 -4,*A*C 

WHENEVER D *L. 0.»TRANSFER TO BETA2 
BETAl PRINT FORMAT REAL, (-B+SQRT. ( D ) ) / ( 2.*A) , ( -B-SQRT . (D) ) /{ 2 .*A ) 

TRANSFER TO GAMMA 
BETA2 PRINT FORMAT COMPLX»-B/ ( 2.*A) ,SQRT. (-D I / ( 2 .*A ) , 
1-B/(2.*A) ,-SQRT»{-D)/{2.*A> 

TRANSFER TO GAMMA 
R 

RFORMAT SPECIFICATIONS 
R 

VECTOR VALUES INPUT = $ 3F10.4*$ 

VECTOR VALUES CHECK = $ 4H0A = F10.4,S8, 
13HB = F10.4,S8,3HC = F10.4*$ 

VECTOR VALUES LINEAR = $21H0LINEAR EQUATION, X = F10.4*S 

VECTOR VALUES REAL = $21H0REAL S0LUTI0NS,X1 = 
1F10.4,S8»4HX2 = F10.4*$ 

VECTOR VALUES COMPLX = S19H0COMPLEX SOLUTIONS, 
1S4,7HR{X1) = F10.4,S8»7HI (Xl) = F10*4,S8, 
27HR{X2) = F10.4,S8»7HI(X2) = F10.4*$ 

END OF PROGRAM 



DATA 






4. 


-8. 


4. 





5. 


10. 


1. 


1. 


1. 
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Exaaple 2 
Problem? A logieaZ. (Boolean) expression such 

T = (P .AND. Q) .6-R. (.I0T. P ,ANB. E .AIS. S) .0R. (R ,0R. P) 

will have a value TRUE or FALSE (represented here by IB and OB^ respectively) 
depending on the "input values" of the variables- involved; PjQ,R,S. Thus^ 
if P = IB^ Q = R = S = OB, then the total expression T will have the value 
IB. The entire table of outputs for all possible inputs would be as follows; 



P 


Q 


R 


r"^ — 
s 


T 


OB 


OB 


OB 


OB 


03 


OB 


OB 


OB 


IB 


03 


OB 


03 


IE 


OB 


IB 


OB 


OB 


IB 


13 


13 


OB 


IE 


OB 


OB 


03 


OB 


IB 


OB 


IB i i OB 


OB 


IB 


IB 


03 1 13 


OB 


IB 


IE 


IB 1 

i 


i ^^ 


IB 


OB . 


OB 


OB 


IB 


IB^ 


03 


OB 


IB 


13 


IB 


OB 


IB 


OB 


13 


IB 


OB 


IB 


IB 


IB 


IB 


IB 


OB 


OB 

1 


IB 


IB 


IB 


OB 


IB 1 


IB 


IB 


IB 


'■ U 


OB ! 


IB i 


IB 


IB 


IB 


IB ' 1 IB 
M 



The problem is to write a program to generate the entire "truth table" 
for the given expression T. 
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THROUGH A, 
FOR VALUES OF 

P = OB, IB 



THROUGH A, 
'FOR VALUES OF 

= OB, IB 





a. 



THROUGH A, \ 
FOR VALUES OF 

R = OB, IB / 



THROUGH A, 
FOR VALUES OF 

S = OB, IB 




(J>- 



PRINT 
P,Q,R,S 



PRINT (P.AHD.i^ 
.OR. (.NOT. P 
.AND. R .AM). 
S) .OR. (R .or 
P) 
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PRINT FORMAT HEADER 
BOOLEAN P.O»RfS 

THROUGH AsFOR VALUES OF P = OBflS 
THROUGH AsFOR VALUES OF Q = 0B»1B 
THROUGH AsFOR VALUES OF R = 0B»1B 
THROUGH A»FOR VALUES OF S = OB. IB 

PRINT FORMAT TABLE »P »Q » R »S ,( P .AND. Q) .OR. (.NOT. P .AND. R 
l.AND.S).OR. {R .OR.P) 

\iCf-Tr\n \/Aiiicc utTAnPD — ciui.Cin-lua.Cin.iUrs-ClA.TUD.cin.iuc 

VLV^iwrV vr*L_wuo !!L_rAU'i— 1\ — »*^J,ttX^vJiW>J.tir TI^JXV'ill^V7vJXV^»XlH\»vJiW5XllsJ 

1S15»1HT*$ 
VECTOR VALUES TABLE = $1H0 »4 ( S 10 . 1 1 ) »S15 » 1 1*$ 

END OF PROGRAM 



Ip .6 

f. /on I'i^o 
-^t ■ — '/ -^ 



SCIENTIFIC 

Note: Although it would have meant only a slight change in the format 
information^ no attempt was uiade here to label the "0" and "1" that 
print as values in the table as Boolean, i.e., "OB" and "IB". This 
points up the fact that internally OB and IB are stored as and 1, 
respectively. Also, the statement 

NS^RMAL M0DE IS BjZ50:LEA3I 

could have been used as the second statement of this program instead 
of the B^^LEA'N declaration . 
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Example j 

Problem; To approximate y f (xj "by Simpson's Rule, for an arliitrary 

' — = a 

interval [a, "b] using 51 equal subintervals (where W is an arbitrary even 

integer ) . , 

Analysis '; By Simpson's Rule, J f(x)dx '^ -^(Yq + %i "^ ^^2 "*" ^"^5 + ••• 

+ S^.l + %^> 

where y^ = f(x^), and a = x^, x.^^, ..., x^ =. b are the partition points of 
the interval [a, bj . 

Method i We shall write the program in the form of an external fionction, so 
that it could be used with any other program. The evaluation of f (x) may be 
accomplished by another external function or an internal function. 



15.8 
6/20/62 



Flow Diagram: ■ 



vo 
o^ o 




o 

n 

M 

t-H 
O 
CO 




T 



'THROUGH ALPHA' 
FOR x=a+h,2hj 

x>b 



S^=S-^+f (x) 



Sg=S2+f(x+h) 




FUNCTION RETURN 
|(f(a)+i^S^+2S2 
^f(b)) 





EXTERNAL FUNCTION 




INTEGER N 




ENTRY TO SIMPS, 




H = {B-A)/N 




SI = 0. 




S2 = 0. 




THROUGH ALPHA»FOR 




SI = SI + F.{X) 


ALPHA 


S2 = S2 + F. (X+H) 



SCIEIWIFIC 



{A,B,N,F.) 



X = A+H, 2.*H» X .G. B 



FUNCTION RETURN H* { F , { A )+4.*Sl+2 .*S2-F. ( B ) ) /3. 
END OF FUNCTION 

If, for some reason, the integral of sin Jx - cos (rx + l) were needed if < r < 3, 
and the integral of sin 5x - cos rx otherwise, the program might then "be as 
follows : 



READ 



INPUTS. A, B,N,R 



R .LE. 3. 
!,R,SIMPS.(A,B,N,F1 



READ FORMAT 

INTEGER N 

WHENEVER 0. .LE. R .AND. 

PRINT FORMAT RESULT, A, B»l 

OTHERWISE 

PRINT FORMAT RESULT , A , B » N ,R ,SI MPS . { A , B ,N , F2 . 

END OF CONDITIONAL 

TRANSFER TO READ 
R 

RDEFINITION OF FUNCTIONS 
R 

INTERNAL FUNCTION F1.{X) 

INTERNAL FUNCTION F2.(X) 
R 

RFORMAT SPECIFICATIONS 
R 

VECTOR VALUES INPUTS = $2F12 .4 . 1 6 ,F12 .4*$ 

VECTOR VALUES RESULT = $23H1 FOR THE INTERVAL FROM 
1F12.4,3H TO F12.4.5H WITH I6,38H EQUAL SUB-INTERVALS AND 
2PARAMETER F12 .4/29H0THE VALUE OF THE INTEGRAL IS F12.4*$ 

END OF PROGRAM 



= SIN. {3.*X)-C0S. 
= SIN. {3.*X)-C0S. 



{R*X+1.) 
(R*X) 



EXTERNAL FUNCTION (A,B,N,F.) 
INTEGER N 
ENTRY TO SIMPS, 
H = (B-A)/N 

51 = 0. 

52 = 0. 

THROUGH ALPHA, FOR X = A+H,2.*H, 



X .G. B 
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SI = S1+F.(X) 
ALPHA S2 * S2+F.{X+H} 

FUNCTION RETURN H» ( F. ( A)+4.*S1+2.»S2^F. ( B) ) /3. 
END OF FUNCTION 

$ DATA 

2. 10 10. 

An alternate way to write the first eight lines of this program, illustrating 
one use of the FUNCTION NAME mode, would be: 

READ READ FORMAT INPUTS .A, B,N,R 

INTEGER N 
WHENEVER 0.«LE. R .AND. R .LE. 3. 

S = Fl. 

OTHERWISE 

S = F2. 

END OF CONDITIONAL 

KRini rur^mAi r\coui_ i »m»o tm trx ♦ on'^rij* ««fa»ri»oj 

TRANSFER TO READ 

FUNCTION NAME S 
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Example k 

Problem; To find a real solution (if it exists) of the equation f (x) =0 
(where f is a continuous fimction) on an arbitrary interval [a^ b], pro- 
vided the roots (if there are more than one) are at least e apart. 
Analysis ; we specify a, o, and e as parameters. The method used will be 
"ha If -interval convergence," in which the function is evaluated at x = a, 
and then the interval is scanned for a change of sign* in the value of 
f(x). If no change of sign is found, the scanning is repeated with a step 
size for searching equal to one-half the previous step size. If the step 
size becomes smaller than €, and no change of sign is found, the process 
is terminated, and comment is printed; N0 S0LUTI0N. 

If a change of sign is found between x and x^, the value of f is 

computed at x„ - ^ ^ , i.e., the midpoint of the interval of uncertainty 

2 

[x, , Xp]. We then determine which of the intervals [xy, x„], [x^, Xj,] now 

contains a change in sign. We then compute the value of f at midpoint 
of that smaller interval, etc., until the interval being considered finally 
has length less than €, at which time either end may be taken as the 
solution with an error less than €. 

The method used here to handle the x.^ computation is perhaps not 
the most obvious one- It consists of a simple loop in which the value 
X is adjusted by h' = ^, then h" = -^ = j-, etc., until h is small enough. 
The adjustment of x is either to the left or right, depending on the 
occurrence or non-occurrence, respectively, of a change of sign between 
f(a) and f(x). 

It should be understood that this method may not find a root which 
is one of a pair of roots which either coincide or are less than € apart. 



*A change of sign is detected when the n\mbers involved have a negative 
product . 
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Print 
a, b, e 










\ 


= f(a) 











'THEOUGH ALPHA) 
FOR 
b-a 



i = - 



2 
h < e 



V [ a 1 



I THROW 
r^\_J ALPHA, 



THROUGH 

V^y \ x=a+h,h,x>b 

\ I 





THROUGH 
SIGMA, FOR 
, h -h 



2' 2^ 

h < e 



f(x):0 



® 




"WO 
SOLUTION" 




Definition 



SIGN.(Z) = Z/lz| 
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(it is assumed here that f (referred to as F. in the program) will he defined 
as an internal f-unction.). 



INTERNAL FUNCTION F.(Z)= Z .P. 2 - 2. 
PSI READ FORMAT ABEPS, A,B,EPS 

PRINT FORMAT INVAL» A»B»tPS 

YA = F.(A) 

THROUGH ALPHAiFOR H = ( B-A ) /2 . »-H/2 . .H .L. EPS 

THROUGH ALPHA, FOR X = A+H,H, X .G. B 

WHENEVER F.(X) .E. 0., TRANSFER TO ETA 
ALPHA WHENEVER YA*F.(X) .L. 0., TRANSFER TO DELTA 

PRINT FORMAT NO ROOT 

TRANSFER TO PSI 

R 

RTHE NEXT SECTION IS ENTERED WHEN A CHANGE 

ROF SIGN IS FOUND 

R 
DELTA THROUGH SIGMA, FOR H=H/2 . ,-H/2, tH .L.EPS 
SIGMA X = X+SIGN. (YA*F.{X) )*H 
ETA PRINT FORMAT ROOT»X 
TRANSFER TO PSI 

R 

RDEFINITION OF SIGN. FUNCTION 

R 

INTERNAL FUNCTION SIGN.tZ) = Z/.ABS.Z 
R 

RFORMAT SPECIFICATIONS 
R 

VECTOR VALUES ABEPS = $ 3F12.^*$ 

VECTOR VALUES INVAL = $18H1 INPUT*VALUES, A = F12.4,S3, 
13HB = F12.4,S3,5HEPS = F12.4*$ 

VECTOR VALUES NO ROOT = S12H0NO*SOLUT I ON *$ 

VECTOR VALUES ROOT = $14H0S0LUTI0N, X = F12.4*$ 

END OF PROGRAM 

$ DATA 

1. 2. .01 



15.1^ 



seisNTiFie 



Example $ 



Problem : Find the transpose A' of an nxn matrix A = (^i^)- 

Amlysis ; If we write A' = (t>j j}i then ^i* = &<< • W® shall interchange 

symmetrically placed pairs of elements, leaving untouched elements on 

the main diagonal. The program will be in the foito of an external function. 
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o H 

ro 



Flow Diagram: 




THROUGH 
BETA, FOR 

K=1,1,K>N 



THROUGH 
BETA, FOR 

I=K+1,1;. 

I > N 





a 



a(i,k)=a(k,i) 



A(K,I) = Z 




/FOTdTior 

RETURN// 



O 

HI 

n 

M 
•^ 
M 
O 



SCIENTIFIC 



EXTERNAL FUNCTION (A»N) 
ENTRY TO TRANS. 

THROUGH BETA, FOR K = !♦!♦ K .GE. N 
THROUGH BETA, FOR I = K+1,1, I .G. N 
Z = A( I ,K) 
A(I»K) = A(K,I) 
BETA A(K»I ) = Z 

FUNCTION RETURN 
INTEGER N,K,I 
END OF FUNCTION 

Wo dimension information is needed for A, since it is an argument in a 
function definition program. This function would be called in a stateme'nt 
of the form EXECUTE TRANS. (A,N). 
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Example b 

Problem: Multiply the matrix A = (a. .) by the matrix B = (b. .) to produce 

"""^ 
the matrix C = {c .), i.e., C = A*B. Assume that A has dimensions m x n 

with m^n < 1500, B has dimensions n x p, with n-p < 1500, and C has 

dimensions m x p, with m-p < 1500. 

Analysis : An element c. of C is computed by the formula 



c. . = 2j T a., b, . 
ij Z=l xk kj 
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Flow Diagram; 



o 

M 



H 
O 
K3 



( Start y 



7!^ 



Read ra,n,p. 



II**" 



3 

np 



Print 
all in- 
put v; " 
-values 



'THROUGH Q, 

FOR 1=1,1, 

1 > m 



cy 
.^ o 

- CM 



^THROUGH Q, 
FOR j=l,l, 

J >P 



^ij=° 




THROUGH Q, 
FOR k = 1,1 

k > n 



SCIENTIFIC 

DIMENSION A( 1500 »ADIM ) ,B { 1500* BDIM) 9C(1500,COIM) 
EQUIVALENCE (N9ADIM{ 2) ) » ( P »BDIM { 2 ) j 
INTEGER I»J»KjM,N»P 
VECTOR VALUES ADIM = 290.0 
VECTOR VALUES BDIM = 2»0,0 
VECTOR VALUES COIM = 2»0*0 
READ READ FORMAT INPUT»M,N »P »A ( 1 » 1 ) o . .A { M,N ) 98 ( I » 1) . . . B( N ♦? ) 
PRINT FORMAT INVAL 9M9N9P » AC 1» 1 ) . ..A (M.N ) ♦Bfl » 1 ) . . .B { N9P) 

r* n t M / o \ — n 

THROUGH Q9 FOR I = 1»1» I .G. M 

THROUGH Q» FOR J = I9I9 J .G« P 

CCIfJ) = 0. 

THROUGH Q» FOR K = I9I9K .G. N 
Q CdsJ) = C(i»J5 + Af I»K)«-B{K9J) 

PRINT FORMAT RESULT9C{ I » 1 ) .. .C ( M9P) 

TRANSFER TO READ 
R 

PFORMAT SPECIFICATIONS 
R 

VECTOR VALUES INPUT = $3 I 4/ ( 6F 12 .4) *$ 

VECTOR VALUES INVAL = $13H0INPUT*VALUES/4H0M = I6»S6»3HN = 16 
1»S693HP = I6//(1H098F13.4)*$ 

VECTOR VALUES RESULT = $9H1C MATRIX// ( 1H0»8F13 .4 )*$ 

END OF PROGRAM 

$ DATA 
2 3 3 

1. 2. 3. 4» 5» 5» 

7. 8. 9. 10. 11. 12. 

13. 14. 15. 



15.20 
6/20/62 



SCIENTIFIC 

Example J 
Problem ; Solve a system of n < 20 simultaneous linear equations in n 
unknowns J assuming that one does not encounter a zero on the main diagonal 
of the coefficient matrix during the solution process. 

Analysis ; We shall use a Jordan EliKdnatlon Method., in which each diajronal 
coefficient is used to "clear" all other coefficients. in its column to 
zero hy appropriate multiplications and subtractions . Since we shall 
divide the "clearing row" by the diagonal element in that row before clear- 
ing the column, we shall finish the process with only a diagonal of ones 
and the solution to the problem as the resulting right hand side of the 
equations . 

We denote the system of equations to be solved by: 

a^TX-, + a^^x^ + . . . + a^ X = a-. , 

11 1 12 2 In n l^n+l 

(l) Sm-Xn + ^^o^o + . ■ . + a^ X = a. , 

^ 21 1 22 2 2n n 2,n+l 



a ^x^ + a ^X- + . . . + a x = a , 
nl 1 n2 2 nn n n,n+l 

we divide the first row by its diagonal element a, , . Then to clear 
api to zero we subtract a^ times the first row from the second row, and 
so on. In general, to clear a., to zero (after row k has been divided 
by ai^v)j ^^ subtract a., times row k from row i (i/k). A typical element 
a. . is thus transformed each time by the formulas: 

^2) \j=^kj/\k 

^J = "ij-"ik^kj ^^^^) 

where the value of a . in (5) is the result of (2). These transformations 

kj 

are performed for k = l,2,...,n. For each (fixed) k, we will let i = 1, 
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2, •..^•'l^Is+l, ...,n,.-so as to operate on all rows exeept i=k. While 
transforming each row we will cycle on j from right to leftj i.e., j = 
n+l,n,n-l, .. -A and we stop at g=K since for j < k there is no change in 
the matrix. 

The array 



V^12 ••• %n+l 



A s (a,..) = 

• ij 



V^n2 



m 

is called the "matrix of coefficients" of the system (l). 

It should be understood that this method, involving the assumption 
of no zeros on the diagonal and not searching for the largest element of 
a row to use as a divisor (to minimize round-off error), is not satis- 
factory from a mathematical point of view. It could serve as a basis 
for a larger, more complete program, however, and serves here only as 
an example problem. 
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Irtint 
input 
values 



THROUGH B, 
FOR 

k = 1,1. 
k > n 



THROUGH C, 
FOR j=n+l. 





THROUGH B, 
FOR i = 1^ 
1, i > n 



h 



{UEzy^ 






THROUGH D 
FOR j = n+1, 




. .=a . .-an a, . 
ij ij zk. kj 







15.25 

6/20/62 



SCIENTIFIC 



DIMENSION A{420»ADIM) 
VECTOR VALUES ADIM = 2.0»0 
DELTA READ FORMAT NVAL,N 

A©rM(2> = H+1 

READ FORMAT INPUT, A( 1 » 1 ) . .,A( N »N+1 5 
PRINT FORMAT INVAL »N,A ( 1 , 1 ) . .. A( N ,N+1 ) 
THROUGH B, FOR K = 1,1»K .G. H 

C A(K»J) = A{K,J)/A(K,K) 

THROUGH B, FOR I = 1,1,1 .G. N 
WHENEVER I .E. K, TRANSFER TO B 
THROUGH D, FOR J = N+1,-1, J .L. K 
D A(I,J) = A(I,J)-A( I,K)*A{K,J) 

B CONTINUE 

THROUGH E, FOR I = 1,1,1 .6. N 
E PRINT FORMAT RESULT, I , A{ I ,N+1 ) 

TRANSFER TO DELTA 
INTEGER I,J,K,N 
R 

RFORMAT SPECIFICATIONS 
R 
VECTOR VALUES NVAL = $ U* S 
VECTOR VALUES INPUT = $ 6F12.4* S 
VECTOR VALUES INVAL = $7H1 INPUT//4H N = 14// 
17H MATRIX//{1H0,8F12.4)*$ 
VECTOR VALUES RESULT = $ 1H0»S20 »2HX ( » 12 ,3H ) = F12.4»$ 
END OF PROGRAM 



DATA 



t ^ 1- 1. A. -1- r*i 

0. -1. -1. -2. -9. -32* 
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15 .2 Business Data Processing Examples 

Example 1 
Problem ; Compute the social security deduction and accumulated gross pay. 
The program should read a card containing: (a) the employee's name, (b) 
his payroll number, (c) his gross pay for the current week, and (d) his 
accumulated gross pay for the current year (but not including item (c)). 
For each card read, the program should print (a) and (b) from the card, 
and, in addition, print (e) the updated gross pay, (f ) the social security 
deduction for the current week, and (g) the net pay for the current week, 
taking into account only the social security deduction. 

Analysis ; The social security deduction is currently % of the gross pay 
until the accumulated gross pay for the year exceeds $14-800.00. The 
updated gross pay can be computed from the formula, (e) = (c) + (d). The 
■social security deduction has already been made on ^dy. xnere are onUS 
three cases to consider: 

(1) (d) > il-800.00i:3n this case (f) = 0| 

(2) (d) < ij-BOO-OO and (c) + (d) > 4800.00, in this case (f ) = 3^0 of 

4800.00 - (d)^ 

(3) (c) + (d) < 4800.00, in this case f = 3^ of (c). 

The information on the cards to be read will be in the following 
format s 

Information 

(a) employee's name 

(b) pajrroll number 
(c )^ gross pay for the ciirrent week in 

the form XXX .XX 
45-52 (d) accumulated gross pay for the 

current year in the form XXXXX.XX 

The printed output will be in the following format: 



Card 


Columns 


1. 


.30 


31- 


^58 


39= 


=44 
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Line Columns 

1 

2-51 
52-3^^ 
55-^2 

46-55 

57-61 

62-64 
65-70 



Information 

Carriage control for printer 

(a) employee's name 
Blank 

(b) payroll number 
Blank 

(e) updated gross pay for current year 
in the form XXXXX .XX 

Blank 

(f ) social security deduction for current 
week in the form XX .XX 

Blank 

(g) net pay for current week in the form 
XXX .XX 



Flow Chart ; We will use the following abbreviations ^ 

KAME for employee's name (a) 

PAYHR for payroll number (b) 

CERj^SSW gross pay for current week (c) 

AQR^Sy accumulated gross pay for ciirrent year (d) 

UCffi^Y updaijed gross pay for current year (e) 

PICA social security deduction for current week (f) 

NET PAY net pay for current week (g) 
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j Start j 



READ NAME, 

PAYIR.GROSSW 

AGROSY 



(AGR0SY;ii-800. "\ 



< 



> 



FICA = 0, 



r 



< 



[ ■ GROSSW+AGROSY ; I+8OO . 



> 



FICA = .03 (i)-800. -AGROSY) 



e 



FICA = .05 GROSSW 



UGROSY=AGROSY+GROSSW 



NETPAY=GROSSW"FICA 



Print NAME, 
PAYHR/iJGROSY 
FICA^NETPAY 
on one line 
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START 



BIGGRS 

BIGTOT 
UPDATE 



READ FORMAT IN, NAME( 1 ) . . .NAME ( 5 ) ,PAYNR »GROSSW»AGROSY 

VECTOR VALUES IN = $5C6, 18 »F6.2 .F8.2*$ 

DIMENSION NAME{5) 

INTEGER PAYNR, NAME 

WHENEVER AGROSY .GE. 4800. , TRANSFER TO BIGGRS 

WHENEVER GROSSW+AGROSY .G. 4800. »TRANSF£R TO BIGTOT 

FICA = .03*GROSSW 



FICA = 0. 

TRANSFER TO UPDATE 

FICA = .03*{4800. -AGROSY) 

UGROSY = AGROSY+GROSSW 

NETPAY = GROSSW-FICA 

PRINT FORMAT OUT,NAME { 1 ) . . .NAME ( 5 ) ♦PAYNR»UGR0SY,FICA,NETPAY 

TRANSFER TO START 

VECTOR VALUES OUT = $1H0»5C6,S3» 18 ,S3»F8.2»S3 ,F5 .2»S3 ,F6.2*$ 

END OF PROGRAM 



$ DATA 
GEORGE WASHINGTON 
JOHN ADAMS 
THOMAS JEFFERSON 
JAMES MADISON 
JOHN QUINCY ADAMS 



12345678 


100. 


4800* 


12345679 


200. 


4900* 


12345680 


200. 


4600* 


12345681 


200. 


4700. 


12345682 


100. 


300. 



Alternate Program: 



START 



READ FORMAT IN,NAME ( 1 ) . . .NAMEf 5 ) ,PAYNR »GROSSW»AGROSY 

VECTOR VALUES IN = $5C6» I8»F6. 2 »F8. 2*$ 

DIMENSION NAME{5) 

INTEGER PAYNR, NAME 

WHENEVER AGROSY .GE. 4800. 

FICA = 0. 

OR WHENEVER GROSSW + AGROSY .G. 4800. 

FICA = .03*(4800.*AGROSY) 

OTHERWISE 

FICA = .03*GROSSW 

END OF CONDITIONAL 

UGROSY = AGROSY+GROSSW 

NETPAY = GROSSW-FICA 

PRINT FORMAT OUT, NAME( 1 ).. .NAME ( 5 ) ,PAYNR, UGROSY , FICA , NETPAY 

TRANSFER TO START 

VECTOR VALUES OUT = $1H0,5C6 »S3 » 18 ,S3,F8.2,S3 ♦F5.2,S3,F6.2*$ 

END OF PROGRAM 



$ DATA 
GEORGE WASHINGTON 
JOHN ADAMS 
THOMAS JEFFERSON 



12345678 
12345679 
12345680 



100. 
2 00. 
200. 



4800. 
4900. 
4600. 
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JAMES MADISON IZS'ifSSSl 200. 4700. 

JOHN QUINCY ADAMS 12345682 lOOo 300« 
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Notes on Example 1 

1. The maximum number of characters which can be stored in one 
machine word is six. Hence, we need five machine words to 
store the 50 characters allowed for the employee's name. We 

actually to be a block and that MME(5) is the last word of 
this block, In the read and print statements we specify that 
theJwhole block is to be read or printed by writing NAME(i)... 
NAME (5) and giving the format specification 5C6j i.e., 5 words 
of 6 characters . 

integer) we give an integer mode declaration stating that 
PAYNR is an integer. Similarly, since alphabetic information is 
assumed to be in the integer mode, NAME is also declared to be 
integer. 
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Example 2 



Problem ; Assiome that a master tape is available containing basic informa- 
tion for each employee; (l) tl^e employee number^ (2) his hourly rate, 
(5) gross pay to fiate, (^) amount of withholding tax withheld to date, 
(5) social security deduction withheld to date^ (6) net pay to date^ and 
(7) the number of exemptions. Input will be in the form of m cards re- 
presenting the current pay record, containing the employee's number and 
the number of hours worked during the current week. Pay is to be computed 
at time and a half for any hours worked over forty. (We shall assume 
that the input deck is already sorted according to increasing employee 
niimber J but we shall provide^ for cards which may be out of order . ) The 
last input card must have an employee number greater than the last 
employee number of the master tape . 

The withholding tax W is to be computed by the formula: 

¥ = .18 (Gross pay -= 15 n) 

where n is the number of exemptions . If n is negative, we set ¥=0 (see 
Note 2 below). The social security deduction PICA is 3 percent of gross 
pay up to $ij-800, with no deduction i'or gross pay over $4800 . 

A program is desired which will produce a listing (for each input 
card) of (a) employee number, (b) gross pay this week, (c) withholding 
tax, (d) PICA, (e) net pay for the week. Moreover, a new updated master 
tape should be prepared, with provision for saving the previous master 
tape as well. As much checking as possible should be incorporated, 
including specifying to the operator the number of the master tape needed, 
and the number to be assigned to the new tape produced by the program, and 
the automatic checking that the correct tape has been mounted on the 
unit. 
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Notes on Example 2 
Note 1: 



Note 2: 



Note 3: 



Abbreviations used here are outlined in Example 1, except 
for the following new terms: 

AWITHY accumulated withholding tax for year 

AFICAY accumulated social security deduction for year 

ANETY accumulated net pay for year 

EKEMPT number of exemptions 

In the computation of the gross pay for the current week we 
shall find it useful to be able to compute a function (which 
we shall call EXCESS.) of two numbers, say a and b, whose 
value is if a < bj and a-b if a > b. A formula for this 
function is 



EXCESS, (a, b) 



a - b + a 



where | | denote the usual "absolute value". In fact, by using 
this fixnction, a simple one-line formula for this function is: 

FICA = .03 EXCESS. ((GRjzfsSW - EXCESS . (AGR^SY, i^800.)),0) 

where AGR^SY is ass\jmed to already contain GE^SSW, i.e., to 
have been updated already. We shall also apply this function 
in the case of the withholding tax to guarantee that we do 
not make & negative deduction. Thus 

W = .18*EXCESS. (GR0SSW, 13*EXEMPT) 

To check the order of input cards (normally in order of increasing 
employee number with a large 'Employee number greater than the 
last employee number on the master tape) the program uses the 
subroutine SETE0F. (LABEL), where LABEL is the statement label 
of a statement tq be executed if -lan .end- of- file, condition is 
detected during reading. 
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Since the last input card has a large employee number, 
the first end-of-file condition is normally detected at the 
end of processing, "but an illegal input card may also exist 
with a high employee number. After the first end-of-file 
is detected the end-of-file return is changed and the input 
tape checked for end-of-file. If no end-of-file exists a 
comment is printed to change tapes and processing begins again. 



15.35 
6/20/62 



o U4- 

o\ 
ro 



[ Start y 



Read 

TAPEWO 



Print 
TAPEWO 
for oper- 
ator 




a 



IF 



Pause while 
operator 
hangs master 
tape 




Initialize 

grand 
totals 



Set End 
of File 
to 8 





Read Input 
€ard 



PAYto, HOURS 



Copy- 
Master 
Tape 



READ 
MASTER 
RECORD 
WUMB.^ 

ETC, 



e- 



< 



^IMB'.B 



AYWR 



1^ 




Compute !^ 
Accumulate 
GROSSW 



> 



Compute 

withhold- 
ing tax 

W(^ PICA 



WETPAY = 
GROSSW -W 
-PICA 



Accumulate 
NETPAY,Wj 
PICA \ 
check tot- 
als 



Ki) 



PRINT 
EimOR 
COMMENT 






CO 

o 
o 

CO 
CO 
H 
SSI 
Q 




Accumiilate 

^Jrand 

Totals 



Print 
Man's 
Totals 



Copy 

lyfester 
Record 



HD 



VD 

ro o 
• OJ 

H VO 



H 
CQ 
CD 

O 

o 

l1 

m 



B 




Set End 
of File 

to 6 




Print 
Comment to 

Operator 



Read Input 
Card. 



Print 
Comment 
to 
Operator 



Pause 
while oper- 
ator re- 
selects 
tape units 



Backspace 

Input 

Tape 




Pause while 
operator 
removes 
master tapes 



Print 
Grand 
Totals 



e 



BUSINESS DATA FROCTSSING 



START 



M FILE 



READ FORMAT IDENT,TAPENO 
INTEGER TAPENO,PAYNR,NUMB,J,OLTAPE 
PRINT ON LINE FORMAT OPER,TAPENO 
PAUSE NO. 1 





REWIND TAPE 4 


TEST 


REWIND TAPE 3 




READ BINARY TAPE 3»0LTAPE 




WHENEVER OLTAPE .E. TAPENO, TRANSFER TO MAIN 




PRINT ON LINE FORMAT WRONG 




PAUSE NO. 3 




TRANSFER TO TEST 


MAIN 


CUMGRS = 0. 




CUMFIC = 0. 




CUMNET = 0. 




CUMW = 0. 


REDO 


EXECUTE SETEOF.{M FILE) 




WRITE BINARY TAPE 4,TAPEN0+1 


READ(l) 


READ FORMAT EMPLOY jPAYNR , HOURS 


READ(2) 


READ BINARY TAPE 3» NUMB , RATE, AGR0SY»AWITHY, 



'0) 



1AFICAY,ANETY, EXEMPT 
WHENEVER NUMB.E. PAYNR 

GROSSW = RATE*HOURS+.5*RATE»EXCESS.(HOURS,40.) 

AGROSY = AGROSY+GROSSW 

W = .18*EXCES5. (GROSSW, 13. *EXEMPT) 

FICA = .03*EXCESS. ( (GROSSW-EXCESS.fAGROSY»4800.)) 

NETPAY = GROSSW-W-FICA 

AWITHY = AWITHY+W 

AFICAY = AFICAY+FICA 

ANETY = ANETY+NETPAY 

CUMGRS = CUMGRS+GROSSW 

CUMFIC = CUMFIC+FICA 

CUMNET = CUMNET+NETPAY 

CUMW = CUMW+W 

WHENEVER .ABS. ( AGROSY-ANETY-AFlCAY-AWITHY ) 
IFORMAT ERROR, PAYNR 

PRINT FORMAT OUTPUT, PAYNR , GROSSW, W, FICA, NETPAY 

J = 1 

OR WHENEVER NUMB. G. PAYNR 

PRINT FORMAT ORDER, PAYNR 

BACKSPACE RECORD OF TAPE 3 

TRANSFER TO READ(l) 

OTHERWISE 

J = 2 

END OF CONDITIONAL 

WRITE BINARY TAPE 4, NUMB, RATE, AGROSY, AWITHY, AFICAV, ANETY, 
lEXEMPT •• t- T , 

TRANSFER TO READ{J) 
END OF FILE TAPE 4 
REWIND TAPE 3 



GE. .005, PRINT 
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REWIND TAPE 4 

EXECUTE SETEOF.iC FILE) 

READ FORMAT EMPLOY » DUMMY » DUMMY 

PRINT FORMAT NOMAN jPAYNR »TAPENO 

PRINT ON LINE FORMAT NOMAN ,PAYNR »TAPENO 

PAUSE NO. 4 

TAPEN0=TAPEN0+1 

BACKSPACE RECORD OF TAPE 7 

READ BINARY TAPE 3»DUMMY 

TRANSFER TO REDO 
C FILE PRINT ON LINE FORMAT OFF ,TAPENO »TAPEN0+1 

PAUSE NO. 2 

PRINT FORMAT TOTALS9CUMGRS »CUMFIC »CUMNET ♦CUMW 

EXECUTE SYSTEM. 
R 

INTERNAL FUNCTION EXCESS. { X >Y) = ( X-Y+ .ABS. (X-'Y))/2. 
R 

RFORMAT SPECIFICATIONS 
R 

VECTOR VALUES IDENT = $I8*$ 

VECTOR VALUES OPER = $15H4M0UNT TAPE NO. I8,S2,30HON TAPE UNI 
IT NO. 3»PRESS START*$ 

VECTOR VALUES WRONG = $48H4THE WRONG TAPE HAS BEEN USED. PLEA 
ISE TRY AGAIN. *$ 

VECTOR VALUES EMPLOY = $I8>F10.2*S 

VECTOR VALUES ERROR = $37H0ERROR IN CHECKING TOTALS FOR MAN N 
10. I8*$ 

VECTOR VALUES OUTPUT = $1H0 s 18 »4F20 . 2*$ 

VECTOR VALUES OFF = $24H4REM0VE TAPE 3»LABEL IT 18 »S4 »23HREM0 
IVE TAPE 4» LABEL IT I8*$ 

VECTOR VALUES NOMAN = $38H0THERE IS NO MASTER RECORD FOR MAN 
1NO.I8/22H0PULL TAPE 3. LABEL IT I 8/51H0RESELECT TAPE 4 AS TAP 
2E 3 AND HANG BLANK TAPE ON 4/16H0THEN PUSH START*$ 

VECTOR VALUES ORDER = $8H0MAN NO.I8»43H IS OUT OF ORDER OR NO 
1 MASTER RECORD EXISTS*$ 

VECTOR VALUES TOTALS = $13H1CUM. GROSS =F10.2/12H0CUM. FICA = 
1F10.2/11H0CUM. NET = F10.2/23H0CUM. WITHHOLDING TAX = F10.2*J 

END OF PROGRAM 
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Example 3 



Protilem ; Mortgage Payment. The type of mortgage we consider here Is the 
fixed principal type for which each installment consists of an interest 
payment, a fixed amount to be deducted from the outstanding principal, 
and an additional amount to be placed in escrow, to be used to make 

Assume that a master card file is available containing the following 
information for each mortgage: (l) the mortgage number; (2) amount of 
outstanding principal; (5) annual payment on principal; (h) interest 
rate; (5) annual escrow payment; and (6) current escrow balance. There 
is also a file of cards available containing the current payment record 
consisting of mortgage number and amount of pajonent received. The master 
file and current payment file are assumed to be in order of increasing 
mortgage number. 

The program is to read a card from the current payment record and 
check to see if it is acceptable. A payment is deemed acceptable if 
it consists of a single normal payment (i.e., a payment consisting of 
a single principal pasraient, a single escrow payment, and an interest 
payment for a single period) or if it consists of exactly two normal 
payments and any number (i=0,l,2...) of principal payments. 

Note Is Observe the use of the alphabetic constant 

to print the character "$" using a C-field specification (see . 
section 2.1.3, Alphabetic Constants). 
Note 2: The current payments are processed until the file is exhaxisted. 
The detection of the end-of-file on reading transfers control 
to the section of the program which punches the new master file. 



15.58 
6/20/62 



( start V 



v^ 



READ 
N 



/THROUGH A, 
FOR 
i=l,l,i>N 



^Read Master 
File RECORD 

RECORD (1,6) 




A 




a 



OJ 

fO O 

• OJ 

H VD 



H 
CQ 
CQ 

O 
O 

■a 

m 

CQ 




a 



READ Current 

File 

IDEET^AMOUNT 




THROUGH B, 
FOR 
1=1,1 1>N 



/F^ 



M 




;ortgage:idf™^> 



DEm 



Connect 

OUT OF 

order 




Connect 
WO 

MASTER 
RECORI 




s 

pq 




Compute 

Single 

Payment 

due:) 



-Q 



AMOUNT :DUE1 



;^ 



Compute 
Double Pay- 
ment DUE2 



Update Out- 
standing 
Principal 
and Amount 
In Escrow 




O -f^ 

ro 




THROUGH C^ \ 

FOR PAY=^,HMJAD. 

PAYMENT, 

ANNUAL PAY- 
mh:nt 




Print 
Unsatis- 
factory 
Payment 







©v 



Update 

Outstanding 

Principal 



Update 
Amount in 
Escrow 



Set indica- 
tor for 
punching new 
master recorc 




B 



CO 



O 

o 

CO 
M 

B 



On end of file 



'THROUGH D, 
FOR 
i=l,l, i>N 



<; 



Indicator i 



i^ 



> 



< 



Punch new 

master 

cards 




BUSINESS DATA EROCESSING 



DIMENSION RECORD( 1400»DIM) 

INTEGER I»NUMB 

VECTOR VALUES DIM = 2,1»7 

VECTOR VALUES DOLLAR = -$=$ 
START READ FORMAT SIZE»NUMB 

THROUGH A» FOR I = 1»1»I .G. NUMB 
A READ FORMAT MASTER »RECORD ( I , 1 ) . . .RECORD ( I ,6 ) 

EXECUTE SETEOF. (UPDATE) 

I = 1 
CARDS READ FORMAT PAYMT, IDENT»AMOUNT 

THROUGH B* FOR I = I -. 1 » I .G. NUMB 

WHENEVER REC0RD(I,1) .E. IDENT 

DUEl = RECORD( I ,3)+REC0RD(I»5)+REC0RD{ I»4)*REC0RD( I»2) 

WHENEVER ,ABS. ( AMOUNT-DUEl ) .L. .005 

RECORD (I»2) = RECORD{ I>2}-RECORD(I,3) 

RECORD ( I»6)=REC0RD(I ,6)+REC0RD{I »5) 

TRANSFER TO CODE 

OTHERWISE 

DUE2 = 2.*DUE1 - RECORD ( I »4 ) *RECORD ( I » 3 ) 

END OF CONDITIONAL 

WHENEVER .ABS. ( AM0UNT-DUE2 ) .L. .005 

r\L-<^wKL>v i 9 ^ } -rM-v_wn.u \ j. 9^ ; £.«^Ki-\-wr\L/i i. Slj t 

TRANSFER TO ESCROW 
OR WHENEVER AMOUNT .G. DUE2 

THROUGH C, FOR PAY = RECORD( I , 3 ) »RECORD (I » 3 ) » 
lAMOUNT .L. DUE2+PAY 
C WHENEVER .ABS. ( AM0UNT-DUE2-PAY ) .L. .005* TRANSFER TO 

IPAID 
TRANSFER TO OVRPAY 
PAID RECORD { l92)=RECORD{ I , 2 ) -2 .^RECORD { I ♦3) -PAY 
ESCROW RECORD ( I ,6)=REC0RD( I , 6 ) +2 .*RECORD{ I . 5 ) 
CODE RECORDC I 97)=1. 

OTHERWISE 
OVRPAY PRINT FORMAT REJECT » IDENT »DOLLAR »AMOUNT 
END OF CONDITIONAL 

OR WHENEVER REC0RD(I,1) .G. IDENT 
PRINT FORMAT ORDER « IDENT sDOLLAR »AM0UNT 
OTHERWISE 
B CONTINUE 

1 = 1 

PRINT FORMAT NONE, IDENT 
END OF CONDITIONAL 
TRANSFER TO CARDS 
UPDATE THROUGH D, FOR I=1»1,I .G. NUMB 

D WHENEVER RECORD ( I , 7 ) .G. 0., PUNCH FORMAT MASTER , RECORD ( I , 1 ) 

1REC0RD{I»6) 
R 

RFORMAT SPECIFICATIONS 
R 
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VECTOR VALUES SIZE=$I10*$ 

VECTOR VALUES MASTER=$F10. »5F 10. 2*$ 

VECTOR VALUES PAYMT = SFIO.O » F 10 . 2*S 

VECTOR VALUES REJECT = $20H0PAYMENT ON MORTGA6E»F10. ,3H, »C 
11»F10.2»19H IS UNSATISFACTORY, *$ 

VECTOR VALUES ORDER = $26H0PAYMENT CARD FOR MORTGAGE »F10.0 » 3H 
Is ClsF10«2»4^H IS OUT OF ORDER OR NO MASTER RECORD EXISTS. *$ 
"vector values none = $41H0NO MASTER RECORD EXISTS FOR MORTGAG 
IE NO., FIG. 0*$ 

END OF PROGRAM 
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Example k 
Problem : Computation of actuarial commutation colimins based on an 
arbitrary set of mortality rates and an interest rate, as an external 
fimction to be used by another program. 

Analysis ; Commutation columns , which are very important tools in 
actuarial problems are generated very easily by means of the formulas 
given below. The quantities M , N j and D in these formulas occur most 
often in combination, as in the computation of P . Assuming a population 
of some initial size (at x = bQ) (hsre 1,000,000),^ is the number living 
at age x {so that ^ ^ = 1,000,000), q is the mortality rate , and d is 
the number of deaths at age x . Thus d = Q' Jc . The quantity D is 

- X "^ X X 

computed by the formula D , =x (l+i) , where i is the interest rate. 

XX f 1 ) 

Another quantity, C is given by the formula C = d (l+i)"^'^ . It „ 

X X X o 

can be used, for example, to compute the cost of term insurance, since «- 

X 

is the premium for one year term Insurance of $1 at age x. 
The sums M and W are obtained by the formulas 

XX 

CO CD 

M=ZG,n=ZD 
X y X y 

y=x "^ y=x ^ 

We note that for some w, we al'Ways have q^ = 1, so that /t = 

(since Jc , = i -d = ^ - jl =0), therefore D , = 0, d , = 0, C -, = 0, 

^ w+1 '^w w '^w w " w+1 ' w+1 ' w+1 ' 

and the sums for M and W are actually finite sums . 

The three most useful quantities computed here are (l);, P = M^/N , 

.X ..-^ X 

which is the annual premium payable for an entire life for $1 of whole 
life insurance, (2) A = M /d ^ which is the single: premium payable. at 
age X for $1 of whole life insurance, and (3) A = N /D, which- is the 

X - X X 

present value at age k of a whole life annuity of $1, first payment at 
age X. 

Printing of results is under control of an input variable PRINT. 
Certain relationships must hold between some independently computed 
values, and these are used as checks on the computations 

^0 ^0- %'^^ \ 
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P^ - 1/a^ - i/(l+i) 



These cannot be expected to come out exactly equal, because of round-off, 
but they should differ by very little. 
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(Entry Y 




•b =1,000,0001 



D =ji -v 
X x 



-b 



k. 



THROUGH A, 
FOR 



\ x=b j\, 
\ x > w 



/ 



V = v/(l+i) 



C = d -v 
x x 



d =q^.i 



X 



X 



X+1 X X 






vo 

^ o 
• t\j 



CO 

w 
o 
o 

m 
m 



H 





N = D 
w w 



M = C 
w w 



THROUGH B, 
FOR 

X=W-l,"lj 



X < b. 




-y y__y"M^N Check ;1 \_^ 



> 



Error 
Comment 



THROUGH G, 
FOR 

x-b ^ 

X > w 



X 

x+1 X 



M = 

X 

M ,+C 
x+1 X 
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R 

RSAMPLE CALLING PROGRAM 

R 

READ FORMAT IN, Q{ ) . . .Q{ 100 ) 

VECTOR VALUES IN =$12F6.5* $ 

DIMENSION Q(120) ,L ( 120 ) tSMALLD ( 120) fBlGD ( 120 ) .C ( 120 ) ,N(120)» 

EXECUTE COMFCN.{Q,0»099».03.L» SMALLD . BIGOtC »N sM.BlGA »SMALLA » 
IP.l) 
INTEGER PRINT 
END OF PROGRAM 

The External Fiinctlons 

RCOMMUTATION TABLE FUNCTION 
RIF PRINT = 0»SUPPRESS PRINTING 

EXTERNAL FUNCTION{Q,BZERO,OMEGA,I ♦L»SMALL0,8IGD,C.N, 
1M»BIGA»SMALLA»P,PRINT) 

ENTRY TO COMFCN. 

INTEGER BZEROjOMEGA, PRINT, X 

LfBZERO) = 1E6 

V = (i. + n .p. -ezERo 

THROUGH A»FOR X = BZER0.1,X .G. OMEGA 
SMALLDCXl = Q(X)*L(X) 
BIGDCX) = L(X)*V 

V = V/(le+I ) 

CfX) = SMALLD{X)*V 

L(X+1) = L(X)-SMALLD{X) 

N{ OMEGA) = BIGDC OMEGA) 

M( OMEGA) = C{ OMEGA) 

THROUGH B» FOR X = OMEGA-1,-1, X .L. BZERO 

N(X) = N{X+1 ) + BIGD(X) 

M(X; = M(X+1 ) + CfX) 

WHENEVER .ABS. ( M{ BZERO )+N ( BZERO+1 )-N {BZERO )/( I+l .) ) .G. I.9 
ITRANSFER TO MNERR 

THROUGH G ,FOR X = BZERO, 1,X .G. OMEGA 

BIGA(X) = M(X)/BIGD(X) 

SMALLA(X) = N(X)/BIGD(X) 

PCX) = M{X)/N(X) 

WHENEVER .ABS. ( P< X )-l . /SMALLA ( X ) + I/{I+1.)) .6. 
llE-4, TRANSFER TO PERROR 

CONTINUE 

WHENEVER PRINT .E. 0, FUNCTION RETURN 
R 

ROUTPUT GENERATOR 
R 

PRINT FORMAT HEAD0I9I 

VECTOR VALUES HEADOl = $1H1,4HI = F5.4// 
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1 4H X»S13»4H0(X) »S18»4HL(X) ,S15,10HSMALL D{XJ» 

2 S8»1HX*$ 

THROUGH BETA, FOR X = BZER0,1,X .G. OMEGA 
BETA PRINT FORMAT Fl ,X ,0 ( X ) » L { X ) ,SMALLD{ X ) ,X 

VECTOR VALUES Fl = $1H0 , 13 ,3E22 . 9 , 1 7*$ 

PRINT FORMAT HEAD02 

VECTOR VALUES HEAD02 = $1H1,4H X,S11,8HBIG D{X), 
1 S16,4HC{X) »S18,4HM{X) ,S18,4HN{X) ♦S11,1HX*$ 

THROUGH GAMMA.FOR X=BZER0,1,X .G. OMEGA 
GAMMA PRINT FORMAT F2 »X» BIGD ( X ) ,C ( X ) ,M ( X ) ,N ( X ) ,X 

VECTOR VALUES F2 = $1H0, I494E22 . 9 , I 7*$ 

VECTOR VALUES HEAD03 = $4H1 X,S11,8HBIG A{X)» 
1S13»10HSMALL A(X), SI 5 ,4HP ( X ) ,S1 1 , IHX *$ 

THROUGH DELTA, FOR X = BZER0,1,X .G. OMEGA 
DELTA PRINT FORMAT F3 »X , BIGA ( X ) ,SMALLA ( X ) ,P ( X ) ,X 

VECTOR VALUES F3 = $1H0, I 3 ,3E22 . 9 , 1 7*$ 

FUNCTION RETURN 
PERROR PRINT FORMAT PERR , P ( X ) ,SMALLA ( X ) , I 

VECTOR VALUES PERR = $27H0ERROR ON P CHECK. P(X) = E18.9, 
1S10,13HSMALL A(X) = E18.9,S10,4HI = F5.4*$ 

TRANSFER TO G 
MNERR PRINT FORMAT MNERRl 

VECTOR VALUES MNERRl = $19H0ERROR ON M,N CHECK*$ 

TRANSFER TO E 

END OF FUNCTION 

$ DATA 

2258 577 414 338 299 276 261 247 231 212 197 191 

192 198 207 215 219 225 230 237 243 251 259 268 

277 288 299 311 325 340 356 373 392 412 435 459 

486 515 546 581 618 659 703 751 804 861 923 991 

1064 1145 1232 1327 1430 1543 1665 1798 1943 2100 2271 2457 

2659 2878 3118 3376 3658 3964 4296 4656 5046 5470 5930 6427 

6966 7550 8181 8864 9602 10399 11259 12186 13185 14260 15416 16657 

17988 19413 20937 22563 24300 26144 28099 30173 32364 34666 37100 39621 

44719 54826 72467100000 
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15.5 Symbol Manipulation and Reciirsive Function Examples 

Example 1 
Problems Find the first occurrence of an arbitrary word in a given text. 

* -» ■ ^ -e- I -nTT ji T. _J!» _V *- -ff J_1__ J_ L. 3 m /t \ 

T(n) be the text stored one character per word. Let L be the number of 
letters in the word which is stored one character per word in W(1)...w(L), 



1 c I.Q 
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{ Start ) 



Read N, T(i; 

...t(n) 



f Read L, W(l) 

... w(l) 



/ THROUGH 

/scan, for 

\I =1,1,1 > 
\ N-L+1 / 



THROUGH 
TST, FOR 

J = 0,1, 
J > L 



(t(1+J);¥(J+1)\. 



( TST j 



PRINT w(i; 
...W(L) 
"f oxmd at 
character", 

I 





PRINT 
"Word not 
Found" 
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ALPHA 

TST 
SCAN 



DIMENSION T{720)»W{30) 

NORMAL MODE IS INTEGER 

READ FORMAT 

READ FORMAT 

READ FORMAT 

READ FORMAT 

TWDniif;H £;rAN 



CNT»N 

TXT»T{ 1)...T(N) 

CNT»L 

TXT»W{ 1)...W{L) 

. FOR T=1.1.T -R- N-i +1 



THROUGH TST» FOR J=0.1»J .GE. L 

WHENEVER T(I+J) .NE. W(J+1)» TRANSFER TO SCAN 

PRINT FORMAT OUT, I »W( 1 ) . . .W( L ) 

TRANSFER TO ALPHA 

CONTINUE 

PRINT FORMAT NOT 

TRANSFER TO ALPHA 

VECTOR VALUES CNT=$I3*$ 

VECTOR VALUES TXT=$72C1*$ 

VECTOR VALUES OUT=S11HOCHARACTER I3,13H 

VECTOR VALUES NOT5=$15HOWORD NOT FOUND*$ 

END OF PROGRAM 



IS START OF 30C1«$ 
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Example 2 

Problem : Evaluate the recursive function, 

f (0) = 1 

f (n) = f (n-1) * n 

Analysis ; This is the deflaition of nl. Although nl can be evaluated 
directly using a THRjZiUGH statement, in this example it will be evaluated 
using its recursive definition to illustrate how recursive functions can 
be handled in MAD. 
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EXTERNAL FUNCTION (N) 

NORMAL MODE IS INTEGER 

ENTRY TO FACT. 

WHENEVER N .E. 0. FUNCTION RETURN 1 

SAVE RETURN 

SAVE DATA N 

TI = FACT.(N-l) 

RESTORE DATA N 

RESTORE RETURN 

FUNCTION RETURN T1*N 

END OF FUNCTION 



In order to use this f 'unction the calling program would have to specify 
a list for use in the SAVE and RESTj^isRE statements . The following is 
an example of a program which uses FACT.. 

DIMENSION LIST (100) 
NORMAL MODE IS INTEGER 
SET LIST TO LIST 
BACK READ FORMAT IN, NR 

PRINT FORMAT OUT, NR, FACT.(NR) 

TRANSFER TO 3ACK 

VECTOR VALUES IN = $I2*$ 

VECTOR VALUES OUT = $4H0N= I3,14HN FACT0RIAL= 111*$ 

END OF PROGRAM 
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Example 3 

Problem ; To find the greatest common divisor of two integers Z and Y. 
Analysis ; The greatest common divisor is defined recursively by three 
equations : 

I Y > Z, -»GCD.(y,Z) 
GCD.(Z,Y) = \ REM.(Z,Y) = ^Y 

\^ otherwise -> GCD. (REM. (Z,Y),y) 

where REM.(A,B) is the remainder of A/B. This function expects the 
arguments to be found on the temporary storage list as the two most recent 
additions. The use of the list as a parameter list makes the establishment 
of dummy variables unnecessary. This is less efficient than the usual 
way of defining functions but serves to remove many pitfalls encountered 
in using dummy variables with recursive functions. 
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EXTERNAL FUNCTION 

INTERNAL FUNCTION REM,(A,B)= A - {A/B)*B 

ENTRY TO GCD. 

NORMAL MODE IS INTEGER 

RESTORE DATA Z,Y 

WHENEVER Y .G. Z 

SAVE RETURN 

SAVE DATA Z,Y 

X = GCD.(O) 

RESTORE RETURN 

FUNCTION RETURN X 

OR WHENEVER REM. (Z»Y) .E. 

FUNCTION RETURN Y 

END OF CONDITIONAL 

SAVE RETURN 

SAVE DATA REM* (Z,Y5 »Y 

X = GCD. (0) 

RESTORE RETURN 

FUNCTION RETURN X 

END OF FUNCTION 



When called upon for a value, a functioasucli as GCD. laust; have-.at 
least one argument (in this example a dummy argument of zero Is used) even 
^tiough the ai^umenti is never called upon. ; -SEhis is. because .GCiD. :ls;:.thej 
name of the function while GCD. .(..•) is the value of the function. 

The SET LIST Tp statement need be executed once, either in the main 
program or in a subprogram (but before any use of SAVE or RESTjZ5EE), since 
the SAVE and RESTORE statements always refer to the current list. 



NORMAL MODE IS INTEGER 

SET LIST TO LIST 

DIMENSION LIST (50) 

READ FORMAT IN,M,N 

SAVE DATA M,N 

PRINT FORMAT OUT »M sN »GCD . i ) 

TRANSFER TO S 

VECTOR VALUES IN = $216*$ 

VECTOR VALUES OUT = $1H0.3HM= »I7,S10»3HN= ♦ 1 7 ,S10 . 5HGCD = 17 

END OF PROGRAM 
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SYMBOL MNIPUIATIOW AND RECURSIVE FUNCTION 

Example h 

Problem ; To evaluate Tschebychev polynomials . 

Analysis : The Tschebychev polynomial T(N,X) is defined recursively as 

follows : 

N = -^1 
T(W,X) = ^ N = 1 ^X 

N > 1 -*2*X^(N-1,X)-T(N-2,X) 

It is important to understand that when an expression is written as 
an argument of a function its value is computed and stored in a temporary 
location. It is this location (or address) which is acti;ially used as the 
argument of the function. The implication of this use of a temporary 
location is that often expressions cannot be used as arguments of re- 
cursive functions . 
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EXTERNAL FUNCTIOINJ (N»X) 

ENTRY TO TSCHEB. 

INTEGER N,Z 

WHENEVER N .E. 0» FUNCTION RETURN 1. 

WHENEVER N .E. 1» FUNCTION RETURN X 

SAVE RETURN 

CAWF nATfi M-5 

Z = N-1 

Y = 2**X*TSCHEB.{Z»X) 
RESTORE DATA Z 
SAVE DATA Y 
M=TSCHEB.{Z»X! 
RESTORE DATA Y 
RESTORE RETURN 
FUNCTION RETURN Y-M 
END OF FUNCTION 

A Prograin vhi--h uses TSCHEB. is? 

SET LIST TO LIST 
DIMENSION LISTdOOO) 
BEGIN READ FORMAT INPUT. N»X 

PRINT FORMAT OUTPUT .N » X ♦TSCHEB. ( N »X ) 
TRANSFER TO BEGIN 
RFORMATS 
VECTOR VALUES INPUT=$I 5»F10»2*$ 

VECTOR VALUES 0UTPUT=$1H0»4HN= »I6»4H X« F10*2» 
111HOFUNCTION= F15.6*$ 
END OF PROGRAM 
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SUMMARY OF MAD STATE^ISNTS 

20 Summary of MAD Statements 

I. Declarations 

A. Remark 

R in column 11 

any remark in columns 12-72 

B . Mode 

a. f^ V.f w^j •••i U^ 

where 7)| is one of; 
FL0ATIIG V^im 
INTEGER 
B0^LEAW 
P'UNCTI|Z^N WAME 
STATEMENT LABEL 

and each U. is a variable name or function name. 

i 

b. NiZiRMAL M0DE IS "fYl 

where y/| is one of five listed in a . ebove . 

C. EQUrv\ftLENCE (V-, , V,, .-,. V ), (V , , V _, . . . , V ),-.., 

^ 1^ 2' ' n'-* ^ n+1^ n+2^ ^ n+v ' 

(V ^o V ^^ .... V ) 
n+p+q+l"* n+p+q+2^ ^ n+p+q+r' 

where each V. is a variable name or linearly subscripted variable, 

D . PROGRAM CpMMpN V, , V^ ^ . . . ^ V 

where ea 3h V , is a variable name . 

i 



E . ERASABLE V, , "/^ , . . . , J 



x\ 



where each V. is a siaaple variable or array. 
F. DIMENSION ^ fa,), Vglc^), '--.^Cc^) 

where each V. is a variable name and each a, is 5 

a. for vector DIMEJISII^N^ot^ is one argument, an integer constant 
which is the largest value that the subscript of V. will assumej . 

b. for matrix DIMENSION; of. represents two arguments; the first 
is an integer constant, which is the largest value that the 
subscript of V. will assume ^ the second is the name of the first 



element of a dimension vector for V. ■. 

1 
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G. VEGTjZjR VALUES 1/"= c„^ c 



QJ* -j_.! 



• 3 



n 



or 



VECTjZjR VALUES tr= $GqC, ..., c^$, ..., .$c^^^^^, 
The two above types of values way be intermixed , 
H. INTERKAL FUHCTI^N F.(A^, k^, ..., A^) = E 

(single statement internal function definition) 



n+p+q.^ 



II".» Executable Statements 

A. Substitution 

V = F 

B. TRANSFER T^ ^ 

C . Conditional 

a- Simple Conditionals 

WHENEVER B,Q 
b. Compound Conditional j 
J . WHEHEU'ER B 



i 



1 

a • • 

... \ 

tSR WHEKEVER B, 



a 



*<0. 0E WHEMEv'ER B,, 



£i, ^ EKD 0F C^ITOITIjJJnAL 



k+1 



th 



* The k statement may be replaced bys 

JB. WHEJ3EVSR IB 

ors 
jZ^THERWlSE 
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D . C^MTINUE 

E. THRJZ5UGH (iteration) 

a. THRjZiUGH J, F0E VALUES 0F V = E^, E^^ • • '* \ 

b. THRj/iUGH J, Fj6r V = E^, E^, B 

F. PAUSE N^ n 

G. EXECUTE ,C.(A^, A , ...,A^) 

or 
EXECUTE C . 
H. mD 0F ER0GRAM 



20.2 
6/20/62 



